CN108595772A - A kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor - Google Patents
A kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor Download PDFInfo
- Publication number
- CN108595772A CN108595772A CN201810263411.8A CN201810263411A CN108595772A CN 108595772 A CN108595772 A CN 108595772A CN 201810263411 A CN201810263411 A CN 201810263411A CN 108595772 A CN108595772 A CN 108595772A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- rotor
- permanent magnet
- stator
- motor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Permanent Field Magnets Of Synchronous Machinery (AREA)
- Linear Motors (AREA)
Abstract
A kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor of the present invention, include the division to motor complexity field region and rule field region, and equivalent to two kinds of regions progress diamond shape small grid subdivisions and magnetic circuit magnetic conductance respectively；Magnetic conductance node connection relation at two class zone boundaries is improved, complete 2D magnetic circuits subdivision model is built, establishes magnetic conductance solution matrix equation successively；Build the relational expression of magnetic conductivity convergence factor and motor teeth portion flux density maximum difference；Solve magnetic conductance matrix equation, realize that the rapid solving of matrix obtains the magnetic flux and magnetic potential by each magnetic conductance node using iterative algorithm, and then calculate the magnetic flux density and magnetic conductivity of each magnetic conductance, and it calculates and constrains according to motor electromagnetic, the electromagnetic parameters such as each phase magnetic linkage of motor, winding counter electromotive force can be calculated, are finally compared with the result of finite element.The present invention is directed to the double permanent magnetism vernier motors of rotor and carries out 2D magnetic circuit subdivisions for the first time, and the scheme provided can provide for the double magnetoes of the type with reference to research.
Description
Technical field
The present invention relates to a kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor, belong to Electromagnetic Calculation
Field.
Background technology
The double permanent magnetism vernier motors of rotor, to the modulating action of space magnetic field, have in low-speed run state using modulation tooth
Under the characteristics of capable of exporting big torque, be suitable for agricultural mechanical equipment field.The double permanent magnetism vernier motors of rotor are more traditional forever
Magnetic vernier motor improves torque density, while not only reducing motor cost using mixing magnetic material, also reduces rotor
Part heat demagnetization influences.Traditional permanent magnetism vernier motor is divided into two kinds of stator type and rotor-type by permanent magnetism body position, the permanent magnetism vernier
Motor is the double permanent-magnet types of rotor, i.e., all there is permanent magnet on the stator of motor and on rotor.In the case of same volume,
The torque density of motor is improved, and also improves the space availability ratio inside motor.
It is at present business finite element (Finite Element to the common analysis method of the double permanent magnetism vernier motors of rotor
Analysis, FEA) software, since its high-precision is widely used with easy to operate in design of electrical motor each stage.But FEA software bases
The whole mesh subdivision principle inside the motor, time-consuming that efficiency is low for analysis, is not suitable for the preliminary design phase of motor.And 2D magnetic circuits
Subdivision modeling method has then merged traditional Magnetic Circuit Method and mesh generation method, and Computing Principle is simple, and precision is high, and emulates time-consuming short
Bring higher design efficiency.By carrying out the modeling of 2D magnetic circuit subdivisions to the double permanent magnetism vernier motors of rotor, can solve rapidly
Go out the electromagnetic parameters such as motor magnetic linkage, back-emf and torque, realizes the high efficient design at motor research initial stage.
Invention content
The purpose of the present invention is to provide the analyses of magnetic circuit model in a kind of double permanent magnetism vernier motor complexity magnetic fields of rotor
With method for solving, it is main including the mesh generation of complicated field region magnetic conductance, the approximate magnetic conductance in rule magnetic field in motor substitute with
The foundation of magnetic conductance matrix model.Model is added in the constraint formula that convergence factor is changed, and the Fast Convergent for matrix solves.
To achieve the above object, the technical solution adopted by the present invention is：A kind of 2D magnetic of the double permanent magnetism vernier motors of rotor
Road subdivision modeling method, includes the following steps：
Step 1, motor complexity field region and rule field region are divided using finite element software；
Step 2, the diamond shape small grid subdivision of motor complexity field region is established；
Step 3, the magnetic circuit magnetic conductance for establishing motor rule field region is equivalent；
Step 4, magnetic conductance node connection relation at two class zone boundaries is improved, builds complete 2D magnetic circuits subdivision model, successively
Establish magnetic conductance solution matrix equation；
Step 5, the relational expression of magnetic conductivity convergence factor and motor teeth portion flux density maximum difference is built；
Step 6, magnetic conductance matrix equation is solved, realizes that the rapid solving of matrix is obtained through each magnetic conductance section using iterative algorithm
The magnetic flux and magnetic potential of point, and then calculate the magnetic flux density and magnetic conductivity of each magnetic conductance；
Step 7, it according to the magnetic circuit parameters solved in step 6, calculates and constrains according to motor electromagnetic, each phase of motor can be calculated
The electromagnetic parameters such as magnetic linkage, winding counter electromotive force.
Further, the double permanent magnetism vernier motors of the rotor are the three phase electric machine of the pole of 12 slots/28, are divided into stator, air gap, turn
Son and four part of shaft；Include stator yoke, stator teeth, stator slot, armature winding and pole shoe permanent magnet array, electricity in stator
Pivot flute profile is flat-bottom slot, and for armature winding using centralized winding method, span is 4 stator slots；Stator permanent magnet material is
NdFe35 has 2 notches in each pole shoe of stator, and there are one Halbach permanent magnet arrays, each Halbach permanent magnetism battle arrays for each dress
The secondary permanent magnet that row rush the main permanent magnet of magnetic by intermediate radial and both sides tangentially rush magnetic is constituted, and can be effectively improved point of air-gap field
Cloth situation；Rotor is cylindrical shape, and the installation embedding permanent magnet of table of slotting on surface, permanent magnet material is ferrite Y30, the embedding permanent magnetism of table
Body cross section is trapezoidal, is evenly distributed on rotor circumference direction；The material of stator core and rotor core is silicon steel sheet
DW540_50；Air gap is between stator and rotor, air gap thickness 0.5mm；Machine shaft is made of un-conducted magnetic material, is real
The heart is cylindric, and is connect with rotor coaxial.
Further, in the step 1, according to the simulation result of finite element software, motor complexity Distribution of Magnetic Field region is main
It concentrates in a segment limit of rotor two-sided permanent magnet body, including at air gap, stator pole shoes and rotor teeth portion region；Motor rule
Distribution of Magnetic Field region includes stator yoke and teeth portion, rotor yoke region, since both sides permanent magnet magnetizing direction has been fixed, forever
The distribution of magnet internal magnetic field is also that rule is distributed.
Further, in the step 2, the double permanent magnetism vernier motors of rotor by Halbach permanent magnet arrays and single rush magnetic
Direction set of permanent magnets is distributed at bilateral permanent magnet, and structure is complex, thus leads to magnetic in a segment limit of rotor and air gap
Field distribution is complicated.Therefore to entire complicated Distribution of Magnetic Field region progress diamond shape small grid subdivision, each small grid region equivalent at
One magnetic conductance；
Further, it is to the detailed process of complicated Distribution of Magnetic Field domain mesh：According to real electrical machinery design size, consider pole
Air gap, rotor teeth portion use different number of plies nets air gap, rotor leakage field between boots at air gap includes pole shoe between stator pole shoes
The mode of lattice subdivision establishes according to subdivision result the magnetic conductance connection relation in each grid to obtain suitable subdivision model.
Further, the detailed process of the step 3 is：To the magnetic conductance model modeling of stator yoke and teeth portion, rotor yoke
Using traditional approach, by each stator tooth, the stator yoke of between cog, rotor yoke and each permanent magnet blocks, be respectively seen as one it is whole
Body magnetic conductance, and according to Distribution of Magnetic Field rule, it is sequentially connected magnetic conductance node in each region.
Further, the detailed process of the step 4 is：The equivalent magnetic of main permanent magnet in pole shoe Halbach permanent magnet arrays
It leads and is connected respectively with magnetic conductance at both sides pair permanent magnet magnetic conductance, upside pole shoe and downside air gap, secondary permanent magnet magnetic conductance and rotor
Permanent magnet magnetic conductance is corresponding with adjacent subdivision grid magnetic conductance to be connected；The magnetic conductance of stator and air gap intersection divided region is according to quantity
Relationship is correspondingly connected with；Rotor teeth portion magnetic conductance and the connection relation of air-gap permeance needs are constantly updated in motor rotation process.
Further, in the step 5, in the step 5, when solution, is modified the magnetic conductivity of iron core magnetic conductance,
Iterative formula is：μ(k)=k1×μ(k-1)+(1-k1)×μ(k), wherein k1Meet 0<k1<1.And due to the iterative calculation meeting of magnetic conductivity
Generate the difference △ B=max { B of teeth portion flux densityi (k)-Bi-1 (k), wherein i is the motor stator number of teeth, is 12.k1Value according to △
The calculated value of B changes and changes, k1=k2× △ B, k2Take 0.008.
Further, in the step 6, magnetic circuit magnetic conductance model constraint formulations are Ф=FG；To the upper of permanent magnet and winding
The matrix of vertical Ф is set up in bottom, while establishing G matrix to magnetic circuit magnetic conductance in motor, and uses iterative method, inquiry BH curve ginseng
Number, when being iterated to calculate to stationary value, i.e. △ B≤0.5% using the update of magnetic conductivity, you can obtain the magnetic potential of each node.
Further, in the step 7, according to each node magnetic potential for solving the acquisition of magnetic circuit matrix, electrical angle can be obtained
The motor stator tooth at a moment flows through magnetic flux in period.Again next rotor-position is calculated, so obtains one
The magnetic linkage Ф of electrical angle period each tooth obtains the electromagnetic parameters such as motor three-phase magnetic flux, induced electromotive force accordingly.
The invention has the advantages that：
1, the present invention is in modeling, and magnetic saturation and leakage field are complete caused by the stator permanent magnet of air gap both sides with rotor permanent magnet
It takes into account, is allowed to approach with FEA Method And Principles, helps to improve model accuracy in computation.
2, in the present invention, rule magnetic circuit carries out the equivalent of magnetic circuit magnetic conductance according to the rule range of the magnetic line of force, simultaneously because electric
Machine internal magnetic field has symmetry, and symmetric position magnetic circuit magnetic conductance also has very high similarity, while not reducing precision
Simplification is programmed with lift scheme.
3, in the present invention, complicated magnetic field divides magnetic conductance size of mesh opening according to core structure parameter, not due to electric machine rotor
It is all different with positional structure, therefore uses three kinds of different grids and carry out subdivision, made to effectively reduce grid mismatch
At error.
4, the program that convergence factor is adjusted with the variation of worst error value is completed, it is unnecessary compared on traditional algorithm to make always
It uses fixed constant as convergence factor, the iteration of Matrix Solving can be completed more quickly.
Description of the drawings
Fig. 1 is the 2D structure charts of motor used in the present invention；
Fig. 2 is the magnetic line of force analysis chart of motor used in the present invention；
Fig. 3 is the structure chart of motor permanent magnet magnetizing direction；
Fig. 4 is the magnetic network model structure chart of the invention with mesh generation；
Fig. 5 is the magnetic network model partial enlargement structure chart of the invention with mesh generation；
Fig. 6 is the magnetic network model magnetic field complex region structure chart of the invention with mesh generation；
Fig. 7 is finite element software and the single-phase unloaded magnetic flux contrast schematic diagram with mesh generation magnetic network model.
Fig. 8 is the flow chart of the modeling method of the present invention.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes.
Legibly illustrate beneficial effects of the present invention in order to simpler, it is double with reference to a specific rotor
Permanent magnetism vernier motor is described in detail：Fig. 1 opens up benefit structure chart for the motor, and 1 is stator yoke in figure, and 2-1 is fixed
Sub- teeth portion；2-2 is stator teeth pole shoe, and 2-3 air gaps between stator teeth pole shoe, 3 be winding, and 4-1 is rotor yoke, and 4-2 is
Rotor permanent magnet, 4-3 are rotor teeth portion, and 5-1 is stator Halbach the first permanent magnets of permanent magnet array, and 5-2 is stator Halbach
The second permanent magnet of permanent magnet array, 5-3 be stator Halbach permanent magnet array third permanent magnets, 6 between rotor air gap；The present invention
Embodiment is the three phase electric machine of the pole of 12 slots/28, is divided into stator, four part of rotor, air gap and shaft；It include stator yoke in stator
Portion, stator teeth, stator slot, armature winding and pole shoe permanent magnet array, armature flute profile are flat-bottom slot, and armature winding is using concentration
Formula winding method, span are 4 stator slots；Stator permanent magnet material is NdFe35, has 2 notches in each pole shoe of stator,
For each dress there are one Halbach permanent magnet arrays, each Halbach permanent magnet arrays are rushed main permanent magnet and two lateral incisions of magnetic by intermediate radial
It is constituted to the secondary permanent magnet for rushing magnetic, the distribution situation of air-gap field can be effectively improved；Rotor is cylindrical shape, peace of slotting on surface
The embedding permanent magnet of table is filled, permanent magnet material is ferrite Y30, and the embedding permanent magnet cross section of table is trapezoidal, is evenly distributed on rotor circumference
Direction；The material of stator core and rotor core is silicon steel sheet DW540_50；
Flow chart as shown in Figure 8 is divided into following steps realization：
Step 1, motor complexity field region and rule field region are divided using finite element software.
Fig. 2 is magnetic line of force analysis chart in the magnetic field of motor of the embodiment of the present invention.In motor model, the area of motor-field complexity
Domain is concentrated mainly on three stator teeth pole shoe, air gap and rotor tooth parts, and the region shape is without consistency, magnetic line of force tendency
It is complicated irregular, it is also easy to produce between cog leakage field and the winding magnetic line of force；The simple region of motor-field is concentrated mainly on stator yoke, determines
Three parts of sub- teeth portion and rotor yoke, the position shape is relatively regular, and magnetic line of force tendency rule is substantially towards the same side
To and magnetic line of force quantity remains unchanged on a magnetic circuit, while stator armature slot leakage is led minimum, is not likely to produce leakage field.
Step 2, the diamond shape small grid subdivision of motor complexity field region is established.
Fig. 4 is the complete 2D magnetic circuits subdivision model of the embodiment of the present invention, and Fig. 6 is the diamond shape that motor amplifies complicated field region
Small grid subdivision schematic diagram；Include pole to stator pole shoes according to the leakage field between real electrical machinery design size, consideration pole shoe at air gap
Air gap, rotor teeth portion use the mode of different number of plies mesh generations to obtain suitable subdivision model between air gap, rotor between boots,
And the magnetic conductance connection relation in each grid is established according to subdivision result.The region of stationary part magnetic field complexity is stator teeth pole
Air gap 2-3 between boots 2-2 and pole shoe.In the stator pole shoes 2-2 of each tooth of motor, all it is embedded in two groups and uses Halbach array
Permanent magnet, be made of the first permanent magnet 5-1, the second permanent magnet 5-2 and third permanent magnet 5-3.Forever due to embedded Halbach
Magnetic array, grid opsition dependent is distributed in pole shoe, forever with the left side of the second permanent magnet 5-2, the top of the first permanent magnet 5-1 and third
The right side of magnet 5-3 is connected.Simultaneously because whole magnetizing direction is downward, the second permanent magnet 5-2 and third permanent magnet 5-3, with the
The upper point of one permanent magnet 5-1 is connected.There are air gap 2-3 between pole shoe and pole shoe simultaneously.Due to pole shoe magnetic field complexity, the magnetic line of force
Trend is mixed and disorderly and there are leakage fields, therefore carries out small grid subdivision to pole shoe, and equal grid height-width is l1.From pole shoe 2-2
Bottom to top is divided into 5 layers of rectangular mesh, and since permanent magnet array magnetizing direction is fixed, by net at permanent magnetism body position
Grillages removes, and traditional approach is used to permanent magnet, regards each piece of permanent magnet all as a whole magnetic conductance, permanent magnet array height is 4l1。
Therefore pole shoe 2-2 mesh generations share 5 rows, and the first row is made of 28 grids, and following four rows is made of per a line 12 grids,
It is divided into three sections by permanent magnet, each section there are 4 row grids.Air gap 2-3 between pole shoe equally uses an equal amount of 5 row grid, often
A line is made of 3 grids, and both sides are connected respectively at adjacent pole shoe grid.
Air gap is the important place of energy exchange and the region of magnetic field complexity the most.By rationally designing, to air gap 6
Using three layers of subdivision, all there are 1800 grid protocol rings per a line, grid height-width is equal, is l3.Air gap 6 simultaneously
Sizing grid is approximately the 1/5 of stator pole shoes 2-2 sizing grids, that is, therefore each air gap 2-3 between pole shoe 2-2 and pole shoe is most
Lower layer's grid connects 5 air gap grids within the scope of mesh width.For Halbach permanent magnet arrays, due to tangentially filling in practice
The permanent magnet of magnetic will not be formed down the magnetic line of force, thus the second permanent magnet 5-2 and third permanent magnet 5-3 not with 6 grid phase of air gap
Even, only the first permanent magnet 5-1 magnetic conductances lower end connects 6 grid of air gap in all first permanent magnet 5-1 width ranges.
The region of rotor portion magnetic field complexity is rotor teeth portion 4-3 and rotor permanent magnet 4-2.For rotor teeth portion 4-3, by
In partial magnetic field complexity, magnetic line of force distribution is irregular, therefore carries out mesh generation, mesh width height to rotor teeth portion 4-3
It is equal, it is l2.Inside rotor teeth portion 4-3, the mesh generation mode arranged using 4 rows 5, due to the boundary of rotor teeth portion 4-3
It is not completely radial, there is certain deviation angle, therefore increase one on the boundary of every row magnetic conductance grid, according to rotor teeth portion side
Cut the calculated equivalent tangential magnetic conductance of area ratio of grid in boundary.Similarly, since rotor permanent magnet 4-2 magnetizing directions are downward,
And magnetizing direction fixes basic zero deflection, and traditional magnetic circuit model, and rotor permanent magnet 4-2 magnetic are used to rotor permanent magnet 4-2
Upper end is led with upper three row of left and right sides rotor teeth portion 4-3 to be connected.Since rotor teeth portion 4-3 and rotor permanent magnet 4-2 are being pasted
Nearly rotor tooth bottom, there are a little flux loops, therefore by design at rotor permanent magnet 4-2 and rotor teeth portion 4-3, after
One layer of grid identical with rotor teeth portion size of mesh opening of continuous addition, to lay out the winding magnetic circuit between teeth portion and permanent magnet.Rotor
The equivalent magnetic conductance lower ends permanent magnet 4-2 are connected with the grid in permanent magnet width range, rotor teeth portion 4-3 grids and next layer of net
Lattice, which are sequentially connected, to be connect.
Permeance calculation formula is as follows：
Step 3, the magnetic circuit magnetic conductance for establishing motor rule field region is equivalent.
For the stationary part of motor, since stator yoke 1 and stator teeth 2-1 magnetic line of force tendencies are smooth, no winding shape
The stator yoke of each stator tooth and between cog is regarded this magnetic conductance model modeling at two using traditional approach by the magnetic line of force respectively
For a whole magnetic conductance, and it is sequentially connected.And since motor stator pole shoe uses mesh generation, stator teeth 2-
1 magnetic conductance model is connected with 28 grids of pole shoe 2-2 the first rows under corresponding tooth.
Rotor yoke portion 4-1 is similar to stator yoke 1, and magnetic line of force path is circulated in tangent line, and smooth nothing is returned
Ring, therefore by the parts rotor yoke 4-1 between each rotor teeth portion 4-1 and adjacent rotor permanent magnet 4-2, using traditional magnetic
Road model is equivalent to a whole magnetic conductance, and one, which shares 56 yoke portion magnetic conductances, connects to form rotor yoke annulus.
Step 4, magnetic conductance node connection relation at two class zone boundaries is improved, builds complete 2D magnetic circuits subdivision model, successively
Establish magnetic conductance solution matrix equation.
In pole shoe Halbach permanent magnet arrays the equivalent magnetic conductance of main permanent magnet respectively with both sides pair permanent magnet magnetic conductance, upside
Pole shoe is connected with magnetic conductance at the air gap of downside, secondary permanent magnet magnetic conductance and rotor permanent magnet magnetic conductance with adjacent subdivision grid magnetic conductance
It is corresponding to be connected；For air gap 6, since rotor is rotated relative to stator, i.e. stator pole shoes 2-2, stator pole shoes air gap 2-3 and first
Permanent magnet 5-1 and the connection type of air gap 6 immobilize, and rotor changes from the connection of air gap with the different of rotation angle
Become.It is thus determined that the connection method of rotor and air gap, is the key that establish rotation magnetic network model.Rotor is in different location
When, 4 groups of grids of the first row of rotor teeth portion 4-3 are connected with the air gap grid in corresponding range respectively, in rotor tooth 4-3 nets
About 6 to 7 air gap grids are attached thereto within the scope of lattice.Rotor permanent magnet 4-2 also with the air gap grid phase in corresponding range
Connection, due to being modeled using traditional magnetic circuit to rotor permanent magnet 4-2, covering air gap range is larger, in rotor permanent magnet 4-2
About 32 6 grids of air gap are connected in range.Whenever rotor turns over a set angle, rotor permanent magnet 4-2 and
All rotor tooth 4-3 grids can all rotate the angle along Z axis, and footprint size is constant, but coverage area connected to it
Interior 6 mesh node of air gap changes.When rotor turns over different angle, the magnetic conductance of stator and air gap intersection divided region
It is correspondingly connected with according to quantitative relation；Rotor teeth portion magnetic conductance and the connection relation of air-gap permeance need in motor rotation process constantly
Update.Rotor teeth portion 4-3 grids and rotor permanent magnet 4-2 to the different node order of connection of the grid of air gap 6, update rotor with
The Node connectedness relationship of air gap.
Stator, air gap, each magnetic circuit magnetic conductance node of rotor are connected by step 1, step 2, step 3 with step 4, can be built
Complete 2D magnetic circuit subdivision models, and magnetic conductance solution matrix equation is established successively, calculation formula is as follows：
GF=Q (2)
In formula：
And
F=[F (1) ... F (7672)]T
Q=[Q (1) ... Q (7672)]T
Node magnetic potential F meets：
F=G-1·Q (3)
Under MATLAB environment, directly solved with following formula：
F=G Q (4)
Flux density calculation formula is：
B indicates that flux density, s, t are the node at magnetic conductance both ends in formula, and S is that the effective cross section in magnetic conductance region is accumulated.
Step 5, the relational expression of magnetic conductivity convergence factor and motor teeth portion flux density maximum difference is built；
The solution procedure of matrix equation (2) is nonlinear iterative process, need to be solved, be used with non-linear iterative
Newton iteration method formula be：
μ(k)=k1×μ(k-1)+(1-k1)×μ(k) (6)
In formula, magnetic permeability μ is iterative calculation parameter, k1For coefficient, k is iterations；Wherein k1Meet 0<k1<1.And by
The difference △ B=max { B of teeth portion flux density are will produce in the iterative calculation of magnetic conductivityi (k)-Bi-1 (k), wherein i is motor stator tooth
Number is 12.k1Value according to the calculated value of △ B change and change, k1=k2× △ B, k2Take 0.008.By by teeth portion flux density
Calculating difference △ B and coefficient k1Connection generates the functional relation changed with iteration, and when difference is larger, coefficient value is larger,
Result of calculation is set to be affected by last time calculating；When difference is smaller, coefficient value is smaller, and result of calculation is made to be counted by the last time
It is smaller to calculate influence.Accelerate the iteration speed of entire program, improves efficiency.
Step 6, magnetic conductance matrix equation is solved, realizes that the rapid solving of matrix is obtained through each magnetic conductance section using iterative algorithm
The magnetic flux and magnetic potential of point, and then calculate the magnetic flux density and magnetic conductivity of each magnetic conductance；Magnetic circuit magnetic conductance model constraint formulations be Ф=
F·G；The matrix of Ф is established to the upper and lower position of permanent magnet and winding, while G matrix is established to magnetic circuit magnetic conductance in motor, and is adopted
With iterative method, BH curve parameter is inquired, is iterated to calculate to stationary value using the update of magnetic conductivity.
Linear interpolation is used according to the magnetic conductance matrix equation established in step 4, while to BH curve：
When teeth portion flux density B is satisfied the following conditional expression, then represents an iteration and calculate completion：
ζ is worst error value in formula, and value takes 0.5%；
Step 7, it according to the magnetic circuit parameters solved in step 6, calculates and constrains according to motor electromagnetic, each phase of motor can be calculated
The electromagnetic parameters such as magnetic linkage, winding counter electromotive force and output torque.
The each node magnetic potential obtained according to magnetic circuit matrix is solved, can obtain the motor at a moment in the electrical angle period
Stator tooth flows through magnetic flux, calculates again next rotor-position, so obtains the magnetic of electrical angle period each tooth
Chain Ф obtains the electromagnetic parameters such as motor three-phase magnetic flux, induced electromotive force accordingly, can be used for calculating motor if band load situation defeated
Go out torque.In order to which the 2D magnetic circuit subdivision modeling methods for verifying a kind of double permanent magnetism vernier motors of rotor proposed by the invention are accurate
Reliably, Fig. 7 gives simulation result, and the result obtained with finite element business software carries out contrast verification.
It is illustrated in figure 7 using business finite element analysis software and using the equivalent magnetic network model point of multi grid subdivision
The single-phase magnetic flux of analysis compares figure, and A1 is the single-phase unloaded flux waveforms that business finite element analysis software obtains in figure, and B1 is multiple
The single-phase unloaded flux waveforms that mesh generation equivalent magnetic network model obtains, it can be seen that in the low calculating analysis consumption of equivalent magnetic network
In the case of, the characteristics of ensure that high actuarial precision and accuracy.
To sum up, the 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of a kind of rotor of the invention, including to motor
The division of complicated field region and rule field region, and diamond shape small grid subdivision and magnetic circuit magnetic conductance are carried out to two kinds of regions respectively
It is equivalent；Magnetic conductance node connection relation at two class zone boundaries is improved, complete 2D magnetic circuits subdivision model is built, magnetic conductance is established successively and asks
Dematrix equation；The relational expression for building magnetic conductivity convergence factor and motor teeth portion flux density maximum difference, to improve Matrix Solving
Speed, so as to optimization design；Magnetic conductance matrix equation is solved, realizes that the rapid solving of matrix is obtained through each magnetic using iterative algorithm
The magnetic flux and magnetic potential of node are led, and then calculates the magnetic flux density and magnetic conductivity of each magnetic conductance, and calculates and constrains according to motor electromagnetic,
The electromagnetic parameters such as each phase magnetic linkage of motor, winding counter electromotive force can be calculated, are finally compared with the result of finite element.It is of the invention first
Secondary to carry out 2D magnetic circuit subdivisions for the double permanent magnetism vernier motors of rotor, the scheme provided can be that the double magnetoes of the type carry
It studies for reference.
Although the present invention has been described by way of example and in terms of the preferred embodiments, embodiment is not for the purpose of limiting the invention.Not
It is detached from the spirit and scope of the present invention, any equivalent change or retouch done belongs to the application appended claims institute
The protection domain of restriction.
Claims (9)
1. a kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor, which is characterized in that include the following steps：
Step 1, motor complexity field region and rule field region are divided using finite element software；
Step 2, the diamond shape small grid subdivision of motor complexity field region is established；
Step 3, the magnetic circuit magnetic conductance for establishing motor rule field region is equivalent；
Step 4, magnetic conductance node connection relation at two class zone boundaries is improved, complete 2D magnetic circuits subdivision model is built, establishes successively
Magnetic conductance solution matrix equation；
Step 5, the relational expression of magnetic conductivity convergence factor and motor teeth portion flux density maximum difference is built；
Step 6, magnetic conductance matrix equation is solved, realizes that the rapid solving of matrix is obtained through each magnetic conductance node using iterative algorithm
Magnetic flux and magnetic potential, and then calculate the magnetic flux density and magnetic conductivity of each magnetic conductance；
Step 7, it according to the magnetic circuit parameters solved in step 6, calculates and constrains according to motor electromagnetic, each phase magnetic of motor can be calculated
The electromagnetic parameters such as chain, winding counter electromotive force.
2. a kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor according to claim 1, feature
It is, the double permanent magnetism vernier motors of the rotor are the three phase electric machine of the pole of 12 slots/28, are divided into stator, air gap, rotor and shaft four
Part；Include stator yoke, stator teeth, stator slot, armature winding and pole shoe permanent magnet array in stator, armature flute profile is flat
Kerve, for armature winding using centralized winding method, span is 4 stator slots；Stator permanent magnet material is NdFe35, in stator
Each pole shoe has 2 notches, and there are one Halbach permanent magnet arrays, each Halbach permanent magnet arrays are rushed each dress by intermediate radial
The secondary permanent magnet that the main permanent magnet of magnetic and both sides tangentially rush magnetic is constituted；Rotor is cylindrical shape, and installation table of slotting on surface is embedding forever
Magnet, permanent magnet material are ferrite Y30, and the embedding permanent magnet cross section of table is trapezoidal, is evenly distributed on rotor circumference direction；Stator
The material of iron core and rotor core is silicon steel sheet DW540_50；Between stator and rotor, air gap thickness is air gap
0.5mm；Machine shaft is made of un-conducted magnetic material, is solid cylindrical, and connect with rotor coaxial.
3. a kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor according to claim 2, feature
It is, in the step 1, according to the simulation result of finite element software, motor complexity Distribution of Magnetic Field region focuses primarily upon fixed turn
In one segment limit of sub- two-sided permanent magnet body, including at air gap, stator pole shoes and rotor teeth portion region；Motor rule Distribution of Magnetic Field area
Domain includes stator yoke and teeth portion, rotor yoke region, since both sides permanent magnet magnetizing direction has been fixed, permanent magnet inner magnet
Field distribution is also that rule is distributed.
4. a kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor according to claim 2, feature
It is, the detailed process to complicated Distribution of Magnetic Field domain mesh is：According to real electrical machinery design size, consider between pole shoe at air gap
Leakage field, air gap air gap, rotor, rotor teeth portion use the side of different number of plies mesh generations including pole shoe between stator pole shoes
Formula establishes according to subdivision result the magnetic conductance connection relation in each grid to obtain suitable subdivision model.
5. a kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor according to claim 2, feature
It is, the detailed process of the step 3 is：Tradition side is used to the magnetic conductance model modeling of stator yoke and teeth portion, rotor yoke
Each stator tooth, the stator yoke of between cog, rotor yoke and each permanent magnet blocks are respectively seen as a whole magnetic conductance, and press by formula
According to Distribution of Magnetic Field rule, it is sequentially connected magnetic conductance node in each region.
6. a kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor according to claim 2, feature
It is, the detailed process of the step 4 is：The equivalent magnetic conductance of main permanent magnet is respectively with two in pole shoe Halbach permanent magnet arrays
Side pair permanent magnet magnetic conductance, upside pole shoe are connected with magnetic conductance at the air gap of downside, secondary permanent magnet magnetic conductance and rotor permanent magnet magnetic conductance
It is corresponding with adjacent subdivision grid magnetic conductance to be connected；The magnetic conductance of stator and air gap intersection divided region is corresponded to according to quantitative relation to be connected
It connects；Rotor teeth portion magnetic conductance and the connection relation of air-gap permeance needs are constantly updated in motor rotation process.
7. a kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor according to claim 1, feature
It is, in the step 5, when solution is modified the magnetic conductivity of iron core magnetic conductance, and iterative formula is：μ(k)=k1×
μ(k-1)+(1-k1)×μ(k), wherein k1Meet 0<k1<1, and since the iterative calculation of magnetic conductivity will produce the difference △ of teeth portion flux density
B=max { Bi (k)-Bi-1 (k), wherein i is the motor stator number of teeth, is 12, k1Value according to the calculated value of △ B change and change,
k1=k2× △ B, k2Take 0.008.
8. a kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor according to claim 1, feature
It is, in the step 6, magnetic circuit magnetic conductance model constraint formulations are Ф=FG；The upper and lower position of permanent magnet and winding is established
The matrix of Ф, while G matrix is established to magnetic circuit magnetic conductance in motor, and iterative method is used, BH curve parameter is inquired, magnetic conductance is utilized
When the update of rate is iterated to calculate to stationary value, i.e. △ B≤0.5%, you can obtain the magnetic potential of each node.
9. a kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor according to claim 1, feature
It is, in the step 7, according to each node magnetic potential for solving the acquisition of magnetic circuit matrix, one can be obtained in the electrical angle period
The motor stator tooth at moment flows through magnetic flux, calculates again next rotor-position, so obtains an electrical angle period
The magnetic linkage Ф of each tooth, obtains the electromagnetic parameters such as motor three-phase magnetic flux, induced electromotive force accordingly, and carrying situation if band can be used for
Calculate motor output torque.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810263411.8A CN108595772B (en) | 2018-03-26 | 2018-03-26 | 2D magnetic circuit subdivision modeling method for stator and rotor double-permanent-magnet vernier motor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810263411.8A CN108595772B (en) | 2018-03-26 | 2018-03-26 | 2D magnetic circuit subdivision modeling method for stator and rotor double-permanent-magnet vernier motor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108595772A true CN108595772A (en) | 2018-09-28 |
| CN108595772B CN108595772B (en) | 2020-03-31 |
Family
ID=63624725
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810263411.8A Active CN108595772B (en) | 2018-03-26 | 2018-03-26 | 2D magnetic circuit subdivision modeling method for stator and rotor double-permanent-magnet vernier motor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108595772B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110705088A (en) * | 2019-09-27 | 2020-01-17 | 东南大学 | Permanent magnet motor modeling and electromagnetic performance calculating method based on magnetic network |
| CN111190128A (en) * | 2018-11-15 | 2020-05-22 | 北京自动化控制设备研究所 | Detection algorithm for BH characteristics of ferromagnetic material of reluctance motor |
| EP3660522A1 (en) * | 2018-11-30 | 2020-06-03 | Zhejiang University | Method for evaluating electromagnetic performace of permanent-magnet machines |
| CN112491169A (en) * | 2020-11-10 | 2021-03-12 | 东南大学 | Stator magnetic-gathering type bilateral permanent magnet motor |
| CN113285536A (en) * | 2021-07-19 | 2021-08-20 | 北京航空航天大学 | Novel composite rotor high-speed permanent magnet motor equivalent magnetic circuit model construction method |
| CN113746385A (en) * | 2021-07-28 | 2021-12-03 | 华中科技大学 | Permanent magnet motor topology construction method based on discrete magnetic conductance model |
| US11574082B1 (en) * | 2019-12-27 | 2023-02-07 | Ansys, Inc. | Partial simulation model creation, simulation and post-processing for rotational machines |
| CN116542065A (en) * | 2023-05-16 | 2023-08-04 | 江苏大学 | Electromagnetic force and vibration modeling method for surface-mounted permanent magnet motor |
| CN116992723A (en) * | 2023-07-31 | 2023-11-03 | 重庆理工大学 | Motor dynamic magnetic network modeling method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103795159A (en) * | 2014-01-03 | 2014-05-14 | 东南大学 | Stator-rotor double-permanent-magnet-type-vernier motor |
| CN103853891A (en) * | 2014-03-21 | 2014-06-11 | 西北工业大学 | Finite element analysis-based variable-element permanent magnet synchronous motor modeling method |
| CN104063556A (en) * | 2014-07-07 | 2014-09-24 | 东南大学 | Modeling method of axial permanent magnetic motor equivalent magnetic circuit model |
| CN106682375A (en) * | 2017-03-17 | 2017-05-17 | 福州大学 | Bearing-free flux-switching motor magnetic circuit analysis method considering rotor eccentricity |
| CN107565782A (en) * | 2017-08-31 | 2018-01-09 | 江苏大学 | A kind of equivalent magnetic network analysis method of mixed rotor permagnetic synchronous motor |
-
2018
- 2018-03-26 CN CN201810263411.8A patent/CN108595772B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103795159A (en) * | 2014-01-03 | 2014-05-14 | 东南大学 | Stator-rotor double-permanent-magnet-type-vernier motor |
| CN103853891A (en) * | 2014-03-21 | 2014-06-11 | 西北工业大学 | Finite element analysis-based variable-element permanent magnet synchronous motor modeling method |
| CN104063556A (en) * | 2014-07-07 | 2014-09-24 | 东南大学 | Modeling method of axial permanent magnetic motor equivalent magnetic circuit model |
| CN106682375A (en) * | 2017-03-17 | 2017-05-17 | 福州大学 | Bearing-free flux-switching motor magnetic circuit analysis method considering rotor eccentricity |
| CN107565782A (en) * | 2017-08-31 | 2018-01-09 | 江苏大学 | A kind of equivalent magnetic network analysis method of mixed rotor permagnetic synchronous motor |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111190128A (en) * | 2018-11-15 | 2020-05-22 | 北京自动化控制设备研究所 | Detection algorithm for BH characteristics of ferromagnetic material of reluctance motor |
| CN111190128B (en) * | 2018-11-15 | 2022-10-18 | 北京自动化控制设备研究所 | Detection algorithm for BH characteristics of ferromagnetic material of reluctance motor |
| EP3660522A1 (en) * | 2018-11-30 | 2020-06-03 | Zhejiang University | Method for evaluating electromagnetic performace of permanent-magnet machines |
| US12040659B2 (en) | 2018-11-30 | 2024-07-16 | Zhejiang University | Method for evaluating electromagnetic performance of electric machines in particular of permanent-magnet machines |
| CN110705088B (en) * | 2019-09-27 | 2022-11-18 | 东南大学 | Permanent magnet motor modeling and electromagnetic performance calculating method based on magnetic network |
| CN110705088A (en) * | 2019-09-27 | 2020-01-17 | 东南大学 | Permanent magnet motor modeling and electromagnetic performance calculating method based on magnetic network |
| US11574082B1 (en) * | 2019-12-27 | 2023-02-07 | Ansys, Inc. | Partial simulation model creation, simulation and post-processing for rotational machines |
| CN112491169A (en) * | 2020-11-10 | 2021-03-12 | 东南大学 | Stator magnetic-gathering type bilateral permanent magnet motor |
| CN113285536A (en) * | 2021-07-19 | 2021-08-20 | 北京航空航天大学 | Novel composite rotor high-speed permanent magnet motor equivalent magnetic circuit model construction method |
| CN113285536B (en) * | 2021-07-19 | 2021-11-02 | 北京航空航天大学 | Method for constructing equivalent magnetic circuit model of composite rotor high-speed permanent magnet motor |
| CN113746385A (en) * | 2021-07-28 | 2021-12-03 | 华中科技大学 | Permanent magnet motor topology construction method based on discrete magnetic conductance model |
| CN113746385B (en) * | 2021-07-28 | 2023-10-24 | 华中科技大学 | Permanent magnet motor topology construction method based on discrete magnetic conduction model |
| CN116542065A (en) * | 2023-05-16 | 2023-08-04 | 江苏大学 | Electromagnetic force and vibration modeling method for surface-mounted permanent magnet motor |
| CN116542065B (en) * | 2023-05-16 | 2024-09-03 | 江苏大学 | Electromagnetic force and vibration modeling method for surface-mounted permanent magnet motor |
| CN116992723A (en) * | 2023-07-31 | 2023-11-03 | 重庆理工大学 | Motor dynamic magnetic network modeling method |
| CN116992723B (en) * | 2023-07-31 | 2024-01-16 | 重庆理工大学 | Motor dynamic magnetic network modeling method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108595772B (en) | 2020-03-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108595772A (en) | A kind of 2D magnetic circuit subdivision modeling methods of the double permanent magnetism vernier motors of rotor | |
| WO2021237848A1 (en) | Parametric equivalent magnetic network modeling method for multi-objective optimization of permanent magnet electric motor | |
| CN109600006B (en) | Solving method for electromagnetic design of surface-mounted permanent magnet motor | |
| CN107565782B (en) | A kind of equivalent magnetic network analysis method of mixed rotor permanent magnet synchronous motor | |
| Zhang et al. | Investigation of auxiliary poles optimal design on reduction of end effect detent force for PMLSM with typical slot–pole combinations | |
| CN110705088B (en) | Permanent magnet motor modeling and electromagnetic performance calculating method based on magnetic network | |
| Sun et al. | Electromagnetic analysis on novel rotor-segmented axial-field SRM based on dynamic magnetic equivalent circuit | |
| CN103116704B (en) | A kind of based on the parameterized continuous tide computing method of local geometric | |
| CN108563912A (en) | A kind of analytic method of durface mounted permanent magnet synchronous motor air-gap field | |
| CN109684775B (en) | Online magnetic flux regulation performance prediction and optimization design method of magnetic flux controllable memory motor based on nonlinear equivalent variable magnetic network model | |
| CN107786062A (en) | A kind of segmented cylinder-type transverse-flux linear motor | |
| Ghods et al. | Equivalent magnetic network modeling of dual-winding outer-rotor Vernier permanent magnet machine considering pentagonal meshing in the air-gap | |
| Ghods et al. | A mesh design technique for double stator linear PM Vernier machine based on equivalent magnetic network modeling | |
| Huang et al. | An improved magnetic equivalent circuit method for segmented-halbach axial-flux permanent magnet machines | |
| Gauthey et al. | Multi-material topology optimization with continuous magnetization direction for permanent magnet synchronous reluctance motors | |
| Li et al. | An improved equivalent magnetic network model of modular IPM machines | |
| CN107786061B (en) | A kind of plate transverse flux linear motor | |
| Poudel et al. | Analytical investigation and heuristic optimization of surface mounted permanent magnet machines with hybrid magnetic structure | |
| CN103699752A (en) | Coupling method for processing moving boundary problems in electromagnetic field based on edge element method | |
| Zhao et al. | A flux-intensifying permanent magnet synchronous motor and its magnetic equivalent circuit model | |
| Ruzbehi et al. | Topology optimization of a nonlinear electromagnetic actuator based on a local importance measure | |
| Yang et al. | A fast nonlinear equivalent magnetic network model for magnetic jack type control rod drive mechanism in reactor | |
| CN104616863A (en) | Orthogonal iron core type controllable electric reactor with wedge-shaped groove and equivalent magnetic path length calculating method of orthogonal iron core type controllable electric reactor with wedge-shaped groove under direct-current magnetic bias | |
| Gheisari et al. | A general analytical no-load magnetic model for axial flux permanent magnet synchronous machines | |
| CN118428296A (en) | Permanent magnet flat wire motor topology meshing optimization method based on micro-reluctance unit |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhao Wenxiang Inventor after: Cao Donghui Inventor after: Liu Tong Inventor before: Cao Donghui Inventor before: Liu Tong Inventor before: Zhao Wenxiang |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |