CN108494186A - A kind of optimization method improving ferrite assist type synchronous magnetic resistance motor power factor - Google Patents
A kind of optimization method improving ferrite assist type synchronous magnetic resistance motor power factor Download PDFInfo
- Publication number
- CN108494186A CN108494186A CN201810361597.0A CN201810361597A CN108494186A CN 108494186 A CN108494186 A CN 108494186A CN 201810361597 A CN201810361597 A CN 201810361597A CN 108494186 A CN108494186 A CN 108494186A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- rotor
- diameter
- power factor
- magnetic resistance
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02K—DYNAMO-ELECTRIC MACHINES
- H02K15/00—Methods or apparatus specially adapted for manufacturing, assembling, maintaining or repairing of dynamo-electric machines
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Power Engineering (AREA)
- Computational Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Iron Core Of Rotating Electric Machines (AREA)
Abstract
The present invention relates to a kind of optimization methods improving ferrite assist type synchronous magnetic resistance motor power factor, this method is on the basis of the rotor magnetic of designed synchronous magnetic resistance motor hinders structure, add ferrite, by analyzing the equivalent magnetic loop of ferrite assist type synchronous magnetic resistance motor, find the relational expression of power factor and d axis and q axis magnetic resistance, and then optimize ferrite structure using genetic algorithm and Taguchi methods, higher power factor and smaller ferrite volume structure parameter are found, while in view of splitting the influence than coefficient to power factor.This method can greatly reduce occupied calculator memory resource, while reducing the difficulty of mathematical operation compared with FInite Element, and the power factor and realization that can significantly improve Fa SynRM optimize the structural parameters of Fa SynRM.
Description
Technical field
The present invention relates to a kind of optimization methods improving ferrite assist type synchronous magnetic resistance motor power factor.
Background technology
In recent years, with the rise of the prices of raw materials and the increase of product cost pressure, people are in concern motor performance
While, it is also considered that the cost problem of motor.Synchronous magnetic resistance motor (Synchronous Reluctance Machine, below
Abbreviation SynRM) it is low with its manufacturing cost, output-constant operation range is wide, overload capacity is strong and reliability is high the features such as and
It is concerned.However, SynRM power factors are relatively low, the inverter of external bigger is needed, cost is added additional.To improve
The power factor of SynRM, some scholars are inserted into a small amount of permanent magnet in the magnet isolation tank of rotor.Compared to traditional Permanent Magnet and Electric
The dosage of machine, permanent magnet is relatively fewer, therefore this motor is referred to as permanent magnet assist type synchronous magnetic resistance motor
(Permanent Magnet-assisted Synchronous Reluctance Machine, hereinafter referred to as PMaSynRM).
The two is compared, and PMaSynRM average torques, efficiency and power factor all increase.Therefore, permanent magnet assist type synchronous reluctance
Motor, especially with the ferrite assist type synchronous magnetic resistance motor (Ferrite- of the more cheap Ferrite Material of price
Assisted Synchronous Reluctance Machine, hereinafter referred to as Fa-SynRM), it is a kind of inexpensive, high-performance
Motor.Keep Fa-SynRM power factors big so how to design ferrite structure, becomes the difficult point studied in recent years and hot spot
Problem.
Since ferrite assist type synchronous magnetic resistance motor rotor structure is similar with synchronous magnetic resistance motor, in design,
Usually elder generation's optimization design synchronous magnetic resistance motor, then it is put into ferrite in certain magnet isolation tanks of rotor.Document [Wenliang
Zhao,Dezhi Chen,Thomas A.Lipo,Byung-Il Kwon.Performance Improvement of
Ferrite-Assisted Synchronous Reluctance Machines Using Asymmetrical Rotor
Configurations[J].IEEE Transactions on Magnetics,2015,51(11):8108504.] using new
Type asymmetrical rotor structure, improves torque, efficiency and the power factor of ferrite assist type synchronous magnetic resistance motor, and
Under the auxiliary of JMAG-Designer, it is based on two-dimensional finite element method, is detached reluctance torque and electromagnetic torque using osmosis is freezed
It comes, by making electromagnetic torque and reluctance torque reach maximum value under the same current phase angle, and then is significantly promoted
Power factor and efficiency, still, this method increase the workload and production cost of designer, also increase the torque of motor
Pulsation.
Invention content
In view of the deficiencies of the prior art, the technical issues of present invention intends to solve, which is to provide, improves the synchronization of ferrite assist type
The optimization method of reluctance motor power factor.On the basis of the rotor magnetic of designed synchronous magnetic resistance motor hinders structure, add
Add ferrite, by analyzing the equivalent magnetic loop of ferrite assist type synchronous magnetic resistance motor, finds power factor and d axis and q axis magnetic
The relational expression of resistance, and then optimize ferrite structure using genetic algorithm and Taguchi methods finds higher power factor and more
Small ferrite volume structure parameter, while in view of splitting the influence than coefficient to power factor.This method compared with FInite Element,
Occupied calculator memory resource can be greatly reduced, while reducing the difficulty of mathematical operation, Fa- can be significantly improved
The power factor of SynRM, and can realize and the structural parameters of Fa-SynRM are optimized.
Technical solution is used by the present invention solves the technical problem:A kind of raising ferrite assist type synchronization is provided
The optimization method of reluctance motor power factor, this method comprise the concrete steps that:
The first step:Determine the structural parameters of synchronous magnetic resistance motor (SynRM):
Measurement obtains the stator structure parameter of SynRM, including number of stator slots ns, stator rabbet width Sso, stator slot-pitch angle
β;Secondly, measurement obtains the Structural Parameters of its Rotor of SynRM, including rotor number of pole-pairs P, rotor layers of flux barriers number nl, rotor slot opening
Spend Rso, iron core effective length l, the gas length δ between diameter of stator bore and rotor diameter;
The number of stator slots obtained according to above-mentioned measurement is needed according to formula (1) to determine rotor slot number nr,
Determining that rotor slot number is nrAfterwards, rotor slot-pitch angle α is calculated according to formula (2),
Second step:Taguchi optimizes SynRM rotor insulation magnetic and hinders structure:
Gas length δ between diameter of stator bore and rotor diameter and every layer of magnetic barrier width hi, it optimizes, i=1,
2,···,nl:The level range for providing the gas length δ between diameter of stator bore and rotor diameter first, then in each level
Under, change every layer of magnetic barrier width to compare the size of torque pulsation, finds every layer of magnetic barrier width h respectivelyiThe level of Optimal Parameters
Range;Secondly within the scope of the level of gas length δ between diameter of stator bore and rotor diameter respectively to each diameter of stator bore with turn
Magnetic barrier width h under gas length δ between sub- outer diameteriTaguchi optimization designs are carried out, the torque pulsation value under different δ is found,
Obtain the gas length δ between the diameter of stator bore and rotor diameter after Taguchi optimizations;Recycle the stator after Taguchi optimizations
Gas length δ between internal diameter and rotor diameter, in every layer of magnetic barrier width hiRespectively to every layer of magnetic within the scope of the level of Optimal Parameters
Hinder width hiTaguchi optimization designs are carried out, different h are foundiUnder torque pulsation value, obtain Taguchi optimization after every layer of magnetic
Hinder width hi;Finally determine every layer of magnetic barrier width h after corresponding Taguchi optimizationsiUnder rotor insulation magnetic hinder end spacing tj
And rotor magnetic hinders following arc length xi, j=1,2, nl+1;
Third walks:The rotor structure of ferrite assist type synchronous magnetic resistance motor (Fa-SynRM) optimizes:
The total magnetic resistance R of 3-1, d, q axis magnetic circuitd、RqEquivalent Calculation:
It hinders to meet parallel relationship between each magnetic resistance of each layer magnetic conduction silicon steel sheet magnetic flux circulation, passes through formula (3) equivalent meter
Calculate the total magnetic resistance R of d axis magnetic circuitsd:
Rd=Rd1//Rd2//…//Rdj//… (3)
Wherein, j expressions jth layer magnetic conduction silicon steel sheet, j=1,2, nl+1;RdjIt indicates to hinder jth layer magnetic conduction silicon steel sheet magnetic
Through-flow logical air-gap reluctance value, RdjIt is obtained by formula (6):
In formula, μrFor ferritic relative permeability, δ is outside diameter of stator bore and rotor after second step Taguchi optimizations
Gas length between diameter;
Recycle that ' relationship that △-Y ' become the magnetic resistance for equivalent q axis magnetic circuit of changing commanders is simplified, with formula (11) Equivalent Calculation q
The total magnetic resistance R of axis magnetic circuitq:
R=(R11+…+Rk1+ ...) // (R '11+ ...+R 'k1+…) (11)
In formula, k is the corresponding layers of flux barriers number of q axis magnetic circuits, wherein 1 arrives nlWith nlTo 2nlCorrespondent equal, k=2nl;Rk1It is to close
Magnetic flux is hindered to carry over the magnetic resistance for the air gap that rotor magnetic conduction silicon steel sheet passes through between diameter of stator bore and rotor diameter after ferrite is added
RyjWith obstruction across the ferritic magnetic resistance R of each layermiWith obstruction across the magnetic resistance R of each layer barrierδiFunction, with formula (10) indicate;
RyjIt is indicated with formula (7), RmiIt is indicated with formula (8), RδiIt is indicated with formula (9);Similarly, another portion that q axis is symmetrically distributed
The R ' dividedyj, R 'mi, R 'δiFormula (7), (8), (9) are used to indicate respectively;
Wherein, k1It is magnetic flux area correction value, hmiFor ferrite width, bmiFor ferrite length;
3-2, with the total magnetic resistance R of d, q axis magnetic circuitd、RqIndicate power factor:
The total magnetic resistance R of d axis magnetic circuits that step 3-1 is obtaineddWith the total magnetic resistance R of q axis magnetic circuitsqIt brings into formula (12), obtains by d, q
Power factor calculating formula represented by axis equivalent magnetic resistance, the rotor structure which establishes power factor PF and be determined above
Parameter rotor insulation magnetic hinders end spacing tj, insulation magnetic barrier width h between adjacent two layers magnetic conduction silicon steel sheeti, rotor magnetic hinder following arc
Long xi, ferrite width hmi, ferrite length bmi, relationship between gas length δ between diameter of stator bore and rotor diameter, that is, build
The relationship of vertical power factor and Fa-SynRM dependency structure parameters:
3-3, genetic algorithm and Taguchi methods optimize Fa-SynRM rotor structures:
According to formula (12) by genetic algorithm, Fa-SynRM rotor structures are optimized, power factor is obtained
Then the value of each variable under PF maximums carries out local optimum to variable using Taguchi methods, and then is meeting power factor
Ferritic minimum volume value is obtained under conditions of PF is constant;On this basis, further consider to split than coefficient to Fa-SynRM
The influence of power factor further increases the power factor of motor.
Compared with existing optimization method, the present invention has following advantageous effect:
1. the method for existing optimization power factor is to need to combine to realize power factor with FInite Element
Optimization, and FInite Element needs to expend prodigious memory source when calculating, while having sternly to the length for choosing mesh generation
The requirement of lattice.The method of the present invention only needs certain structural parameters sizes of known original motor body structure to be obtained with higher
The Fa-SynRM structures of power factor, therefore the method for the present invention is more simple and practicable, greatly reduces the memory money to computer
The occupancy in source.
2. can more targetedly optimize the power factor of Fa-SynRM structures by the method for the invention, realize power because
Several significantly improves.By the equivalent magnetic circuit for d, q axis established, can accurately and efficiently sum up power factor with it is to be optimized
Structural parameters between relationship, relationship accordingly can treat optimum structural parameter choosing using genetic algorithm in less time
Optimal value is taken, greatly reduces the execution cycle of the power factor of optimization Fa-SynRM structures, and in power factor maximum
In the case of, reduce ferritic volume as far as possible, experiment shows that the optimization method of the application can make ferrite volume ratio
Script it is small by 35.7%, save considerably production cost.
It, will be uneven present in space 3. the method for the present invention establishes the equivalent magnetic circuit of d, q axis using the method for " Chang Hualu "
The distribution shifts of even magnetic flux are equivalent multistage magnetic circuit, and the magnetic flux being approximately considered in every section of magnetic circuit is equal along section and length
Even distribution, and magnetic circuit and circuit are closely similar, i.e., corresponding voltage, electric current, the resistance in circuit is magnetomotive force in magnetic circuit, magnetic
Logical, magnetic resistance.' △-Y ' converter techniques are further simplified equivalent magnetic circuit to recycling, finally establish power factor knot related to Fa-SynRM
The relationship of structure parameter.
The method of the present invention can be adapted for a plurality of types of synchronous magnetic resistance motors or permanent magnet assist type used at present
The optimization of synchronous magnetic resistance motor power factor.
Description of the drawings
Fig. 1 is SynRM overall structure diagrams in the prior art.
Fig. 2 is SynRM stators slot-pitch angle α and rotor slot-pitch angle beta structure schematic diagram in the prior art.
Fig. 3 is outside a kind of embodiment diameter of stator bore and rotor for the method that the present invention improves synchronous magnetic resistance motor torque pulsation
Gas length δ between diameter and torque pulsation and power factor graph of relation.
Fig. 4 is a kind of d of the embodiment for the method that the present invention improves ferrite assist type synchronous magnetic resistance motor power factor
Axis magnetic flux flows to distribution schematic diagram.
Fig. 5 be the present invention improve ferrite assist type synchronous magnetic resistance motor power factor method a kind of embodiment etc.
Imitate d axis magnetic circuits.
Fig. 6 is a kind of q of the embodiment for the method that the present invention improves ferrite assist type synchronous magnetic resistance motor power factor
Axis magnetic flux flows to distribution schematic diagram.
Fig. 7 be the present invention improve ferrite assist type synchronous magnetic resistance motor power factor method a kind of embodiment etc.
Imitate q axis magnetic circuits.
Fig. 8 is that through a step, ' the equivalent q axis obtained after △-Y ' transformation simplifies magnetic circuit to Fig. 7.
Fig. 9 is that through cycle, ' the equivalent q axis obtained after △-Y ' transformation simplifies magnetic circuit to Fig. 7.
Figure 10 is a kind of embodiment for the method that the present invention improves ferrite assist type synchronous magnetic resistance motor power factor
Design sketch after using power factor as object function genetic algorithm optimization.
Figure 11 is a kind of embodiment for the method that the present invention improves ferrite assist type synchronous magnetic resistance motor power factor
Four layers of ferrite width and power factor graph of relation.
Figure 12 is a kind of embodiment for the method that the present invention improves ferrite assist type synchronous magnetic resistance motor power factor
Motor is split than coefficient and power factor and relationship between efficiency curve graph.
Specific implementation mode
The present invention is described in detail below in conjunction with drawings and examples, it should be understood that embodiment described herein is only used for
The description and interpretation present invention, does not limit the protection domain of the application.
The present invention improves the optimization method of ferrite assist type synchronous magnetic resistance motor power factor, the specific steps of this method
It is:
The first step:Determine the structural parameters of synchronous magnetic resistance motor (SynRM):
Measurement obtains the stator structure parameter of SynRM, including number of stator slots ns, stator rabbet width Sso, stator slot-pitch angle
β;Secondly, measurement obtains the Structural Parameters of its Rotor of SynRM, including rotor number of pole-pairs P, rotor layers of flux barriers number nl, rotor slot opening
Spend Rso, iron core effective length l, the gas length δ between diameter of stator bore and rotor diameter;
The number of stator slots obtained according to above-mentioned measurement is needed according to formula (1) to determine rotor slot number nr, to can guarantee
Low torque ripple characteristic is obtained,
Determining that rotor slot number is nrAfterwards, rotor slot-pitch angle α is calculated according to formula (2), ensure that the both ends alignment of magnetic barrier, and
Correspondingly it is aligned with stator slot on the pitch angle of rotor slot;
Second step:Taguchi optimizes SynRM rotor insulation magnetic and hinders structure:
Gas length δ between diameter of stator bore and rotor diameter and every layer of magnetic barrier width hi(i=1,2, nl) it is shadow
The controllable parameter for ringing motor torque ripple, below optimizes the two:The gas between diameter of stator bore and rotor diameter is provided first
The level range of gap length δ changes every layer of magnetic barrier width to compare the size of torque pulsation, respectively then under each level
Find every layer of magnetic barrier width hiThe level range of Optimal Parameters;Secondly gas length δ between diameter of stator bore and rotor diameter
Magnetic barrier width h under gas length δ between each diameter of stator bore and rotor diameter respectively within the scope of leveliCarry out Taguchi
Optimization design finds the torque pulsation value under different δ, obtains the gas between the diameter of stator bore and rotor diameter after Taguchi optimizations
Gap length δ;The gas length δ between the diameter of stator bore and rotor diameter after Taguchi optimizations is recycled, in every layer of magnetic barrier width hi
Respectively to each every layer of magnetic barrier width h within the scope of the level of Optimal ParametersiTaguchi optimization designs are carried out, different h are foundiUnder
Torque pulsation value obtains every layer of magnetic barrier width h after Taguchi optimizationsi;It finally determines every after corresponding Taguchi optimizations
Layer magnetic barrier width hiUnder rotor insulation magnetic hinder end spacing tj(j=1,2, nl+ 1) and rotor magnetic hinders following arc length
xi;
Third walks:The rotor structure of Fa-SynRM optimizes:
The total magnetic resistance R of 3-1, d, q axis magnetic circuitd、RqEquivalent Calculation:
It hinders to meet parallel relationship between each magnetic resistance of each layer magnetic conduction silicon steel sheet magnetic flux circulation, passes through formula (3) equivalent meter
Calculate the total magnetic resistance R of d axis magnetic circuitsd:
Rd=Rd1//Rd2//…//Rdj//… (3)
Wherein, j expressions jth layer magnetic conduction silicon steel sheet, j=1,2, nl+1;RdjIt indicates to hinder jth layer magnetic conduction silicon steel
The air-gap reluctance value of piece magnetic flux circulation, RdjIt is obtained by formula (4):
In formula, μ0For space permeability, value is 4 π × 10-7H/m, δ are the diameter of stator bore after second step Taguchi optimizations
Gas length between rotor diameter;
To simplify formula, the Equivalent Calculation value of d, q axis magnetic circuit resistance is indicated with perunit value, wherein a reference value R* is with formula (5)
It indicates:
In formula, μrFor ferritic relative permeability;
That is, formula (4) is expressed as formula (6) with perunit value:
Magnetic flux is hindered to carry over rotor magnetic conduction silicon steel sheet across the air gap between diameter of stator bore and rotor diameter after ferrite is added
Magnetic resistance RyjFormula (7) can be used to indicate;It hinders across the ferritic magnetic resistance R of each layermiFormula (8) can be used to indicate;It hinders exhausted across each layer
The magnetic resistance R of edge barrierδiIt is indicated according to formula (9);Similarly, another part that q axis is symmetrically distributed hinders magnetic flux to carry over rotor magnetic conduction
Silicon steel sheet passes through the magnetic resistance R ' of the air gap between diameter of stator bore and rotor diameteryjIt is calculated according to formula (7);It hinders across each layer iron
The magnetic resistance R ' of oxysomemiWith obstruction across the magnetic resistance R ' of each layer insulation barrierδiIt is indicated respectively according to formula (8) and formula (9):
Wherein, k1It is magnetic flux area correction value, general value existsBetween, hmiFor ferrite width, bmiFor iron
Oxysome length;
Using ' relationship that △-Y ' become the magnetic resistance for equivalent q axis magnetic circuit of changing commanders is simplified, then through ' △-Y ' transformation can obtain
Ry(j+1)、Ryj、Rδi、RmiBetween meet formula (10):
In formula, the value of k is that (wherein 1 arrives n to the corresponding layers of flux barriers number of q axis magnetic circuitslWith nlTo 2nlCorrespondent equal), i.e. Rk1It is
About the magnetic resistance R for hindering magnetic flux to carry over the air gap that rotor magnetic conduction silicon steel sheet passes through between diameter of stator bore and rotor diameteryjIt is worn with obstruction
Cross the ferritic magnetic resistance R of each layermiWith obstruction across the magnetic resistance R of each layer barrierδiFunction, and it is another to be symmetrically distributed with q axis
Part hinders the magnetic resistance R ' that magnetic flux carries over rotor magnetic conduction silicon steel sheet across the air gap between diameter of stator bore and rotor diameteryjIt is worn with obstruction
Cross the ferritic magnetic resistance R ' of each layermiWith obstruction across the magnetic resistance R ' of each layer insulation barrierδiThe same formula of calculating formula (10);
It is then recycled that ' △-Y ' transformation, can be obtained total magnetic resistance calculating formula of q axis magnetic circuits;With formula (11) Equivalent Calculation
The total magnetic resistance R of q axis magnetic circuitsq:
Rq=(R11+…+Rk1+…)//(R′11+…+R′k1+…) (11)
3-2, with the total magnetic resistance R of d, q axis magnetic circuitd、RqIndicate power factor:
The total magnetic resistance R of d axis magnetic circuits that step 3-1 is obtaineddWith the total magnetic resistance R of q axis magnetic circuitsqIt brings into formula (12), obtains by d, q
Power factor calculating formula represented by axis equivalent magnetic resistance, the rotor structure which establishes power factor PF and be determined above
Parameter rotor insulation magnetic hinders end spacing tj, insulation magnetic barrier width h between adjacent two layers magnetic conduction silicon steel sheeti, rotor magnetic hinder following arc
Long xi, ferrite width hmi, ferrite length bmi, relationship between gas length δ between diameter of stator bore and rotor diameter, that is, build
The relationship of vertical power factor and Fa-SynRM dependency structure parameters:
Wherein:
3-3, genetic algorithm and Taguchi methods optimize Fa-SynRM rotor structures:
According to formula (12) by genetic algorithm, Fa-SynRM rotor structures are optimized, power factor is obtained
Then the value of each variable under PF maximums carries out local optimum to variable using Taguchi methods, and then is meeting power factor
PF obtains ferritic minimum volume value under conditions of being basically unchanged;On this basis, further consider to split than coefficient (i.e. stator
Internal-and external diameter ratio) influence to Fa-SynRM power factors, further increase the power factor of motor.
SynRM and Fa-SynRM in following embodiment are selected for lateral lamination armature motor.
Fig. 1 show the rotor structure of synchronous magnetic resistance motor SynRM in the prior art, and d is synchronous magnetic resistance motor in figure
D-axis, q are the quadrature axis of synchronous magnetic resistance motor, on the center line of rotor magnetic pole are d-axis directions in SynRM, are ahead of counterclockwise
The direction that 45 ° of d axis is quadrature axis direction, and u is rotor diameter, and V indicates machine shaft, a1、a2、a3、a4Respectively four layers of rotor absolutely
Edge barrier, b1、b2、b3、b4、b5Magnetic conduction silicon steel sheet respectively between arbitrary neighborhood insulation barrier, wherein b1For rotor most extension
Magnetic conduction silicon steel sheet, rotor are by four layers of insulation barrier a1、a2、a3、a4Magnetic conduction silicon steel sheet b between arbitrary neighborhood insulation barrier1、
b2、b3、b4、b5Alternating is overrided to form, and magnetic flux flowed through the silicon steel sheet closest to the parts machine shaft V is defined as layer 5 and lead
Magnetic silicon steel sheet b5。
Embodiment illustrated in fig. 2 shows that analysis clear and definite, the present embodiment is made only to take 1/4 of SynRM whole machine models
Point, represented is stator slot-pitch angle β, rotor slot-pitch angle α schematic diagrames, and what wherein β/2 was indicated is the half of stator slot-pitch angle, α/
2 indicate the half of rotor slot-pitch angle.
Fig. 3 show a kind of embodiment the problem of improving synchronous magnetic resistance motor torque pulsation that the present invention first has to solve
Gas length δ between diameter of stator bore and rotor diameter and torque pulsation and the graph of relation of power factor, it is profound using field mouthful
The result figure that one optimization algorithm (Taguchi methods) obtains, wherein abscissa are gas length (mm), and what left ordinate indicated is work(
Rate factor (PF), right ordinate expression is torque pulsation (%).
Embodiment illustrated in fig. 4 shows that analysis clear and definite, the present embodiment is made only to take the 1/4 of Fa-SynRM whole machine models
Part, that outermost two circular arc lines represent is motor stator winding part D, and ten are carried in the D of motor stator winding part
Current direction in the symbology electrical-coil that word intersects is that vertical paper flows into.
The d axis total magnetic flux Φ of Fa-SynRMdIt is each after the air gap passed through between diameter of stator bore and rotor diameter through stator core
Layer magnetic conduction silicon steel sheet is split into each layer magnetic conduction silicon steel sheet magnetic flux Φd1、Φd2..., and they along each layer magnetic conduction silicon steel sheet again
It is returned in stator core again after reaching the air gap between diameter of stator bore and rotor diameter.
Embodiment illustrated in fig. 5 shows the present embodiment using Fa-SynRM as the equivalent d axis magnetic circuit of objective for implementation.Due to magnetic conduction
The relative permeability of silicon steel sheet is far longer than air-gap permeance μ0, therefore each layer magnetic conduction silicon steel sheet magnetic flux Φd1、Φd2... in approach
Magnetic resistance is generated in air gap when each layer magnetic conduction silicon steel sheet mainly between diameter of stator bore and rotor diameter, and full between each magnetic resistance
Sufficient parallel relationship, i.e. R in corresponding diagramd1、Rd2、Rd3、Rd4、Rd5Respectively hinder magnetic conduction silicon steel sheet b1、b2、b3、b4、b5Location
The air-gap reluctance value in domain, the magnetic resistance value can be calculated by formula (6),
In formula, nl=4,
By R obtained aboved1、Rd2、Rd3、Rd4、Rd5Value bring into formula (3), you can obtain the d axis magnetic of the Fa-SynRM
Total magnetic resistance R on roadd,
Rd=Rd1//Rd2//…//Rdj//… (3)
Embodiment illustrated in fig. 6 shows the present embodiment in the flow direction distribution for studying q axis flux, with 1/2 Fa-SynRM electricity
Machine model is research object, and the magnetic flux flow direction of q axis magnetic circuits is complex compared to d axis distribution situations.Space is had ignored in analysis
With time harmonic, leakage magnetic flux, core loss and magnetically saturated situation, as q axis total magnetic flux ΦqAlong stator core and in stator
A part of magnetic flux of air gap can carry over rotor most extension magnetic conduction silicon steel panel region b between diameter and rotor diameter1Direction flow out rotor, then
The secondary air gap between diameter of stator bore and rotor diameter returns in stator core, such as the Φ in Fig. 6q1, Φ 'q1It is shown;Another portion
Point magnetic flux can pass through ferrite to reach from a, a ' point (its indicate be that magnetic flux passes through ferrite or flows out the separation of rotor)
Next layer of magnetic conduction silicon steel sheet, magnetic flux surplus is Φ at this timea, Φ 'a, and Φa、Φ′aMagnetic flux flow direction can equally be divided into two parts,
Process above is repeated in this way, and magnetic flux sequentially circulates.The description for the current segment that motor stator winding part is connected in figure is:
What cross represented be electric current be vertical paper inflow direction, and it is direction that electric current vertically flows out to put represent.
Embodiment illustrated in fig. 7 shows according to Fig. 7 to q axis total magnetic flux ΦqFa- can be obtained by flowing to the analysis of distribution situation
The equivalent q axis magnetic circuit of the rotor of Syn RM.Since the distribution of q axis total magnetic flux is symmetrically distributed with q axis, with the equivalent d axis of analysis
Magnetic circuit is similar, to by way of magnetic conduction silicon steel sheet (b1、b2、b3、b4、b5) region when mainly between diameter of stator bore and rotor diameter
Air gap in generate magnetic resistance Ryj, and each magnetic resistance RyjBetween meet parallel relationship, but unlike d axis magnetic flux distributions, meeting
Some magnetic flux passes perpendicularly through insulation barrier or ferrite reaches next layer of magnetic conduction silicon steel sheet and therefore initially passes through air gap, arrives
Magnetic resistance up to first layer silicon steel sheet is R1, then in two magnetic resistance R of arbitrary neighborhoodyiBetween series connection magnetic flux pass through insulation barrier
Magnetic resistance RδiFerritic magnetic resistance R is passed through with magnetic fluxmi, and magnetic flux passes through the magnetic resistance R of insulation barrierδiIt is passed through with magnetic flux ferritic
Magnetic resistance RmiBetween be parallel relationship.
RyiTo hinder magnetic flux to carry over the magnetic resistance for the air gap that rotor magnetic conduction silicon steel sheet passes through between diameter of stator bore and rotor diameter, use
Formula (7) indicates;It is hindered across the ferritic magnetic resistance R of each layer after ferrite is addedmi(8) can be used to indicate, hinder to insulate across each layer
The magnetic resistance R of barrierδiIt is indicated according to (9).
Wherein:I indicates that i-th layer of barrier, j indicate jth layer magnetic conduction silicon steel sheet, k1It is magnetic flux area correction value, general value
Between, hmiFor ferrite width, bmiFor ferrite length.Again since the symmetry of q axis magnetic circuits can obtain, R 'yj、
R′mi、R′δiCalculating formula respectively with Ryj、Rmi、RδiIt is identical.
Embodiment illustrated in fig. 8 shows the case where magnetic flux distribution caused by the coil flowed by the vertical paper of electric current in Fig. 8
The magnetic resistance series and parallel relationship for analyzing the equivalent q axis magnetic circuit of obtained Fa-SynRM rotors is complex, therefore, using ' △-Y '
Transformation simplifies magnetic circuit, and ' △-Y ' transformation is a ripe method for applying the transformation on circuit, has moved to magnetic circuit here
On, for solving magnetic resistance, ' △-Y ' are transformed to known technology for this.By ' the equivalent magnetic resistance value after △-Y ' transformation is R11、
R12、R13, equivalent magnetic resistance value R11、R12、R13Formula (10) is brought by the corresponding magnetic resistance value that formula (7) to formula (9) is calculated
In following formula is calculated, i.e.,
Similarly, it is passed through in the equivalent magnetic circuit that magnetic flux distribution caused by the coil by the vertical paper outflow of electric current is set up
' the equivalent magnetic resistance value after △-Y ' transformation is R'11、R'12、R'13Value can be also calculated by formula (10).
Then in Fig. 7, magnetic flux distribution is set up caused by the coil by the vertical paper outflow of electric current equivalent magnetic
Road repeatedly utilizes ' △-Y ' with the equivalent magnetic circuit that magnetic flux distribution is set up caused by the coil flowed by the vertical paper of electric current
The transformation results of Fig. 9 can be obtained in transformation for mula, and transformation results are brought into total magnetic resistance R of q axisqCalculating formula (11) in:
Rq=(R11+…+Rk1+…)//(R′11+…+R′k1+…) (11)
It enables, k=8 is to get to total magnetic resistance R of the q axis of this Fa-SynRMq。
Wherein, 81 represent in Fig. 8 that finally mark is Ry9, and it ' is then that can become 8 series connection magnetic after △-Y ' to pass through
It hinders (referring to Fig. 9), 81 subsequent 1 be due to ' generally being written as being R by first variable of transformation out in △-Y ' transformation11, this
Place is consistent, without physical meaning.
By formula (7) to formula (10) it is found that the key for solving q axis magnetic resistance is to be that solving magnetic flux passes perpendicularly through i-th layer of magnetic barrier
The axial area for hindering (or ferrite), due to radial reinforcing rib TrThe area very little of region seeks arc using two circle intersections
Long xiMethod approximate calculation, and again due to rotor structure symmetry, meet thus only need to analyze using d axis as symmetry axis wherein
Half of magnetic pole.
Figure 10 show one kind that the present invention improves the optimization method of ferrite assist type synchronous magnetic resistance motor power factor
Embodiment using power factor as the genetic algorithm optimization result of object function.
Figure 11 show one kind that the present invention improves the optimization method of ferrite assist type synchronous magnetic resistance motor power factor
Four layers of ferrite width of embodiment and power factor graph of relation.Abscissa is four layers of ferritic width value hmi, indulge and sit
It is designated as the power factor of Fa-SynRM.
Figure 12 show one kind that the present invention improves the optimization method of ferrite assist type synchronous magnetic resistance motor power factor
Embodiment split than coefficient and Fa-SynRM power because and efficiency graph of relation.Abscissa is main vertical to split than coefficient value
Coordinate is the power factor of Fa-SynRM, and secondary ordinate is the efficiency of Fa-SynRM.
Embodiment
The present embodiment is using SynRM as objective for implementation, followed by using Fa-SynRM as objective for implementation first.Two kinds of motors
Structural material is completely the same, and rotor structure is alternately folded by cold-reduced silicon sheet and ambroin (or the insulating materials such as electrician's paper)
It presses, specifically Fa-SynRM has inserted Ferrite Material appropriate (Ferrite Material, which is one kind, to be had in insulation
Ferromagnetic metal oxide, resistivity ratio metal, alloy magnetic material are much bigger, and also higher dielectric properties, high
There is higher magnetic conductivity, since the magnetic energy stored in ferrite unit volume is relatively low, the conjunction intensity of magnetization of satisfying is relatively low, leads to when frequency
Often there was only the 1/3~1/5 of pure iron), cold-reduced silicon sheet is a kind of magnetic conductive material, and physical characteristic shows as anisotropy, i.e.,
Along the different directions of lattice, the periodicity and density degree of atomic arrangement are not quite similar, and thus cause crystal in different directions
Physicochemical characteristics is also different, therefore, determines that the magnetic line of force for flowing through d axis magnetic circuits shows as paramagnetic, and flows through the magnetic of q axis magnetic circuits
The line of force shows as inverse magnetic.
The method and step for improving the present embodiment ferrite assist type synchronous magnetic resistance motor power factor is as follows:
The first step:Determine the structural parameters of SynRM
1-1, it needs to be determined that SynRM elementary structure parameter
I, measurements obtain the stator structure parameter of SynRM, including number of stator slots ns, stator rabbet width Sso, stator slot away from
Angle beta;II, measurements obtain the Structural Parameters of its Rotor of SynRM, including rotor number of pole-pairs P, rotor layers of flux barriers number nl, rotor slot opening
Spend Rso, iron core effective length l, the gas length δ between diameter of stator bore and rotor diameter;III. the stator structure ginseng of SynRM is determined
Number, the stator structure of SynRM directly use the stator of the threephase asynchronous of model Y160-4M, major parameter:Outer diameter
260mm, internal diameter 170mm, parallel branch number are 1, and coil turn 28, iron core effective length l are 155mm;Specific other sizes
Parameter can refer to that Zhao Bo, Zhang Hongliang are written, be published within 2010 by Chinese Water Conservancy water power publishing house《Ansoft12 is in engineering electricity
Application in magnetic field》One book;
The establishment of 1-2, rotor slot-pitch angle:
To obtain low torque ripple characteristic, need to determine rotor slot number according to given number of stator slots, such as formula (1)
It is shown:
Determine that rotor slot number is nr, magnetic hinders both ends alignment, and is correspondingly aligned with stator slot on the pitch angle of rotor slot;
Shown in the mathematical expression of rotor slot-pitch angle such as formula (2):
Second step:Taguchi optimizes SynRM rotor insulation magnetic and hinders structure:
Define Structural Parameters of its Rotor:Gas length δ between diameter of stator bore and rotor diameter and every layer of magnetic barrier width hi, composition
Influence the controllable parameter of motor torque ripple;Gas length level range is provided first, then at each δ, changes every layer of magnetic
Hinder width to compare the size of torque pulsation, finds the level range of every layer of insulation magnetic barrier width Optimal Parameters respectively, secondly exist
Gas between each diameter of stator bore and rotor diameter respectively within the scope of the level of gas length δ between diameter of stator bore and rotor diameter
Magnetic barrier width h under gap length δiTaguchi optimization designs are carried out, the torque pulsation value under different δ is found, obtains Taguchi
The gas length δ between diameter of stator bore and rotor diameter after optimization;It recycles outside the diameter of stator bore and rotor after Taguchi optimizations
Gas length δ between diameter, in every layer of magnetic barrier width hiRespectively to every layer of magnetic barrier width h within the scope of the level of Optimal ParametersiIt carries out
Taguchi optimization designs find different hiUnder torque pulsation value, obtain every layer of magnetic barrier width h after Taguchi optimizationsi;Most
Every layer of magnetic barrier width h after corresponding Taguchi optimizations is determined afterwardsiUnder rotor insulation magnetic hinder end spacing tj(j=1,
2,···,nl+ 1) and rotor magnetic hinders following arc length xi;
Table 1 is with Y160-4M threephase asynchronouses with elementary structure parameter in the SynRM of capacity, nsFor number of stator slots;
SsoFor stator rabbet width;β is stator slot-pitch angle;α is rotor slot-pitch angle;nlFor rotor layers of flux barriers number;nrFor rotor slot number;Rso
For rotor width of rebate;L is iron core effective length;Gas lengths of the δ between diameter of stator bore and rotor diameter;hiFor adjacent two layers
Insulation magnetic between magnetic conduction silicon steel sheet hinders width;tjHinder end spacing for rotor insulation magnetic;xiHinder following arc length for rotor magnetic;
Table 2 is that Taguchi optimizes SynRM rotor insulation magnetic barrier structure results value;
Third walks:The rotor structure of Fa-SynRM optimizes:
The total magnetic resistance R of 3-1, d, q axis magnetic circuitd、RqEquivalent Calculation:
It hinders to meet parallel relationship between each magnetic resistance of each layer magnetic conduction silicon steel sheet magnetic flux circulation, passes through formula (3) equivalent meter
Calculate the total magnetic resistance R of d axis magnetic circuitsd:
Rd=Rd1//Rd2//…//Rdj//… (3)
Wherein, j is expressed as jth layer magnetic conduction silicon steel sheet, RdjIndicate the air gap magnetic for hindering jth layer magnetic conduction silicon steel sheet magnetic flux to circulate
Resistance value, RdjIt is obtained by formula (4):
In formula, μ0For space permeability, value is 4 π × 10-7H/m;
To simplify formula, the Equivalent Calculation value of d, q axis magnetic circuit resistance is indicated with perunit value, wherein a reference value is formula (5):
In formula, μrFor ferritic relative permeability;
That is, formula (4) is expressed as formula (6) with perunit value:
Magnetic flux is hindered to carry over rotor magnetic conduction silicon steel sheet across the air gap between diameter of stator bore and rotor diameter after ferrite is added
Magnetic resistance RyjFormula (7) can be used to indicate;It hinders across the ferritic magnetic resistance R of each layermiFormula (8) can be used to indicate;It hinders exhausted across each layer
The magnetic resistance R of edge barrierδiIt is indicated according to formula (9);Similarly, another part that q axis is symmetrically distributed hinders magnetic flux to carry over rotor magnetic conduction
Silicon steel sheet passes through the magnetic resistance R ' of the air gap between diameter of stator bore and rotor diameteryjIt is calculated according to formula (7);It hinders across each layer iron
The magnetic resistance R ' of oxysomemiWith obstruction across the magnetic resistance R ' of each layer insulation barrierδiIt is indicated respectively according to formula (8) and formula (9):
Wherein, k1It is magnetic flux area correction value, general value existsBetween, hmiFor ferrite width, bmiFor iron oxygen
Body length;
Using ' relationship that △-Y ' become the magnetic resistance for equivalent q axis magnetic circuit of changing commanders is simplified, then through ' △-Y ' transformation can obtain
Ry(j+1)、Ryj、Rδi、RmiBetween meet formula (10):
In formula, the value of k is that (wherein 1 arrives n to the corresponding layers of flux barriers number of q axis magnetic circuitslWith nlTo 2nlCorrespondent equal), i.e. Rk1It is
About the magnetic resistance R for hindering magnetic flux to carry over the air gap that rotor magnetic conduction silicon steel sheet passes through between diameter of stator bore and rotor diameteryjIt is worn with obstruction
Cross the ferritic magnetic resistance R of each layermiWith obstruction across the magnetic resistance R of each layer barrierδiFunction, and it is another to be symmetrically distributed with q axis
Part hinders the magnetic resistance R ' that magnetic flux carries over rotor magnetic conduction silicon steel sheet across the air gap between diameter of stator bore and rotor diameteryjIt is worn with obstruction
Cross the ferritic magnetic resistance R ' of each layermiWith obstruction across the magnetic resistance R ' of each layer insulation barrierδiThe same formula of calculating formula (10);
It is then recycled that ' △-Y ' transformation, can be obtained total magnetic resistance calculating formula of q axis magnetic circuits;With formula (11) Equivalent Calculation
The total magnetic resistance R of q axis magnetic circuitsq:
Rq=(R11+…+Ek1+…)//(R′11+…+R′k1+…) (11)
3-2, power factor and the total magnetic resistance R of d, q axis magnetic circuit are establishedd、RqRelationship:
By total magnetic resistance R of the obtained d axis magnetic circuits of 3-1dWith total magnetic resistance R of q axis magnetic circuitsqIt brings into formula (12), obtains d, q axis
The calculating formula of the power factor PF of equivalent magnetic circuit, the Structural Parameters of its Rotor which establishes power factor PF and be determined above
Rotor insulation magnetic hinders end spacing tj, insulation magnetic barrier width h between adjacent two layers magnetic conduction silicon steel sheeti, rotor magnetic hinder following arc length
xi, ferrite width hmi, ferrite length bmi, relationship between gas length δ between diameter of stator bore and rotor diameter, that is, establish
The relationship of power factor and Fa-SynRM dependency structure parameters:
Wherein:
3-3, genetic algorithm and Taguchi methods optimize Fa-SynRM rotor structures:
Relevant genetic algorithm program is write according to formula (12), Fa-SynRM rotor structures are optimized, mesh
Be to obtain the value range of each variable under power factor PF maximums, then using Taguchi methods to variable carry out part
Optimization, and then obtain ferritic minimum volume value under conditions of meeting power factor PF and being basically unchanged;On this basis, into
One step considers to split the influence than coefficient to Fa-SynRM power factors, obtains optimal electric machine structure parameter, obtain final electricity
Machine Optimal Structure Designing value;Table 3 is the front and back ferrite structure of Fa-SynRM optimizations and motor performance contrast table.
The initial value of each structural parameters to be optimized of SynRM in table 1 is substituted into respectively in the present embodiment and calculates the total magnetic resistance of d axis
In the formula (3) and formula (11) of the total magnetic resistance of q axis, obtain that the total magnetic resistance of q axis under perunit value is 0.0514, the total magnetic resistance of d axis is 0.0166.
So far, the total magnetic resistance of equivalent d axis and the total magnetic resistance of equivalent q axis magnetic circuit are calculated and are terminated.By the total magnetic resistance of obtained d axis and the total magnetic resistance of q axis
It is updated in calculating formula (12), obtains initial power factor of the method for the present invention embodiment after the calculating of d, q axis equivalent magnetic resistance
It is 0.512, and is directly 0.498 using the power factor that FInite Element obtains, torque pulsation 8.31%, the method for the present invention obtains
To initial power factor compared with the power factor obtained by FInite Element, the two differs only by 0.014, as a result connects very much
Closely, this show do not consider room and time harmonic wave, leakage magnetic flux, core loss and it is magnetically saturated under the conditions of, by the method for the present invention
The equivalent magnetic circuit modeling of the SynRM set up and the method for calculating the power factor of SynRM are rationally effective, Ke Yitong
It ought to be used in Fa-SynRM.
Ensure that the structure torque pulsation of SynRM is minimum first, and then suitable ferrite is added in insulation magnetic barrier, leads to
The relationship of the power factor and Fa-SynRM dependency structure parameters established is crossed, it is excellent to the carry out of the dependency structure parameter of Fa-SynRM
Change adjustment, comprises the concrete steps that:
A. the gas length δ and every layer of magnetic barrier width h between SynRM diameter of stator bore and rotor diameter are choseni
The essence of torque pulsation is that stator current and rotor field interact and generate, and passes through and changes SynRM rotors
Magnetic, which hinders structure and gas length, can be effectively improved torque pulsation problem.The magnetic resistance of motor is mainly in the diameter of stator bore of motor and rotor
Between outer diameter, therefore, under conditions of other original dimensions are constant, change the gas length δ between diameter of stator bore and rotor diameter,
D axis magnetic resistance and q axis magnetic resistance can be had an impact.
The level ranging from 0.3mm-0.8mm of gas length δ in the present embodiment between diameter of stator bore and rotor diameter, from figure
Comparison can obtain in 3, and with the variation of gas length, torque pulsation and power factor also change therewith, and show certain rule
Rule, when the gas length δ between diameter of stator bore and rotor diameter is 0.5mm, torque pulsation is minimum, and the results are shown in Table 2, power
Though factor is not maximum, it is not much different.Gas length δ between diameter of stator bore and rotor diameter is bigger, the torque pulsation of motor
First reduce and increases afterwards.Therefore, shown by experimental data the diameter of stator bore of SynRM rotor structures that the method for the present invention optimizes with
Gas length δ between rotor diameter is 0.5mm, and on this basis to every layer of magnetic barrier width hiDegree carries out Taguchi optimizations and sets
Meter.
B. ferrite size is chosen
In conjunction with above-mentioned obtained experimental data, it is 0.5mm to take the gas length δ between diameter of stator bore and rotor diameter, other knots
Structure parameter is still identical as original dimension, and here, d, q axis magnetic resistance change rate in equivalent magnetic circuit are all made of per unit value expression.
The width h of four layers of magnetic barrieriIt is 4.8mm, 6.3mm, 8.5mm, 11mm respectively, is inserted into ferrite at its center, leads to
It crosses genetic algorithm and show that four layers of ferrite length are essentially the same, be 22.8mm.Fa-SynRM is found out by the curve in Figure 11 again
Power factor changes with the variation of ferrite width in insulation magnetic barrier, every curve all there are one relatively large peak value,
Can obtain every layer of ferritic length difference by emulation influences less, so ensureing every layer of ferrite herein power factor
Length it is constant and minimum as possible, therefore, increases of ferrite width means the increase of cost, so in consideration power factor
While promotion, cost is reduced as far as possible.Due to being limited by the size of rotor diameter and mechanical stress, ferrite is wide
Degree can not possibly infinitely increase.Therefore, each of Fa-SynRM rotor structures that the method for the present invention optimizes is shown by experimental data
Layer ferrite width is 3.4mm.Ferrite is equal sized in each layer insulation magnetic barrier of Fa-SynRM rotor structures.
C. searching is best splits than coefficient, further increases power factor
According to above-mentioned experimental result, the ferrite size that each layer insulation magnetic barrier is inserted into is long 22.8mm, width 3.4mm, stator
Gas length δ between internal diameter and rotor diameter is 0.5mm, under conditions of keeping above-mentioned other structures parameter constant, only changes electricity
Machine is split than coefficient, the present embodiment change first the method split and taken than the size of coefficient be compared to it is initial split done than coefficient it is suitable
When increase and reduction, that is, split and change within the scope of 0.5-0.8 than coefficient, as seen in Figure 12, in a certain range,
Power factor is to increase with splitting than the increase of coefficient, but corresponding efficiency value can reduce, so declining not in guaranteed efficiency
Under the premise of big, increase power factor value than coefficient by increasing to split, optimal split is than coefficient value as seen from Figure 12
0.7。
The method that the present invention improves ferrite assist type synchronous magnetic resistance motor power factor show the power of Fa-SynRM because
The magnetic resistance of d, q axis of several sizes and equivalent magnetic circuit has relationship, and then related with insulation magnetic barrier and ferrite size.Herein
In relationship, by the optimization to parameters, power factor is stepped up, finally in splitting than coefficient more one by change motor
The raising power factor of step, to obtain final motor model.
The practical initial value of structural parameters basic 1 SynRM of table
2 SynRM of table utilizes the result table of Taguchi method optimization designs δ
δ(mm) | h1(mm) | h2(mm) | h3(mm) | h4(mm) | Torque pulsation (%) | Power factor |
0.3 | 3.6 | 4.7 | 5 | 9 | 11.50 | 0.441 |
0.4 | 4.7 | 5.3 | 8 | 9.8 | 9.53 | 0.506 |
0.5 | 4.8 | 6.3 | 8.5 | 11 | 8.31 | 0.498 |
0.6 | 4.9 | 6.2 | 8.7 | 11.4 | 9.39 | 0.472 |
0.7 | 4.2 | 6 | 8.2 | 11.3 | 10.31 | 0.433 |
0.8 | 3.8 | 5.3 | 7.3 | 10.5 | 12.04 | 0.412 |
3 Fa-SynRM of table optimizes comparing result table
Before optimization | After optimization | |
b1 | 22.67 | 22.8 |
b2 | 22.73 | 22.8 |
b3 | 23.14 | 22.8 |
b4 | 22.95 | 22.8 |
h1 | 1.25 | 3.6 |
h2 | 3.68 | 3.4 |
h3 | 6.85 | 3.4 |
h4 | 9.55 | 3.4 |
Power factor | 0.724 | 0.826 |
Torque pulsation (%) | 14.31 | 10.87 |
Efficiency (%) | 78.27 | 92.91 |
The present invention does not address place and is suitable for the prior art.
Claims (2)
1. a kind of optimization method improving ferrite assist type synchronous magnetic resistance motor power factor, this method comprise the concrete steps that:
The first step:Determine the structural parameters of synchronous magnetic resistance motor (SynRM):
Measurement obtains the stator structure parameter of SynRM, including number of stator slots ns, stator rabbet width Sso, stator slot-pitch angle β;Its
Secondary, measurement obtains the Structural Parameters of its Rotor of SynRM, including rotor number of pole-pairs P, rotor layers of flux barriers number nl, rotor width of rebate Rso、
Gas length δ between iron core effective length l, diameter of stator bore and rotor diameter;
The number of stator slots obtained according to above-mentioned measurement is needed according to formula (1) to determine rotor slot number nr,
Determining that rotor slot number is nrAfterwards, rotor slot-pitch angle α is calculated according to formula (2),
Second step:Taguchi optimizes SynRM rotor insulation magnetic and hinders structure:
Gas length δ between diameter of stator bore and rotor diameter and every layer of magnetic barrier width hi, optimize, i=1,2 ..., nl:It is first
The level range for first providing the gas length δ between diameter of stator bore and rotor diameter changes every layer of magnetic barrier then under each level
Width compares the size of torque pulsation, respectively finds every layer of magnetic barrier width hiThe level range of Optimal Parameters;Secondly in stator
Air gap within the scope of the level of gas length δ between internal diameter and rotor diameter respectively between each diameter of stator bore and rotor diameter is long
Spend the magnetic barrier width h under δiTaguchi optimization designs are carried out, the torque pulsation value under different δ is found, obtain Taguchi optimizations
The gas length δ between diameter of stator bore and rotor diameter afterwards;Between diameter of stator bore and rotor diameter after recycling Taguchi to optimize
Gas length δ, in every layer of magnetic barrier width hiRespectively to every layer of magnetic barrier width h within the scope of the level of Optimal ParametersiIt carries out
Taguchi optimization designs find different hiUnder torque pulsation value, obtain every layer of magnetic barrier width h after Taguchi optimizationsi;Most
Every layer of magnetic barrier width h after corresponding Taguchi optimizations is determined afterwardsiUnder rotor insulation magnetic hinder end spacing tjAnd rotor magnetic barrier
Following arc length xi, j=1,2, nl+1;
Third walks:The rotor structure of ferrite assist type synchronous magnetic resistance motor (Fa-SynRM) optimizes:
The total magnetic resistance R of 3-1, d, q axis magnetic circuitd、RqEquivalent Calculation:
It hinders to meet parallel relationship between each magnetic resistance of each layer magnetic conduction silicon steel sheet magnetic flux circulation, passes through formula (3) Equivalent Calculation d axis
The total magnetic resistance R of magnetic circuitd:
Rd=Rd1//Rd2//…//Rdj//… (3)
Wherein, j indicates jth layer magnetic conduction silicon steel sheet, j=1,2 ..., nl+1;RdjIt indicates to hinder jth layer magnetic conduction silicon steel sheet magnetic flux stream
Logical air-gap reluctance value, RdjIt is obtained by formula (6):
In formula, μrFor ferritic relative permeability, δ is between diameter of stator bore and rotor diameter after second step Taguchi optimizations
Gas length;
Recycle that ' relationship that △-Y ' become the magnetic resistance for equivalent q axis magnetic circuit of changing commanders is simplified, with formula (11) Equivalent Calculation q axis magnetic
The total magnetic resistance R in roadq:
Rq=(Rll+…+Rkl+…)//(R′ll+…+R′kl+…) (11)
In formula, k is the corresponding layers of flux barriers number of q axis magnetic circuits, wherein 1 arrives nlWith nlTo 2nlCorrespondent equal, k=2nl;Rk1It is about adding
Enter the magnetic resistance R for hindering magnetic flux to carry over the air gap that rotor magnetic conduction silicon steel sheet passes through between diameter of stator bore and rotor diameter after ferriteyjWith
It hinders across the ferritic magnetic resistance R of each layermiWith obstruction across the magnetic resistance R of each layer barrierδiFunction, with formula (10) indicate;
RyjIt is indicated with formula (7), RmiIt is indicated with formula (8), RδiIt is indicated with formula (9);Similarly, another part that q axis is symmetrically distributed
R′yj, R 'mi, R 'δiFormula (7), (8), (9) are used to indicate respectively;
Wherein, k1It is magnetic flux area correction value, hmiFor ferrite width, bmiFor ferrite length;
3-2, with the total magnetic resistance R of d, q axis magnetic circuitd、RqIndicate power factor:
The total magnetic resistance R of d axis magnetic circuits that step 3-1 is obtaineddWith the total magnetic resistance R of q axis magnetic circuitsqIt brings into formula (12), obtains by d, q axis etc.
Imitate the power factor calculating formula represented by magnetic resistance, the Structural Parameters of its Rotor which establishes power factor PF and be determined above
Rotor insulation magnetic hinders end spacing tj, insulation magnetic barrier width h between adjacent two layers magnetic conduction silicon steel sheeti, rotor magnetic hinder following arc length
xi, ferrite width hmi, ferrite length bmi, relationship between gas length δ between diameter of stator bore and rotor diameter, that is, establish
The relationship of power factor and Fa-SynRM dependency structure parameters:
3-3, genetic algorithm and Taguchi methods optimize Fa-SynRM rotor structures:
According to formula (12) by genetic algorithm, Fa-SynRM rotor structures are optimized, obtain power factor PF most
Then the value of each variable under greatly carries out local optimum to variable using Taguchi methods, and then is meeting power factor PF not
Ferritic minimum volume value is obtained under conditions of change;On this basis, further consider to split than coefficient to Fa-SynRM power
The influence of factor further increases the power factor of motor.
2. the optimization method according to claim 1 for improving ferrite assist type synchronous magnetic resistance motor power factor, special
Sign is that the Fa-SynRM selects lateral lamination armature motor.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810361597.0A CN108494186B (en) | 2018-04-20 | 2018-04-20 | A kind of optimization method improving ferrite assist type synchronous magnetic resistance motor power factor |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810361597.0A CN108494186B (en) | 2018-04-20 | 2018-04-20 | A kind of optimization method improving ferrite assist type synchronous magnetic resistance motor power factor |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108494186A true CN108494186A (en) | 2018-09-04 |
CN108494186B CN108494186B (en) | 2019-06-18 |
Family
ID=63313867
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810361597.0A Active CN108494186B (en) | 2018-04-20 | 2018-04-20 | A kind of optimization method improving ferrite assist type synchronous magnetic resistance motor power factor |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108494186B (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110556979A (en) * | 2019-09-05 | 2019-12-10 | 合肥工业大学 | Method for optimizing split ratio and magnetization angle of Halbach array permanent magnet motor |
CN110739821A (en) * | 2019-11-06 | 2020-01-31 | 天津工业大学 | Method for designing robustness of low-iron-loss variable-flux permanent magnet memory motor for electric automobile |
CN110851949A (en) * | 2019-09-01 | 2020-02-28 | 天津工业大学 | Method for analyzing electromagnetic performance of multilayer magnetic barrier permanent magnet auxiliary type synchronous reluctance motor |
CN111859574A (en) * | 2020-07-22 | 2020-10-30 | 合肥工业大学 | Synchronous reluctance motor rotor optimization design method for reducing torque ripple |
CN112541154A (en) * | 2020-11-26 | 2021-03-23 | 东南大学 | Method for calculating magnetic circuit power |
CN113098170A (en) * | 2021-03-31 | 2021-07-09 | 合肥工业大学 | Optimization method of built-in permanent magnet motor air gap field based on Taguchi method |
CN113708520A (en) * | 2021-09-18 | 2021-11-26 | 广东美芝制冷设备有限公司 | Motor, compressor and refrigeration plant |
CN114142688A (en) * | 2021-12-02 | 2022-03-04 | 常州工学院 | Reluctance motor magnetic barrier end deflection method for inhibiting torque pulsation |
CN114997012A (en) * | 2022-06-14 | 2022-09-02 | 福州大学 | Ferrite magnetic field optimization device and method based on genetic algorithm |
CN116670973A (en) * | 2020-11-05 | 2023-08-29 | 伊恩运动美国有限责任公司 | Generator for driving a predefined load |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101425726A (en) * | 2008-11-18 | 2009-05-06 | 天津大学 | Motor optimized design method based on fuzzy expert system multi-target particle team |
CN104901452A (en) * | 2015-05-12 | 2015-09-09 | 上海吉亿电机有限公司 | Permanent magnet assisted synchronous reluctance motor rotor which can be used in high-speed situation |
CN106169896A (en) * | 2016-07-13 | 2016-11-30 | 阳光电源股份有限公司 | The controller of magneto alternator and rotor flux on-line amending method and apparatus |
CN106528948A (en) * | 2016-10-19 | 2017-03-22 | 江苏航天动力机电有限公司 | Teaching and learning algorithm-based calculation method for permanent magnet synchronous motor design |
-
2018
- 2018-04-20 CN CN201810361597.0A patent/CN108494186B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101425726A (en) * | 2008-11-18 | 2009-05-06 | 天津大学 | Motor optimized design method based on fuzzy expert system multi-target particle team |
CN104901452A (en) * | 2015-05-12 | 2015-09-09 | 上海吉亿电机有限公司 | Permanent magnet assisted synchronous reluctance motor rotor which can be used in high-speed situation |
CN106169896A (en) * | 2016-07-13 | 2016-11-30 | 阳光电源股份有限公司 | The controller of magneto alternator and rotor flux on-line amending method and apparatus |
CN106528948A (en) * | 2016-10-19 | 2017-03-22 | 江苏航天动力机电有限公司 | Teaching and learning algorithm-based calculation method for permanent magnet synchronous motor design |
Non-Patent Citations (1)
Title |
---|
蔡顺: "同步磁阻电机性能分析与结构优化", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110851949B (en) * | 2019-09-01 | 2023-08-25 | 天津工业大学 | Method for analyzing electromagnetic performance of multilayer magnetic barrier permanent magnet auxiliary synchronous reluctance motor |
CN110851949A (en) * | 2019-09-01 | 2020-02-28 | 天津工业大学 | Method for analyzing electromagnetic performance of multilayer magnetic barrier permanent magnet auxiliary type synchronous reluctance motor |
CN110556979A (en) * | 2019-09-05 | 2019-12-10 | 合肥工业大学 | Method for optimizing split ratio and magnetization angle of Halbach array permanent magnet motor |
CN110739821A (en) * | 2019-11-06 | 2020-01-31 | 天津工业大学 | Method for designing robustness of low-iron-loss variable-flux permanent magnet memory motor for electric automobile |
CN110739821B (en) * | 2019-11-06 | 2024-04-30 | 天津工业大学 | Method for designing robustness of low-iron-loss variable-flux permanent magnet memory motor for electric automobile |
CN111859574A (en) * | 2020-07-22 | 2020-10-30 | 合肥工业大学 | Synchronous reluctance motor rotor optimization design method for reducing torque ripple |
CN111859574B (en) * | 2020-07-22 | 2024-05-03 | 合肥工业大学 | Synchronous reluctance motor rotor optimization design method for reducing torque pulsation |
CN116670973A (en) * | 2020-11-05 | 2023-08-29 | 伊恩运动美国有限责任公司 | Generator for driving a predefined load |
CN112541154A (en) * | 2020-11-26 | 2021-03-23 | 东南大学 | Method for calculating magnetic circuit power |
WO2022110528A1 (en) * | 2020-11-26 | 2022-06-02 | 东南大学 | Method for calculating power of magnetic circuit |
US11709211B2 (en) | 2020-11-26 | 2023-07-25 | Southeast University | Power calculation method of magnetic circuit |
CN112541154B (en) * | 2020-11-26 | 2021-10-08 | 东南大学 | Method for calculating magnetic circuit power |
CN113098170A (en) * | 2021-03-31 | 2021-07-09 | 合肥工业大学 | Optimization method of built-in permanent magnet motor air gap field based on Taguchi method |
CN113708520A (en) * | 2021-09-18 | 2021-11-26 | 广东美芝制冷设备有限公司 | Motor, compressor and refrigeration plant |
CN114142688A (en) * | 2021-12-02 | 2022-03-04 | 常州工学院 | Reluctance motor magnetic barrier end deflection method for inhibiting torque pulsation |
CN114997012A (en) * | 2022-06-14 | 2022-09-02 | 福州大学 | Ferrite magnetic field optimization device and method based on genetic algorithm |
Also Published As
Publication number | Publication date |
---|---|
CN108494186B (en) | 2019-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108494186B (en) | A kind of optimization method improving ferrite assist type synchronous magnetic resistance motor power factor | |
Dhulipati et al. | Torque performance enhancement in consequent pole PMSM based on magnet pole shape optimization for direct-drive EV | |
Zhang et al. | A hybrid analytical model for open-circuit field calculation of multilayer interior permanent magnet machines | |
Zhao et al. | Optimal design of a spoke-type permanent magnet motor with phase-group concentrated-coil windings to minimize torque pulsations | |
CN107786062A (en) | A kind of segmented cylinder-type transverse-flux linear motor | |
Zhao et al. | Optimal design and experimental test of a SPM motor with cost-effective magnet utilization to suppress torque pulsations | |
Wu et al. | Axial‐flux permanent‐magnet synchronous generator with coreless armature and non‐integral coil–pole ratio | |
CN104656016A (en) | Method for analyzing stable-state performance of non-sine power-supply multiphase induction motor | |
CN109768645B (en) | Permanent magnet auxiliary type synchronous reluctance motor rotor structure and design method thereof | |
CN106295004B (en) | Permanent magnet motor optimization design method considering disturbance type design variable interval | |
Ullah et al. | Analysis of a discrete stator hybrid excited flux switching linear machine | |
Chen et al. | Multi-objective optimum design of five-phase squirrel cage induction motor by differential evolution algorithm | |
CN105356804B (en) | Improve synchronous magnetic resistance motor salient pole than method | |
Yu et al. | Analytical model for magnetic field calculation of SPMSM with chamfered pole considering iron core saturation | |
Xu et al. | Torque ripple improvement for ferrite-assisted synchronous reluctance motor by using asymmetric flux-barrier arrangement | |
Xu et al. | Optimum design of interior permanent magnet synchronous motor using Taguchi method | |
CN110460185A (en) | Rotor and motor | |
CN107786061B (en) | A kind of plate transverse flux linear motor | |
CN108733963A (en) | A kind of threephase asynchronous permanent-magnetic energy-saving transformation rotor design method | |
Kazemi Sangdehi et al. | Analysis and optimization of a 12/14 double‐stator flux‐switching machine using low cost magnet | |
CN115842496A (en) | High-precision magnetic network modeling method suitable for hybrid stator permanent magnet vernier motor | |
Bao et al. | A novel computational method of skewing leakage reactance for a doubly skewed rotor induction motor | |
Yang et al. | Optimization of squirrel-cage rotor for amorphous asynchronous motor | |
Deng et al. | Design and finite element analysis of a novel permanent magnet assisted reluctance synchronous motor | |
Zhang et al. | Characteristic analysis of conventional pole and consequent pole IPMSM for electric vehicle application |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |