WO2023103138A1 - Method for collaborative optimization design of permanent magnet-armature double harmonics of magnetic field-modulated permanent magnet motor - Google Patents

Method for collaborative optimization design of permanent magnet-armature double harmonics of magnetic field-modulated permanent magnet motor Download PDF

Info

Publication number
WO2023103138A1
WO2023103138A1 PCT/CN2022/070687 CN2022070687W WO2023103138A1 WO 2023103138 A1 WO2023103138 A1 WO 2023103138A1 CN 2022070687 W CN2022070687 W CN 2022070687W WO 2023103138 A1 WO2023103138 A1 WO 2023103138A1
Authority
WO
WIPO (PCT)
Prior art keywords
magnetic field
harmonics
armature
working
harmonic
Prior art date
Application number
PCT/CN2022/070687
Other languages
French (fr)
Chinese (zh)
Inventor
徐亮
吴文杰
蒋婷婷
赵文祥
Original Assignee
江苏大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 江苏大学 filed Critical 江苏大学
Priority to GB2218908.8A priority Critical patent/GB2612711A/en
Publication of WO2023103138A1 publication Critical patent/WO2023103138A1/en

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K1/00Details of the magnetic circuit
    • H02K1/06Details of the magnetic circuit characterised by the shape, form or construction
    • H02K1/12Stationary parts of the magnetic circuit
    • H02K1/17Stator cores with permanent magnets
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K1/00Details of the magnetic circuit
    • H02K1/06Details of the magnetic circuit characterised by the shape, form or construction
    • H02K1/22Rotating parts of the magnetic circuit
    • H02K1/223Rotor cores with windings and permanent magnets
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K21/00Synchronous motors having permanent magnets; Synchronous generators having permanent magnets
    • H02K21/02Details
    • H02K21/021Means for mechanical adjustment of the excitation flux
    • H02K21/028Means for mechanical adjustment of the excitation flux by modifying the magnetic circuit within the field or the armature, e.g. by using shunts, by adjusting the magnets position, by vectorial combination of field or armature sections
    • H02K21/029Vectorial combination of the fluxes generated by a plurality of field sections or of the voltages induced in a plurality of armature sections
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K7/00Arrangements for handling mechanical energy structurally associated with dynamo-electric machines, e.g. structural association with mechanical driving motors or auxiliary dynamo-electric machines
    • H02K7/04Balancing means
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/60Other road transportation technologies with climate change mitigation effect
    • Y02T10/64Electric machine technologies in electromobility

Definitions

  • the invention relates to a permanent magnet-armature dual-harmonic collaborative optimization design method for a magnetic field modulation permanent magnet motor, which belongs to the field of motor design and is specifically applicable to motors requiring high torque density and high power factor such as electric vehicles, wind power generation, ship propulsion, etc. system.
  • Field modulated permanent magnet motors are characterized by high torque density due to the "field modulation effect".
  • the so-called “magnetic field modulation effect” means that the magnetic field modulation permanent magnet motor modulates and generates a variety of working magnetic field harmonics that can generate torque through the action of its modulation pole, thereby increasing the torque of the motor.
  • Field modulated permanent magnet motors have great development potential in the field of direct drive motors due to their high torque density.
  • the magnetic field modulation permanent magnet motor itself has the problem of high flux leakage, which makes the power factor low and restricts its practical application.
  • the design method to improve the power factor of the magnetic field modulation permanent magnet motor is mainly to increase the permanent magnet flux linkage from the perspective of increasing the permanent magnet magnetic field, or to reduce the armature current from the perspective of reducing the armature magnetic field.
  • the current design method only improves the power factor of the magnetic field modulated permanent magnet motor from a single angle of the permanent magnet magnetic field or the armature magnetic field. From the final result, in some methods, the increase of the power factor will lead to the decrease of the torque density, while some methods Although the power factor of the motor has been improved, the degree of improvement often does not meet expectations, and there is still a large room for improvement.
  • both the permanent magnet magnetic field and the armature magnetic field of the magnetic field modulation permanent magnet motor have an important influence on the power factor of the motor. Ignoring the influence of any magnetic field on the power factor will make the power factor of the magnetic field modulation permanent magnet motor unable to achieve maximum. Therefore, for the magnetic field modulation permanent magnet motor, it is urgent to propose an effective power factor improvement method considering both the permanent magnet and the armature magnetic field.
  • the content of the present invention is to set up the expression between the flux linkage and the magnetic field harmonic according to the harmonic characteristics of the permanent magnet of the magnetic field modulation permanent magnet motor and the armature magnetic field; Expressions for permanent magnet and armature flux linkages to obtain expressions for power factor with respect to permanent magnet and armature field harmonics. According to the vector diagram and the power factor expression, the influence of the flux linkage corresponding to the permanent magnet and armature magnetic field harmonics on the torque and power factor is analyzed, and the permanent magnet-armature double harmonic collaborative optimization design idea is established.
  • Step 1 According to the permanent magnet magnetic field harmonic characteristic formula and the armature magnetic field harmonic characteristic formula, establish the expression between the permanent magnet flux linkage with respect to the permanent magnet magnetic field harmonics, and establish the relationship between the armature flux linkage with respect to the armature magnetic field harmonics According to the flux vector diagram of the magnetic field modulation permanent magnet motor, the expression of the power factor with respect to the permanent magnet and armature flux linkage is established, and the expression of the power factor with respect to the harmonics of the permanent magnet and armature magnetic field is obtained.
  • the permanent magnet and armature flux linkage is established. Magneto-armature double harmonic synergistic optimization design idea to improve the torque density and power factor of the motor;
  • the permanent magnetic field harmonics are divided into two types: permanent magnetic field working harmonics and permanent magnetic field non-working harmonics.
  • the harmonics of the permanent magnetic field that contribute to the torque are the working harmonics of the permanent magnetic field, and the harmonics of the permanent magnetic field that do not contribute to the torque are the non-working harmonics of the permanent magnetic field.
  • the working harmonics of the synthetic permanent magnetic field can be obtained by weighting the working harmonics of the permanent magnetic field.
  • the working harmonic amplitude of the magnetic field can maintain a high torque density while improving the power factor; for the armature magnetic field harmonics, the armature magnetic field
  • the armature magnetic field There are two types of armature magnetic field working harmonics and armature magnetic field non-working harmonics.
  • the armature magnetic field harmonics that contribute to the torque are the armature magnetic field working harmonics, and the armature magnetic field harmonics that do not contribute to the torque are the armature magnetic field non-working harmonics.
  • the flux linkage corresponding to the non-working harmonic of the armature magnetic field is inversely proportional to the power factor, and according to the flux linkage phasor diagram, reducing the non-working harmonic of the armature magnetic field will not affect the torque corresponding to The size of the working area, so reducing the non-working harmonics of the armature field can improve the power factor without losing torque density.
  • the minimum value of the non-working harmonic of the armature magnetic field and the working harmonic of the armature magnetic field are determined.
  • the optimization of the maximum value can realize the collaborative optimization design of the double harmonics of the permanent magnet-armature magnetic field, thereby improving the torque density and power factor of the motor.
  • Step 2 Limit the working harmonic amplitude of the synthetic permanent magnetic field from the perspective of the permanent magnetic field. Taking the minimum value of the working harmonic amplitude of the synthetic permanent magnetic field as a constraint condition, the design parameters that have a greater impact on the working harmonics of the permanent magnetic field are selected through sensitivity analysis.
  • Step 3 Simplify the armature field harmonics optimization objective. Since there are many harmonic orders of the armature magnetic field to be optimized, it is necessary to simplify the optimization objectives of the armature magnetic field harmonics and reduce the number of optimization objectives.
  • First analyze the sensitivity of the armature magnetic field harmonics to the performance of the motor, select the more sensitive armature magnetic field harmonics as the optimization target, and further calculate the distribution of the armature magnetic field harmonics according to the experimental design method. According to the distribution diagram of wave experiment points, the non-operating harmonics of the armature magnetic field with different changing trends of the working harmonics of the armature magnetic field are selected as the harmonics to be optimized, and the simplified optimization target of the armature magnetic field harmonics is obtained by weighting.
  • Step 4 Independent judgment of non-working harmonics of armature magnetic field.
  • By calculating the interaction effect between the armature magnetic field harmonics it is analyzed and judged whether there are relatively independent non-working harmonics in the armature magnetic field harmonics. If there are relatively independent non-working harmonics of the armature magnetic field, go to steps 5.1 and 5.2; if there are no relatively independent non-working harmonics of the armature magnetic field, go to step 5.3
  • Step 5.1 If there are non-operating harmonics of the armature field with relative independence, optimize the non-operating harmonics of the armature field with relative independence separately from the rest of the harmonics of the armature field to reduce the design parameters and optimization objectives dimension to improve the accuracy of the armature field harmonic optimization results.
  • optimize the non-working harmonics of the armature magnetic field with relative independence use sensitivity analysis to select the design parameters that are more sensitive to the non-working harmonics of the armature magnetic field with relative independence, and establish the design parameters with relative Independent Kriging model of non-working harmonics of armature magnetic field. According to the established Kriging model, the optimal design point of non-working harmonics of armature magnetic field with relative independence is selected.
  • Step 5.2 Optimize the working harmonics of the armature magnetic field and other non-working harmonics of the armature magnetic field that do not have relative independence. Taking the design parameter range limited by the working harmonic amplitude of the synthetic permanent magnet magnetic field as the constraint condition, and taking the working harmonic of the armature magnetic field and the simplified non-working harmonic of the armature magnetic field as the optimization objectives, a multi-objective genetic algorithm is adopted Optimize the working harmonics and non-working harmonics of the motor armature magnetic field, and finally determine the motor design scheme with the optimal armature harmonics on the basis of the high synthetic permanent magnet working harmonic amplitude, and realize the permanent magnet-armature The collaborative optimization design of the double harmonics of the magnetic field improves the torque density and power factor of the motor.
  • Step 5.3 Optimize the working harmonics and non-working harmonics of the armature magnetic field. Taking the parameter range limited by the working harmonic amplitude of the synthesized permanent magnet magnetic field as the constraint condition, taking the working harmonic of the armature magnetic field and the simplified non-working harmonic of the armature magnetic field as the optimization target, the multi-objective genetic algorithm is used to optimize the The working harmonics and non-working harmonics of the armature magnetic field of the machine are optimized, and finally the design scheme of the motor with the optimal armature harmonics is determined on the basis of the high synthetic permanent magnet working harmonic amplitude. The coordinated optimization design of harmonics can improve the torque density and power factor of the motor.
  • C m is the Fourier coefficient of the permanent magnetomotive force
  • D i and D j are the Fourier coefficients of the armature magnetomotive force
  • m is the order of the permanent magnetomotive force
  • k is the order of the permeability
  • i and j are the order of the magnetomotive force of the armature
  • P r is the number of pole pairs of the permanent magnet
  • ⁇ r is the mechanical speed of the motor
  • t is the time
  • ⁇ 0 and ⁇ k are the Fourier coefficients of the air gap permeance
  • N s is the electrical Number of pivot slots.
  • the harmonic order of the permanent magnet magnetic field is mP r , mP r ⁇ kN s
  • the harmonic order of the armature magnetic field is i, j, i ⁇ kN s , j ⁇ kN s
  • F m ( ⁇ ,t) is the permanent magnet magnetomotive force expression
  • ⁇ s ( ⁇ ) is the air gap permeance expression, which can be expressed as:
  • r g is the air gap radius
  • l ef is the axial length
  • n c is the number of turns of the winding.
  • step 1 the expression of the armature flux linkage with respect to the harmonics of the armature magnetic field in step 1 is:
  • step 1 the expression of the power factor with respect to the permanent magnet and armature flux linkage in step 1 is:
  • the corresponding flux linkage for the working harmonic of the armature magnetic field is the flux linkage corresponding to the non-working harmonic of the armature magnetic field
  • U and ⁇ r are the phase voltage and frequency respectively
  • E 0 is the back EMF of the permanent magnet.
  • step 2 the expression of the sensitivity calculation formula in step 2 is:
  • Y(x) represents the harmonic amplitude of the permanent magnet and armature magnetic field under different design parameters
  • N is the number of samples
  • x is the motor design parameter
  • the specific steps of the sensitivity analysis in step 2 are as follows: First, use the central composite design sampling method to sample the second-order factor design points, axis points and zero-level center points that meet the second-order regression rotation criterion, providing a basis for the sensitivity calculation. data support. Then, the sensitivity of different design parameters to the working harmonics of the permanent magnetic field is calculated by the sensitivity formula, and the design parameters with greater sensitivity are selected as the design parameters that limit the amplitude of the synthetic permanent magnetic field working harmonics.
  • H sm is the harmonic amplitude of the synthetic permanent magnet magnetic field of the motor
  • x2 represents the design parameter with high sensitivity to the harmonic of the permanent magnet magnetic field of the motor
  • H mh is the working harmonic amplitude of the hth permanent magnet magnetic field value
  • a h is the weighting coefficient corresponding to the working harmonic of the permanent magnetic field
  • g 1 is the minimum constraint value of the synthetic permanent magnetic field working harmonic.
  • the simplified optimization target of armature magnetic field harmonics in step 3 is obtained from the sensitivity analysis and the distribution of experimental point distribution of armature magnetic field harmonics.
  • the central composite design sampling method is used for sampling to provide data support for the sensitivity calculation.
  • the sensitivity of different armature magnetic field harmonics to power factor is calculated by the sensitivity formula, and the magnetic field harmonic with higher sensitivity is selected as the target of armature harmonic optimization.
  • the experimental points in the composite design sampling method draw the distribution map of the experimental points of the armature magnetic field harmonics, and select the armature with a different trend from the working harmonics of the armature magnetic field according to the changing trend of the different armature magnetic field harmonics in the figure
  • the non-working harmonics of the magnetic field are the harmonics that need to be optimized, and the non-working harmonics of the armature magnetic field that have the same change trend as the working harmonics of the armature magnetic field are not considered as the harmonics that need to be optimized.
  • these harmonics that need to be optimized are linearly weighted as optimization targets, reducing the number of optimization targets for armature magnetic field harmonics, and realizing the simplification of armature magnetic field harmonic optimization targets.
  • step 5.1 the optimized model expression of the non-working harmonics of the armature magnetic field with relative independence in step 5.1 is:
  • x3 is a design parameter that has a great influence on the non-operating harmonic of the armature magnetic field with relative independence after sensitivity analysis
  • Hal is the non-operating harmonic of the armature magnetic field with relative independence for the lth time
  • ⁇ l is the weighting coefficient corresponding to the non-working harmonic of the armature magnetic field with relative independence
  • the minimum value of the optimization model is the optimization target.
  • step 5.2 the optimized model expression of the working harmonics of the armature magnetic field and the simplified non-working harmonics of the armature magnetic field that do not have relative independence is:
  • H aPr is the working harmonic of the armature magnetic field of the motor P r , and its maximum value is set as the optimization target, so as to ensure that the optimized motor has a higher torque.
  • the non-operating harmonics of the armature field are synthesized for the simplified motor, and its minimum value is set as the optimization target.
  • x 4 represents the remaining design parameters after excluding the design parameter x 3 from the total design parameters x 1 , H as is the sth non-working harmonic of the armature magnetic field, ⁇ s is the corresponding s-th non-working armature magnetic field Weighting coefficients for harmonics.
  • f 1 (x 2 ) is a function of the value range of the design parameters after the working harmonics of the synthetic permanent magnet magnetic field are limited.
  • H aPr is the working harmonic of the armature magnetic field of the motor P r , and its maximum value is set as the optimization target, so as to ensure that the optimized motor has a higher torque.
  • the non-operating harmonics of the armature field are synthesized for the simplified motor, and its minimum value is set as the optimization target.
  • x1 is the general design parameter of the motor
  • x2 is the design parameter with high sensitivity to the harmonic of the permanent magnet field of the motor
  • H as is the sth non-working harmonic of the armature magnetic field
  • ⁇ s is the corresponding sth electric Weighting factor for non-operating harmonics of the pivot field.
  • f 1 (x 2 ) is a function of the value range of the design parameters after the working harmonics of the synthetic permanent magnet magnetic field are limited.
  • the present invention sets up the expression between the flux linkage and the magnetic field harmonic according to the harmonic characteristics of the permanent magnet and the armature magnetic field; Expressions for the chain to obtain expressions for the power factor with respect to the harmonics of the permanent and armature fields.
  • the influence of the flux linkage corresponding to the permanent magnet and armature magnetic field harmonics on the torque and power factor is analyzed, and the permanent magnet-armature double harmonic collaborative optimization design idea is established, which will be used for the follow-up from permanent magnet
  • the angle of the magnetic and armature double harmonics gives directions for torque density and power factor improvement.
  • the present invention utilizes experimental point distribution calculation, sensitivity analysis and independent judgment of armature magnetic field non-working harmonics to reduce the dimension of optimization objectives and design parameters, further establish the kriging model, and reduce the calculation amount of motor design. Finally, Based on Kriging model and multi-objective optimization algorithm, the optimal design of magnetic field modulation permanent magnet motor is carried out. Compared with the traditional optimization method, the design method of the invention can significantly improve the optimization efficiency and reduce the optimization time.
  • the present invention uses the synthetic permanent magnet magnetic field working harmonic as a constraint condition to obtain the design parameter range satisfying the constraint condition.
  • the armature magnetic field harmonic is used as an optimization target to optimize it to realize permanent magnet-armature
  • the angle synergistic optimization design of double harmonics determines the motor design scheme for the optimal armature magnetic field harmonics based on the high synthetic permanent magnetic field working harmonics. Compared with the current method of optimizing power factor solely from the perspective of permanent magnet or armature, the design method of the invention can further improve the torque density and power factor of the motor.
  • Fig. 1 is a flow chart of a permanent magnet-armature double harmonic collaborative optimization design method for a magnetic field modulation permanent magnet motor in an embodiment of the present invention
  • Fig. 2 is the flux linkage vector diagram of the magnetic field modulation permanent magnet motor of the present invention
  • Fig. 3 is the topological structure and parameter distribution diagram of the magnetic field modulation permanent magnet motor of the present invention.
  • Fig. 4 is the sensitivity analysis result of the design parameters of the magnetic field modulation permanent magnet motor of the present invention about the working harmonic of the synthetic permanent magnet magnetic field;
  • Fig. 5 (a) is the Kriging model calculation result of the working harmonic of the synthetic permanent magnet magnetic field under the change of w m , h m , of the magnetic field modulation permanent magnet motor of the present invention
  • Fig. 5(b) is the Kriging model calculation result of the working harmonics of the synthetic permanent magnet magnetic field under the change of ws , wp and the magnetic field modulation permanent magnet motor of the present invention
  • Fig. 6 (a) is the analysis diagram of the harmonic sensitivity of the armature magnetic field of the magnetic field modulation permanent magnet motor of the present invention
  • Fig. 6 (b) is the experimental point distribution figure of magnetic field modulation permanent magnet motor armature magnetic field harmonic of the present invention.
  • Fig. 7 (a) is the interaction effect diagram between the 1st armature magnetic field non-working harmonic and the 9th armature magnetic field non-working harmonic of the magnetic field modulation permanent magnet motor of the present invention
  • Fig. 7 (b) is the interaction effect diagram between the 1st non-working harmonic of the armature magnetic field and the 11 non-working harmonics of the armature magnetic field of the magnetic field modulation permanent magnet motor of the present invention
  • Fig. 7 (c) is the interaction effect diagram between the 1st non-working harmonic of the armature magnetic field and the 29 non-working harmonics of the armature magnetic field of the magnetic field modulation permanent magnet motor of the present invention
  • Fig. 7 (d) is the interaction effect diagram between the 1st armature magnetic field non-working harmonic and the 31st armature magnetic field working harmonic of the magnetic field modulation permanent magnet motor of the present invention
  • Fig. 8 is an analysis diagram of the design parameters of the magnetic field modulation permanent magnet motor of the present invention about the non-working harmonic sensitivity of the primary armature magnetic field with relative independence;
  • Fig. 9 is a non-working harmonic Kriging model of the primary armature magnetic field with relative independence of the magnetic field modulation permanent magnet motor of the present invention.
  • Fig. 10 is the distribution diagram of the Pareto frontier of the armature harmonics of the magnetic field modulation permanent magnet motor of the present invention.
  • Fig. 11 is a comparison chart of torque density and power factor before and after optimization of the magnetic field modulation permanent magnet motor of the present invention.
  • Fig. 1 is a flow chart of a permanent magnet-armature double harmonic collaborative optimization design method for a magnetic field modulation permanent magnet motor in an embodiment of the present invention. Referring to FIG. 1 , a design method for permanent magnet-armature double harmonic collaborative optimization of a magnetic field modulated permanent magnet motor in this embodiment will be described in detail.
  • a permanent magnet-armature double harmonic collaborative optimization design method for a magnetic field modulation permanent magnet motor according to the present invention, the specific implementation method is shown in Figure 1, including the following steps:
  • Step 1 According to the harmonic characteristics of the permanent magnet and armature magnetic fields, respectively establish the expressions between the flux linkage and the harmonics of the permanent magnet and armature magnetic fields; according to the magnetic field modulation permanent magnet motor flux vector diagram ( Figure 2), establish The expression of the power factor on the permanent magnet and armature flux linkage, according to the change of the size of the torque working area with the permanent magnet and armature flux linkage in the vector diagram, combined with the expression to analyze the corresponding magnetic field of the permanent magnet and armature magnetic field harmonics Based on the influence of the chain on the torque and power factor, the permanent magnet-armature double harmonic collaborative optimization design idea is established to improve the torque density and power factor of the motor.
  • the permanent magnetic field harmonics are divided into two types: permanent magnetic field working harmonics and permanent magnetic field non-working harmonics.
  • the harmonics of the permanent magnetic field that contribute to the torque are the working harmonics of the permanent magnetic field, and the harmonics of the permanent magnetic field that do not contribute to the torque are the non-working harmonics of the permanent magnetic field.
  • the working harmonics of the synthetic permanent magnetic field can be obtained by weighting the working harmonics of the permanent magnetic field.
  • the working harmonic amplitude of the magnetic field can maintain a high torque density while improving the power factor; for the armature magnetic field harmonics, the armature magnetic field
  • the armature magnetic field There are two types of armature magnetic field working harmonics and armature magnetic field non-working harmonics.
  • the armature magnetic field harmonics that contribute to the torque are the armature magnetic field working harmonics, and the armature magnetic field harmonics that do not contribute to the torque are the armature magnetic field non-working harmonics.
  • the flux linkage corresponding to the non-working harmonic of the armature magnetic field is inversely proportional to the power factor, and according to the flux linkage phasor diagram, reducing the non-working harmonic of the armature magnetic field will not affect the torque corresponding to The size of the working area, so reducing the non-working harmonics of the armature field can improve the power factor without losing torque density.
  • the minimum value of the non-working harmonic of the armature magnetic field and the working harmonic of the armature magnetic field are determined.
  • the optimization of the maximum value can realize the collaborative optimization design of the double harmonics of the permanent magnet-armature magnetic field, thereby improving the torque density and power factor of the motor.
  • the expressions of the permanent magnet magnetic field harmonic characteristic formula B m ( ⁇ , t) and the armature magnetic field harmonic characteristic formula B a ( ⁇ , t) in step 1 are:
  • C m is the Fourier coefficient of the permanent magnetomotive force
  • D i and D j are the Fourier coefficients of the armature magnetomotive force
  • m is the order of the permanent magnetomotive force
  • k is the order of the permeability
  • i and j are the order of the magnetomotive force of the armature
  • P r is the number of pole pairs of the permanent magnet
  • ⁇ r is the mechanical speed of the motor
  • t is the time
  • ⁇ 0 and ⁇ k are the Fourier coefficients of the air gap permeance
  • N s is the electrical Number of pivot slots.
  • the harmonic order of the permanent magnet magnetic field is mP r , mP r ⁇ kN s
  • the harmonic order of the armature magnetic field is i, j, i ⁇ kN s , j ⁇ kN s
  • F m ( ⁇ ,t) is the permanent magnet magnetomotive force expression
  • ⁇ s ( ⁇ ) is the air gap permeance expression, which can be expressed as:
  • r g is the air gap radius
  • l ef is the axial length
  • n c is the number of turns of the winding.
  • C 1 is the Fourier coefficient of the fundamental wave of the permanent magnetomotive force.
  • step 1 the expression of the armature flux linkage with respect to the harmonics of the armature magnetic field in step 1 is:
  • step 1 the expression of the power factor with respect to the permanent magnet and armature flux linkage in step 1 is:
  • the corresponding flux linkage for the working harmonic of the armature magnetic field is the flux linkage corresponding to the non-working harmonic of the armature magnetic field
  • U and ⁇ r are the phase voltage and frequency respectively
  • E 0 is the permanent magnet back EMF.
  • Step 2 Limit the working harmonic amplitude of the synthetic permanent magnetic field from the perspective of the permanent magnetic field. Taking the minimum value of the working harmonic amplitude of the synthetic permanent magnetic field as a constraint condition, the design parameters that have a greater impact on the working harmonics of the permanent magnetic field are selected through sensitivity analysis.
  • a magnetic field modulation permanent magnet motor is selected as the implementation object of this optimization design (Fig. 3), and the sensitivity calculation formula in step 2 is as follows:
  • Y(x) represents the harmonic amplitude of the permanent magnet and armature magnetic field under different design parameters
  • N is the sampling number
  • x is the motor design parameters, including h m , w m , w p , w s , w a , h b , and w b .
  • the specific steps of the sensitivity analysis in step 2 are as follows: First, use the central composite design sampling method to sample the second-order factor design points, axis points and zero-level center points that meet the second-order regression rotation criterion, providing a basis for the sensitivity calculation. data support. Then, the sensitivity of different design parameters to the working harmonics of the permanent magnetic field is calculated through the sensitivity formula, and the design parameters with greater sensitivity are selected as the design parameters that limit the amplitude of the synthetic permanent magnetic field working harmonics.
  • the optimized model of the synthetic permanent magnetic field working harmonics established by selecting the design parameters with high sensitivity is expressed as:
  • H sm is the harmonic amplitude of the synthetic permanent magnet magnetic field of the motor
  • x2 represents the design parameter with high sensitivity to the harmonic of the permanent magnet magnetic field of the motor
  • H mh is the working harmonic amplitude of the hth permanent magnet magnetic field value
  • a h is the weighting coefficient corresponding to the working harmonic of the permanent magnetic field
  • g 1 is the minimum constraint value of the synthetic permanent magnetic field working harmonic.
  • Fig. 5 shows the Kriging model of the working harmonics of the synthesized permanent magnet magnetic field of the magnetic field modulated permanent magnet motor according to the embodiment of the present invention.
  • Step 3 Simplify the armature field harmonics optimization objective. Since there are many harmonic orders of the armature magnetic field to be optimized, it is necessary to simplify the optimization objectives of the armature magnetic field harmonics and reduce the number of optimization objectives.
  • First analyze the sensitivity of the armature magnetic field harmonics to the performance of the motor, select the more sensitive armature magnetic field harmonics as the optimization target, and further calculate the distribution of the armature magnetic field harmonics according to the experimental design method. According to the distribution diagram of wave experiment points, the non-operating harmonics of the armature magnetic field, which are different from the changing trend of the working harmonics of the armature magnetic field, are selected as the harmonics to be optimized, and the simplified optimization target of the armature magnetic field harmonics is obtained by weighting.
  • the simplified optimization target of armature magnetic field harmonics in step 3 is obtained from the sensitivity analysis and the distribution of experimental point distribution of armature magnetic field harmonics.
  • the central composite design sampling method is used for sampling to provide data support for the sensitivity calculation.
  • the sensitivity of different armature magnetic field harmonics to power factor is calculated by the sensitivity formula (Fig. 6(a)), and the magnetic field harmonic with greater sensitivity is selected as the target of armature harmonic optimization.
  • draw the distribution diagram of the experimental points of the armature magnetic field harmonics (Fig. 6(b)), as shown in Fig.
  • the 29th non-working harmonic of the armature magnetic The working harmonics of the armature magnetic field have the same change trend, so the 29th non-working harmonics of the armature magnetic field are not used as the optimization target, and the rest of the non-working harmonics of the armature magnetic field are selected as the harmonics that need to be optimized, and then these harmonics are optimized Linear weighting is used as the optimization objective to reduce the number of harmonic optimization objectives of the armature magnetic field, and realize the simplification of the harmonic optimization objectives of the armature magnetic field.
  • Step 4 Independent judgment of non-working harmonics of armature magnetic field.
  • By calculating the interaction effect between the armature magnetic field harmonics it is analyzed and judged whether there are relatively independent non-working harmonics in the armature magnetic field harmonics. If there are relatively independent non-working harmonics of the armature magnetic field, go to steps 5.1 and 5.2; if there are no relatively independent non-working harmonics of the armature magnetic field, go to step 5.3.
  • Figure 7 analyzes the interaction effect between the non-operating wave of the primary armature magnetic field and other harmonics of the armature magnetic field of the magnetic field modulated permanent magnet motor. As shown in the figure, the 1st non-working harmonic of the armature magnetic field is independent of other harmonics, so the 1st non-working wave of the armature magnetic field is a relatively independent non-working harmonic of the armature magnetic field.
  • Step 5.1 Because there is a relatively independent armature magnetic field non-working harmonic, optimize it separately from the rest of the armature magnetic field harmonics to reduce the design parameters and optimize the target dimension, and improve the armature magnetic field harmonic optimization the accuracy of the results. Firstly, optimize the first-order non-working harmonic of the armature magnetic field with relative independence, and use the sensitivity analysis to select the design parameters that are more sensitive to the first-order non-operating harmonic of the armature magnetic field with relative independence, as shown in Fig.
  • the design parameters h b and w b have greater sensitivity, so based on h b and w b , the Kriging model of the design parameters with respect to the 1st order non-working harmonic of the armature magnetic field with relative independence is established. As shown in Figure 9, according to the established Kriging model, the optimal design point of non-working harmonics of the armature magnetic field with relative independence is selected.
  • step 5.1 the optimized model expression of the non-working harmonics of the armature magnetic field with relative independence in step 5.1 is:
  • x 3 is the design parameters h b and w b that have a great influence on the non-operating harmonics of the relatively independent armature magnetic field after the sensitivity analysis
  • H a1 is the first relatively independent armature magnetic field
  • Non-working harmonics ⁇ 1 is the weighting coefficient corresponding to the relatively independent non-working harmonics of the armature magnetic field
  • the minimum value of the optimization model is the optimization target.
  • Step 5.2 Optimize the working harmonics of the armature magnetic field and other non-working harmonics of the armature magnetic field that do not have relative independence. Taking the design parameter range limited by the working harmonic amplitude of the synthetic permanent magnet magnetic field as the constraint condition, and taking the working harmonic of the armature magnetic field and the simplified non-working harmonic of the armature magnetic field as the optimization objectives, a multi-objective genetic algorithm is adopted The working harmonics and non-working harmonics of the motor armature magnetic field are optimized, and the optimal Pareto frontier diagram of the armature harmonics is obtained ( Figure 10), and finally it is determined that there is The motor design scheme with optimal armature harmonics realizes the collaborative optimization design of permanent magnet-armature magnetic field double harmonics, thereby improving the torque density and power factor of the motor. As shown in Figure 11, the torque density of the optimized magnetic field modulation permanent magnet motor increases from 26.2Nm/L to 30.9Nm/L, and the power factor increases from 0.49 to 0.7, which verifies the effectiveness of the proposed optimal design method
  • H aPr is the working harmonic of the armature magnetic field of the motor P r , and its maximum value is set as the optimization target, so as to ensure that the optimized motor has a higher torque.
  • ⁇ 9 H a9 (x 4 )+ ⁇ 11 H a11 (x 4 ) is the simplified non-working harmonic of the synthesized armature magnetic field of the motor, and its minimum value is set as the optimization target.
  • x 4 represents the remaining design parameters w s , w p , w m , w a , h m after excluding the design parameter x 3 from the total design parameters x 1
  • H a9 and H a11 are the 9th and 11th times
  • the non-working harmonics of the armature magnetic field, ⁇ 9 and ⁇ 11 are the weighting coefficients corresponding to the 9th and 11th non-working harmonics of the armature magnetic field.
  • H a31 is the 31st working harmonic of the armature magnetic field
  • ⁇ 31 is the weighting coefficient corresponding to the 31st working harmonic of the armature magnetic field.
  • f 1 (x 2 ) is a function of the value range of the design parameters after the working harmonics of the synthetic permanent magnet magnetic field are limited.
  • references to the terms “one embodiment,” “some embodiments,” “exemplary embodiments,” “example,” “specific examples,” or “some examples” are intended to mean that the implementation A specific feature, structure, material, or characteristic described by an embodiment or example is included in at least one embodiment or example of the present invention.
  • schematic representations of the above terms do not necessarily refer to the same embodiment or example.
  • the specific features, structures, materials or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Permanent Field Magnets Of Synchronous Machinery (AREA)

Abstract

Disclosed is a method for collaborative optimization design of permanent magnet-armature double harmonics of a magnetic field-modulated permanent magnet motor. An expression between a flux linkage and a magnetic field harmonic is established according to a harmonic characteristic formula of permanent magnet and armature magnetic fields; and the influence of the flux linkages corresponding to the permanent magnet and armature magnetic field harmonics on torque and power factor is analyzed according to a vector diagram and a power factor expression, so that a collaborative optimization design thought for permanent magnet-armature double harmonics is established. The dimensions of armature magnetic field harmonic optimization objectives and design parameters are reduced by means of sensitivity analysis, experimental point distribution calculation, and independence judgment for an armature magnetic field non-working harmonic; a parameter range constrained by the permanent magnet magnetic field harmonic is taken as a constraint condition to optimize the armature magnetic field harmonic by using a Kriging model and a multi-objective optimization algorithm, so that collaborative optimization design for the permanent magnet-armature magnetic field double harmonics is achieved, thereby improving the torque density and power factor of the motor.

Description

一种磁场调制永磁电机永磁-电枢双谐波协同优化设计方法A permanent magnet-armature dual-harmonic collaborative optimization design method for magnetic field modulated permanent magnet motors 所属技术领域Technical field
本发明涉及到磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,属于电机设计领域,具体适用于电动汽车、风力发电、船舶推进等要求高转矩密度、高功率因数的电机系统。The invention relates to a permanent magnet-armature dual-harmonic collaborative optimization design method for a magnetic field modulation permanent magnet motor, which belongs to the field of motor design and is specifically applicable to motors requiring high torque density and high power factor such as electric vehicles, wind power generation, ship propulsion, etc. system.
背景技术Background technique
近年来,随着电动汽车、风力发电、船舶推进等领域的迅速发展,目前市场对于驱动效率较高的直驱式电机需求越来越大。磁场调制永磁电机由于“磁场调制效应”拥有高转矩密度的特点。所谓的“磁场调制效应”是指磁场调制永磁电机通过其调制极的作用,调制产生多种能够产生转矩的工作磁场谐波,从而提高电机的转矩。磁场调制永磁电机以其高转矩密度的优点在直驱式电机领域拥有着巨大的发展潜力。但是,磁场调制永磁电机本身同时存在漏磁较高的问题,使得功率因数较低,制约了其实际应用。In recent years, with the rapid development of electric vehicles, wind power generation, ship propulsion and other fields, the current market demand for direct drive motors with higher driving efficiency is increasing. Field modulated permanent magnet motors are characterized by high torque density due to the "field modulation effect". The so-called "magnetic field modulation effect" means that the magnetic field modulation permanent magnet motor modulates and generates a variety of working magnetic field harmonics that can generate torque through the action of its modulation pole, thereby increasing the torque of the motor. Field modulated permanent magnet motors have great development potential in the field of direct drive motors due to their high torque density. However, the magnetic field modulation permanent magnet motor itself has the problem of high flux leakage, which makes the power factor low and restricts its practical application.
目前提高磁场调制永磁电机功率因数的设计方法主要是从提高永磁磁场的角度提高永磁磁链,或者降低电枢磁场的角度降低电枢电流。目前的设计方法仅从永磁磁场或电枢磁场单个角度提升磁场调制永磁电机功率因数,从最终结果来看,在某些方法中功率因数的提升会导致转矩密度的下降,而有些方法虽提升了电机功率因数提升,但提升的程度往往未达到期望,仍有较大的提升空间。这主要是因为磁场调制永磁电机的永磁磁场和电枢磁场均对对电机的功率因数具有重要影响,忽略任何一个磁场对功率因数的影响,都会使得磁场调制永磁电机的功率因数无法达到最大。因此,针对磁场调制永磁电机,亟待从永磁-电枢磁场两方面同时考虑提出行之有效的功率因数提升方法。At present, the design method to improve the power factor of the magnetic field modulation permanent magnet motor is mainly to increase the permanent magnet flux linkage from the perspective of increasing the permanent magnet magnetic field, or to reduce the armature current from the perspective of reducing the armature magnetic field. The current design method only improves the power factor of the magnetic field modulated permanent magnet motor from a single angle of the permanent magnet magnetic field or the armature magnetic field. From the final result, in some methods, the increase of the power factor will lead to the decrease of the torque density, while some methods Although the power factor of the motor has been improved, the degree of improvement often does not meet expectations, and there is still a large room for improvement. This is mainly because both the permanent magnet magnetic field and the armature magnetic field of the magnetic field modulation permanent magnet motor have an important influence on the power factor of the motor. Ignoring the influence of any magnetic field on the power factor will make the power factor of the magnetic field modulation permanent magnet motor unable to achieve maximum. Therefore, for the magnetic field modulation permanent magnet motor, it is urgent to propose an effective power factor improvement method considering both the permanent magnet and the armature magnetic field.
发明内容Contents of the invention
本发明的内容是根据磁场调制永磁电机永磁和电枢磁场的谐波特征,建立磁链与磁场谐波之间的表达式;根据磁场调制永磁电机磁链向量图,建立功率因数关于永磁和电枢磁链的表达式,获得功率因数关于永磁和电枢磁场谐波的表达式。根据向量图和功率因数表达式,分析永磁和电枢磁场谐波所对应磁链对转矩和功率因数的影响,建立出永磁-电枢双谐波协同优化设计思路。建立合成永磁磁场工作谐波幅值关于设计参数的优化模型,确定永磁磁场谐波限制下对应的设计参数的范围;通过敏感度分析、实验点分布计算和电枢磁 场非工作谐波的独立性判断,减小电枢磁场谐波优化目标和设计参数的维度;将经过永磁磁场谐波约束后的参数范围作为约束条件,利用克里格模型和多目标优化算法对电枢磁场谐波进行优化,实现对永磁-电枢磁场双谐波的协同优化设计,从而提升电机的转矩密度和功率因数。The content of the present invention is to set up the expression between the flux linkage and the magnetic field harmonic according to the harmonic characteristics of the permanent magnet of the magnetic field modulation permanent magnet motor and the armature magnetic field; Expressions for permanent magnet and armature flux linkages to obtain expressions for power factor with respect to permanent magnet and armature field harmonics. According to the vector diagram and the power factor expression, the influence of the flux linkage corresponding to the permanent magnet and armature magnetic field harmonics on the torque and power factor is analyzed, and the permanent magnet-armature double harmonic collaborative optimization design idea is established. Establish the optimization model of the working harmonic amplitude of the synthetic permanent magnetic field with respect to the design parameters, and determine the range of the corresponding design parameters under the limit of the permanent magnetic field harmonics; through sensitivity analysis, calculation of experimental point distribution and non-working harmonics of the armature Independent judgment, reducing the dimension of the armature magnetic field harmonic optimization target and design parameters; using the parameter range after the permanent magnet magnetic field harmonic constraints as constraints, using the kriging model and multi-objective optimization algorithm to The wave is optimized to realize the collaborative optimization design of the permanent magnet-armature magnetic field double harmonic, thereby improving the torque density and power factor of the motor.
本发明的技术方案包括以下步骤:Technical scheme of the present invention comprises the following steps:
步骤1:根据永磁磁场谐波特征公式和电枢磁场谐波特征公式,建立永磁磁链关于永磁磁场谐波之间的表达式,以及建立电枢磁链关于电枢磁场谐波之间的表达式;根据磁场调制永磁电机磁链向量图,建立功率因数关于永磁和电枢磁链的表达式,获得功率因数关于永磁和电枢磁场谐波的表达式。根据向量图中转矩工作区域大小随永磁和电枢磁链变化情况,同时结合以上表达式分析永磁和电枢磁场谐波所对应磁链对转矩和功率因数的影响,建立出永磁-电枢双谐波协同优化设计思路,提升电机的转矩密度和功率因数;Step 1: According to the permanent magnet magnetic field harmonic characteristic formula and the armature magnetic field harmonic characteristic formula, establish the expression between the permanent magnet flux linkage with respect to the permanent magnet magnetic field harmonics, and establish the relationship between the armature flux linkage with respect to the armature magnetic field harmonics According to the flux vector diagram of the magnetic field modulation permanent magnet motor, the expression of the power factor with respect to the permanent magnet and armature flux linkage is established, and the expression of the power factor with respect to the harmonics of the permanent magnet and armature magnetic field is obtained. According to the change of the size of the torque working area with the permanent magnet and armature flux linkage in the vector diagram, and combining the above expressions to analyze the influence of the flux linkage corresponding to the harmonics of the permanent magnet and armature magnetic field on the torque and power factor, the permanent magnet and armature flux linkage is established. Magneto-armature double harmonic synergistic optimization design idea to improve the torque density and power factor of the motor;
对于永磁磁场谐波而言,从永磁磁场谐波对于转矩是否有贡献的角度将永磁磁场谐波分为永磁磁场工作谐波和永磁磁场非工作谐波两类。对转矩有贡献的永磁磁场谐波为永磁磁场工作谐波,对转矩没有贡献的永磁磁场谐波为永磁磁场非工作谐波。对各次永磁磁场工作谐波进行加权可得到合成永磁磁场工作谐波,合成永磁磁场工作谐波所对应的永磁磁链与转矩和功率因数都成正比,因此提高合成永磁磁场工作谐波幅值可在提高功率因数的同时维持较高转矩密度;对于电枢磁场谐波而言,从电枢磁场谐波对于转矩是否有贡献的角度将电枢磁场谐波分为电枢磁场工作谐波和电枢磁场非工作谐波两类。对转矩有贡献的电枢磁场谐波为电枢磁场工作谐波,对转矩没有贡献的电枢磁场谐波为电枢磁场非工作谐波。根据磁链和功率因数表达式,电枢磁场非工作谐波所对应的磁链与功率因数成反比,同时根据磁链相量图,降低电枢磁场非工作谐波不会影响转矩所对应的工作区域大小,因此减少电枢磁场非工作谐波可以提升功率因数,同时不损失转矩密度。For the permanent magnetic field harmonics, from the perspective of whether the permanent magnetic field harmonics contribute to the torque, the permanent magnetic field harmonics are divided into two types: permanent magnetic field working harmonics and permanent magnetic field non-working harmonics. The harmonics of the permanent magnetic field that contribute to the torque are the working harmonics of the permanent magnetic field, and the harmonics of the permanent magnetic field that do not contribute to the torque are the non-working harmonics of the permanent magnetic field. The working harmonics of the synthetic permanent magnetic field can be obtained by weighting the working harmonics of the permanent magnetic field. The working harmonic amplitude of the magnetic field can maintain a high torque density while improving the power factor; for the armature magnetic field harmonics, the armature magnetic field There are two types of armature magnetic field working harmonics and armature magnetic field non-working harmonics. The armature magnetic field harmonics that contribute to the torque are the armature magnetic field working harmonics, and the armature magnetic field harmonics that do not contribute to the torque are the armature magnetic field non-working harmonics. According to the flux linkage and power factor expressions, the flux linkage corresponding to the non-working harmonic of the armature magnetic field is inversely proportional to the power factor, and according to the flux linkage phasor diagram, reducing the non-working harmonic of the armature magnetic field will not affect the torque corresponding to The size of the working area, so reducing the non-working harmonics of the armature field can improve the power factor without losing torque density.
根据针对电机向量图和功率因数表达式的分析,确定出在保证合成永磁磁场工作谐波处于较高水平的前提下,通过对电枢磁场非工作谐波最小值和电枢磁场工作谐波最大值的优化,即可实现对永磁-电枢磁场双谐波的协同优化设计,从而提升电机的转矩密度和功率因数。According to the analysis of the vector diagram of the motor and the power factor expression, it is determined that under the premise of ensuring that the working harmonic of the synthetic permanent magnet magnetic field is at a high level, the minimum value of the non-working harmonic of the armature magnetic field and the working harmonic of the armature magnetic field are determined. The optimization of the maximum value can realize the collaborative optimization design of the double harmonics of the permanent magnet-armature magnetic field, thereby improving the torque density and power factor of the motor.
步骤2:从永磁磁场的角度对合成永磁磁场工作谐波幅值进行限制。将合成永磁磁场工作谐波幅值的最小值作为约束条件,通过敏感度分析选择出对永磁磁场工作谐波影响较大的设计参数,基于对永磁磁场工作谐波具有高敏感度的设计参数,建立合成永磁磁场工作谐波的优化模型,利用克里格模型表示合成永磁磁场工作谐波幅值与对永磁磁场工作谐波 具有高敏感度的设计参数之间的关系,根据优化模型中设定的最小合成永磁磁场工作谐波幅值,基于所建立的克里格模型得出对应的对永磁磁场工作谐波具有高敏感度的设计参数的变化范围。Step 2: Limit the working harmonic amplitude of the synthetic permanent magnetic field from the perspective of the permanent magnetic field. Taking the minimum value of the working harmonic amplitude of the synthetic permanent magnetic field as a constraint condition, the design parameters that have a greater impact on the working harmonics of the permanent magnetic field are selected through sensitivity analysis. Based on the high sensitivity to the working harmonics of the permanent magnetic field Design parameters, establish an optimization model for the working harmonics of the synthetic permanent magnetic field, and use the Kriging model to express the relationship between the working harmonic amplitude of the synthetic permanent magnetic field and the design parameters with high sensitivity to the working harmonics of the permanent magnetic field, According to the minimum synthetic permanent magnet magnetic field working harmonic amplitude set in the optimization model, based on the established Kriging model, the corresponding variation range of the design parameters with high sensitivity to the permanent magnetic magnetic field working harmonic is obtained.
步骤3:简化电枢磁场谐波优化目标。由于进行优化的电枢磁场谐波阶次较多,需要对电枢磁场谐波优化目标进行简化,降低优化目标个数。首先,分析电枢磁场谐波关于电机性能的敏感度,选择敏感度较大的电枢磁场谐波作为优化目标,进一步根据实验设计法计算出电枢磁场谐波的分布,通过电枢磁场谐波实验点分布图,选择出电枢磁场工作谐波变化趋势不同的电枢磁场非工作谐波为需要优化的谐波,并采用加权的方式得到简化的电枢磁场谐波优化目标。Step 3: Simplify the armature field harmonics optimization objective. Since there are many harmonic orders of the armature magnetic field to be optimized, it is necessary to simplify the optimization objectives of the armature magnetic field harmonics and reduce the number of optimization objectives. First, analyze the sensitivity of the armature magnetic field harmonics to the performance of the motor, select the more sensitive armature magnetic field harmonics as the optimization target, and further calculate the distribution of the armature magnetic field harmonics according to the experimental design method. According to the distribution diagram of wave experiment points, the non-operating harmonics of the armature magnetic field with different changing trends of the working harmonics of the armature magnetic field are selected as the harmonics to be optimized, and the simplified optimization target of the armature magnetic field harmonics is obtained by weighting.
步骤4:电枢磁场非工作谐波的独立性判断。通过计算电枢磁场谐波之间的交互效应,分析判断电枢磁场谐波中是否存在具有相对独立性的非工作谐波。如果存在相对独立性的电枢磁场非工作谐波,进入步骤5.1和5.2;如果不存在相对独立性的电枢磁场非工作谐波,进入步骤5.3Step 4: Independent judgment of non-working harmonics of armature magnetic field. By calculating the interaction effect between the armature magnetic field harmonics, it is analyzed and judged whether there are relatively independent non-working harmonics in the armature magnetic field harmonics. If there are relatively independent non-working harmonics of the armature magnetic field, go to steps 5.1 and 5.2; if there are no relatively independent non-working harmonics of the armature magnetic field, go to step 5.3
步骤5.1:如果存在具有相对独立性的电枢磁场非工作谐波,将具有相对独立性的电枢磁场非工作谐波与其余电枢磁场谐波分开进行优化,以减小设计参数和优化目标维度,提高电枢磁场谐波优化结果的准确性。首先,对具有相对独立性的电枢磁场非工作谐波进行优化,利用敏感度分析选择对具有相对独立性的电枢磁场非工作谐波敏感度较大的设计参数,建立设计参数关于具有相对独立性的电枢磁场非工作谐波的克里格模型,根据所建立的克里格模型选择出具有相对独立性的电枢磁场非工作谐波的最优设计点。Step 5.1: If there are non-operating harmonics of the armature field with relative independence, optimize the non-operating harmonics of the armature field with relative independence separately from the rest of the harmonics of the armature field to reduce the design parameters and optimization objectives dimension to improve the accuracy of the armature field harmonic optimization results. First, optimize the non-working harmonics of the armature magnetic field with relative independence, use sensitivity analysis to select the design parameters that are more sensitive to the non-working harmonics of the armature magnetic field with relative independence, and establish the design parameters with relative Independent Kriging model of non-working harmonics of armature magnetic field. According to the established Kriging model, the optimal design point of non-working harmonics of armature magnetic field with relative independence is selected.
步骤5.2:对于电枢磁场工作谐波与其余不具备相对独立性的电枢磁场的非工作谐波进行优化。将经过合成永磁磁场工作谐波幅值限制后的设计参数范围作为约束条件,并将电枢磁场工作谐波与经过简化后的电枢磁场非工作谐波作为优化目标,采用多目标遗传算法对电机电枢磁场工作谐波和非工作谐波进行优化,最终确定出在高合成永磁工作谐波幅值基础上具有最优电枢谐波的电机设计方案,实现对永磁-电枢磁场双谐波的协同优化设计,从而提升电机的转矩密度和功率因数。Step 5.2: Optimize the working harmonics of the armature magnetic field and other non-working harmonics of the armature magnetic field that do not have relative independence. Taking the design parameter range limited by the working harmonic amplitude of the synthetic permanent magnet magnetic field as the constraint condition, and taking the working harmonic of the armature magnetic field and the simplified non-working harmonic of the armature magnetic field as the optimization objectives, a multi-objective genetic algorithm is adopted Optimize the working harmonics and non-working harmonics of the motor armature magnetic field, and finally determine the motor design scheme with the optimal armature harmonics on the basis of the high synthetic permanent magnet working harmonic amplitude, and realize the permanent magnet-armature The collaborative optimization design of the double harmonics of the magnetic field improves the torque density and power factor of the motor.
步骤5.3:对于电枢磁场的工作谐波及非工作谐波进行优化。将经过合成永磁磁场工作谐波幅值限制后的参数范围作为约束条件,将电枢磁场工作谐波与经过简化后的电枢磁场非工作谐波作为优化目标,采用多目标遗传算法对电机电枢磁场工作谐波和非工作谐波进行优化,最终确定出在高合成永磁工作谐波幅值基础上具有最优电枢谐波的电机设计方案,实现对永磁-电枢磁场双谐波的协同优化设计,从而提升电机的转矩密度和功率因数。Step 5.3: Optimize the working harmonics and non-working harmonics of the armature magnetic field. Taking the parameter range limited by the working harmonic amplitude of the synthesized permanent magnet magnetic field as the constraint condition, taking the working harmonic of the armature magnetic field and the simplified non-working harmonic of the armature magnetic field as the optimization target, the multi-objective genetic algorithm is used to optimize the The working harmonics and non-working harmonics of the armature magnetic field of the machine are optimized, and finally the design scheme of the motor with the optimal armature harmonics is determined on the basis of the high synthetic permanent magnet working harmonic amplitude. The coordinated optimization design of harmonics can improve the torque density and power factor of the motor.
进一步,步骤1中永磁磁场谐波特征公式B m(θ,t)和电枢磁场谐波特征公式B a(θ,t)的表达式为: Further, the expressions of the permanent magnet magnetic field harmonic characteristic formula B m (θ, t) and the armature magnetic field harmonic characteristic formula B a (θ, t) in step 1 are:
Figure PCTCN2022070687-appb-000001
Figure PCTCN2022070687-appb-000001
Figure PCTCN2022070687-appb-000002
Figure PCTCN2022070687-appb-000002
式中,C m是永磁磁动势傅里叶系数,D i和D j是电枢磁动势傅里叶系数,m是永磁磁动势阶次,k是磁导阶次,i和j为电枢磁动势阶次,P r是永磁体极对数,Ω r是电机机械转速,t为时间,Λ 0和Λ k是气隙磁导傅里叶系数,N s是电枢槽数。根据表达式可确定出永磁磁场谐波阶次为mP r,mP r±kN s,电枢磁场谐波阶次为i,j,i±kN s,j±kN s。F m(θ,t)为永磁磁动势表达式,Λ s(θ)为气隙磁导表达式,分别可以表示为: In the formula, C m is the Fourier coefficient of the permanent magnetomotive force, D i and D j are the Fourier coefficients of the armature magnetomotive force, m is the order of the permanent magnetomotive force, k is the order of the permeability, i and j are the order of the magnetomotive force of the armature, P r is the number of pole pairs of the permanent magnet, Ω r is the mechanical speed of the motor, t is the time, Λ 0 and Λ k are the Fourier coefficients of the air gap permeance, N s is the electrical Number of pivot slots. According to the expression, it can be determined that the harmonic order of the permanent magnet magnetic field is mP r , mP r ±kN s , and the harmonic order of the armature magnetic field is i, j, i±kN s , j±kN s . F m (θ,t) is the permanent magnet magnetomotive force expression, Λ s (θ) is the air gap permeance expression, which can be expressed as:
Figure PCTCN2022070687-appb-000003
Figure PCTCN2022070687-appb-000003
Figure PCTCN2022070687-appb-000004
Figure PCTCN2022070687-appb-000004
进一步,步骤1中永磁磁链关于永磁磁场谐波的表达式为:Further, the expression of the permanent magnet flux linkage on the harmonics of the permanent magnet field in step 1 is:
Figure PCTCN2022070687-appb-000005
Figure PCTCN2022070687-appb-000005
式中,r g为气隙半径,l ef为轴向长度,n c为绕组匝数。其中,基波永磁磁链幅值表达式为: In the formula, r g is the air gap radius, l ef is the axial length, n c is the number of turns of the winding. Among them, the amplitude expression of the fundamental permanent magnet flux linkage is:
Figure PCTCN2022070687-appb-000006
Figure PCTCN2022070687-appb-000006
进一步,步骤1中电枢磁链关于电枢磁场谐波的表达式为:Further, the expression of the armature flux linkage with respect to the harmonics of the armature magnetic field in step 1 is:
Figure PCTCN2022070687-appb-000007
Figure PCTCN2022070687-appb-000007
其中,基波电枢磁链幅值表达式为:Among them, the amplitude expression of the fundamental armature flux linkage is:
Figure PCTCN2022070687-appb-000008
Figure PCTCN2022070687-appb-000008
进一步,步骤1中功率因数关于永磁和电枢磁链的表达式为:Further, the expression of the power factor with respect to the permanent magnet and armature flux linkage in step 1 is:
Figure PCTCN2022070687-appb-000009
Figure PCTCN2022070687-appb-000009
式中,
Figure PCTCN2022070687-appb-000010
为漏磁,
Figure PCTCN2022070687-appb-000011
为永磁磁链,
Figure PCTCN2022070687-appb-000012
为电枢磁场工作谐波对应磁链,
Figure PCTCN2022070687-appb-000013
为电枢磁场非工作谐波对应磁链,U和ω r分别为相电压和频率,E 0为永磁反电势。
In the formula,
Figure PCTCN2022070687-appb-000010
For flux leakage,
Figure PCTCN2022070687-appb-000011
is the permanent magnet flux linkage,
Figure PCTCN2022070687-appb-000012
The corresponding flux linkage for the working harmonic of the armature magnetic field,
Figure PCTCN2022070687-appb-000013
is the flux linkage corresponding to the non-working harmonic of the armature magnetic field, U and ω r are the phase voltage and frequency respectively, and E 0 is the back EMF of the permanent magnet.
进一步,步骤2中敏感度计算公式表达式为:Further, the expression of the sensitivity calculation formula in step 2 is:
Figure PCTCN2022070687-appb-000014
Figure PCTCN2022070687-appb-000014
式中,Y(x)表示的是不同设计参数下永磁和电枢磁场谐波幅值,N为采样数目,x为电机设计参数。In the formula, Y(x) represents the harmonic amplitude of the permanent magnet and armature magnetic field under different design parameters, N is the number of samples, and x is the motor design parameter.
步骤2中敏感度分析的具体步骤在于:首先,运用中心复合设计采样方法,对满足二阶回归旋转准则的二阶因子设计点,轴点以及零水平中心点进行采样,为进行敏感度计算提供数据支撑。然后,通过灵敏度公式计算不同设计参数对于永磁磁场工作谐波的敏感度, 选择敏感度较大的设计参数作为限制合成永磁磁场工作谐波幅值的设计参数。The specific steps of the sensitivity analysis in step 2 are as follows: First, use the central composite design sampling method to sample the second-order factor design points, axis points and zero-level center points that meet the second-order regression rotation criterion, providing a basis for the sensitivity calculation. data support. Then, the sensitivity of different design parameters to the working harmonics of the permanent magnetic field is calculated by the sensitivity formula, and the design parameters with greater sensitivity are selected as the design parameters that limit the amplitude of the synthetic permanent magnetic field working harmonics.
进行敏感度分析后建立的合成永磁磁场工作谐波的优化模型的表达式为:The expression of the optimization model of the working harmonic of the synthetic permanent magnetic field established after the sensitivity analysis is:
Figure PCTCN2022070687-appb-000015
Figure PCTCN2022070687-appb-000015
Constraint:H sm(x 2)>g 1 Constraint:H sm (x 2 )>g 1
式中,H sm为电机合成永磁磁场工作谐波幅值,x 2表示的是对电机永磁磁场谐波具有高敏感度的设计参数,H mh为第h次永磁磁场工作谐波幅值,a h为对应永磁磁场工作谐波的加权系数,g 1为合成永磁磁场工作谐波的最小约束值。 In the formula, H sm is the harmonic amplitude of the synthetic permanent magnet magnetic field of the motor, x2 represents the design parameter with high sensitivity to the harmonic of the permanent magnet magnetic field of the motor, and H mh is the working harmonic amplitude of the hth permanent magnet magnetic field value, a h is the weighting coefficient corresponding to the working harmonic of the permanent magnetic field, g 1 is the minimum constraint value of the synthetic permanent magnetic field working harmonic.
进一步,步骤3中简化后电枢磁场谐波优化目标由敏感度分析和电枢磁场谐波实验点分布得出。首先,运用中心复合设计采样方法进行采样,为敏感度计算提供数据支撑。然后,通过灵敏度公式计算不同电枢磁场谐波对于功率因数的敏感度,选择敏感度较大的磁场谐波作为电枢谐波优化目标。进一步根据复合设计采样方法中的实验点,绘制电枢磁场谐波实验点分布图,根据图中不同电枢磁场谐波的变化趋势,选择出与电枢磁场工作谐波变化趋势不同的电枢磁场非工作谐波为需要优化的谐波,而与电枢磁场工作谐波变化趋势相同的电枢磁场非工作谐波不作为需要优化的谐波。然后将这些需要优化的谐波进行线性加权作为优化目标,减少电枢磁场谐波优化目标的个数,实现对电枢磁场谐波优化目标的简化。Further, the simplified optimization target of armature magnetic field harmonics in step 3 is obtained from the sensitivity analysis and the distribution of experimental point distribution of armature magnetic field harmonics. First, the central composite design sampling method is used for sampling to provide data support for the sensitivity calculation. Then, the sensitivity of different armature magnetic field harmonics to power factor is calculated by the sensitivity formula, and the magnetic field harmonic with higher sensitivity is selected as the target of armature harmonic optimization. Further, according to the experimental points in the composite design sampling method, draw the distribution map of the experimental points of the armature magnetic field harmonics, and select the armature with a different trend from the working harmonics of the armature magnetic field according to the changing trend of the different armature magnetic field harmonics in the figure The non-working harmonics of the magnetic field are the harmonics that need to be optimized, and the non-working harmonics of the armature magnetic field that have the same change trend as the working harmonics of the armature magnetic field are not considered as the harmonics that need to be optimized. Then these harmonics that need to be optimized are linearly weighted as optimization targets, reducing the number of optimization targets for armature magnetic field harmonics, and realizing the simplification of armature magnetic field harmonic optimization targets.
进一步,步骤5.1中具有相对独立性的电枢磁场非工作谐波的优化模型表达式为:Further, the optimized model expression of the non-working harmonics of the armature magnetic field with relative independence in step 5.1 is:
Figure PCTCN2022070687-appb-000016
Figure PCTCN2022070687-appb-000016
式中,x 3为经过敏感度分析后对具有相对独立性的电枢磁场非工作谐波影响较大的设计参数,H al为第l次具有相对独立性的电枢磁场非工作谐波,λ l为对应具有相对独立性的电枢磁场非工作谐波的加权系数,优化模型的最小值为优化目标。 In the formula, x3 is a design parameter that has a great influence on the non-operating harmonic of the armature magnetic field with relative independence after sensitivity analysis, and Hal is the non-operating harmonic of the armature magnetic field with relative independence for the lth time, λ l is the weighting coefficient corresponding to the non-working harmonic of the armature magnetic field with relative independence, and the minimum value of the optimization model is the optimization target.
步骤5.2中电枢磁场工作谐波与简化后的其余不具备相对独立性的电枢磁场非工作谐波的优化模型表达式为:In step 5.2, the optimized model expression of the working harmonics of the armature magnetic field and the simplified non-working harmonics of the armature magnetic field that do not have relative independence is:
Objectives:
Figure PCTCN2022070687-appb-000017
Objectives:
Figure PCTCN2022070687-appb-000017
Constraint:f 1(x 2)>0 Constraint: f 1 (x 2 )>0
式中,H aPr为电机P r次电枢磁场工作谐波,其最大值被设定为优化目标,从而保证优化后电机具有较高转矩。
Figure PCTCN2022070687-appb-000018
为经过简化后的电机合成电枢磁场非工作谐波,其最小值被设定为优化目标。x 4表示的是在总设计参数x 1中排除掉设计参数x 3后所剩余的设计参数,H as为第s次电枢磁场非工作谐波,μ s为对应s次电枢磁场非工作谐波的加权系数。f 1(x 2)为经过合成永磁磁场工作谐波限制后的设计参数取值范围函数。
In the formula, H aPr is the working harmonic of the armature magnetic field of the motor P r , and its maximum value is set as the optimization target, so as to ensure that the optimized motor has a higher torque.
Figure PCTCN2022070687-appb-000018
The non-operating harmonics of the armature field are synthesized for the simplified motor, and its minimum value is set as the optimization target. x 4 represents the remaining design parameters after excluding the design parameter x 3 from the total design parameters x 1 , H as is the sth non-working harmonic of the armature magnetic field, μ s is the corresponding s-th non-working armature magnetic field Weighting coefficients for harmonics. f 1 (x 2 ) is a function of the value range of the design parameters after the working harmonics of the synthetic permanent magnet magnetic field are limited.
进一步,步骤5.3中电枢磁场工作谐波与经过简化后的电枢磁场非工作谐波的优化模型表达式为:Further, the optimized model expressions of the working harmonics of the armature magnetic field and the simplified non-working harmonics of the armature magnetic field in step 5.3 are:
Objectives:
Figure PCTCN2022070687-appb-000019
Objectives:
Figure PCTCN2022070687-appb-000019
Constraint:f 1(x 2)>0 Constraint: f 1 (x 2 )>0
式中,H aPr为电机P r次电枢磁场工作谐波,其最大值被设定为优化目标,从而保证优化后电机具有较高转矩。
Figure PCTCN2022070687-appb-000020
为经过简化后的电机合成电枢磁场非工作谐波,其最小值被设定为优化目标。x 1为电机总的设计参数,x 2表示的是对电机永磁磁场谐波具有高敏感度的设计参数,H as为第s次电枢磁场非工作谐波,μ s为对应s次电枢磁场非工作谐波的加权系数。f 1(x 2)为经过合成永磁磁场工作谐波限制后的设计参数取值范围函数。
In the formula, H aPr is the working harmonic of the armature magnetic field of the motor P r , and its maximum value is set as the optimization target, so as to ensure that the optimized motor has a higher torque.
Figure PCTCN2022070687-appb-000020
The non-operating harmonics of the armature field are synthesized for the simplified motor, and its minimum value is set as the optimization target. x1 is the general design parameter of the motor, x2 is the design parameter with high sensitivity to the harmonic of the permanent magnet field of the motor, H as is the sth non-working harmonic of the armature magnetic field, μ s is the corresponding sth electric Weighting factor for non-operating harmonics of the pivot field. f 1 (x 2 ) is a function of the value range of the design parameters after the working harmonics of the synthetic permanent magnet magnetic field are limited.
有益效果Beneficial effect
本发明采用上述设计方案后,可以具备如下有益效果:After the present invention adopts the above design scheme, it can have the following beneficial effects:
1)本发明根据永磁和电枢磁场的谐波特征,建立磁链与磁场谐波之间的表达式;根据磁场调制永磁电机磁链向量图,建立功率因数关于永磁和电枢磁链的表达式,获得功率因数关于永磁和电枢磁场谐波的表达式。根据向量图和功率因数表达式,分析永磁和电枢磁场谐波所对应磁链对转矩和功率因数的影响,建立出永磁-电枢双谐波协同优化设计思路,为后续从永磁和电枢双谐波的角度对转矩密度和功率因数提升指明了方向。1) The present invention sets up the expression between the flux linkage and the magnetic field harmonic according to the harmonic characteristics of the permanent magnet and the armature magnetic field; Expressions for the chain to obtain expressions for the power factor with respect to the harmonics of the permanent and armature fields. According to the vector diagram and the power factor expression, the influence of the flux linkage corresponding to the permanent magnet and armature magnetic field harmonics on the torque and power factor is analyzed, and the permanent magnet-armature double harmonic collaborative optimization design idea is established, which will be used for the follow-up from permanent magnet The angle of the magnetic and armature double harmonics gives directions for torque density and power factor improvement.
2)本发明利用实验点分布计算、灵敏度分析和电枢磁场非工作谐波的独立性判断,降低优化目标和设计参数的维度,进一步建立克里格模型,减少电机设计的计算量,最后,基于克里格模型和多目标优化算法对磁场调制永磁电机进行优化设计。与传统的优化方法相比,采用本发明设计方法可显著提高优化效率,减少优化时间。2) The present invention utilizes experimental point distribution calculation, sensitivity analysis and independent judgment of armature magnetic field non-working harmonics to reduce the dimension of optimization objectives and design parameters, further establish the kriging model, and reduce the calculation amount of motor design. Finally, Based on Kriging model and multi-objective optimization algorithm, the optimal design of magnetic field modulation permanent magnet motor is carried out. Compared with the traditional optimization method, the design method of the invention can significantly improve the optimization efficiency and reduce the optimization time.
3)本发明将合成永磁磁场工作谐波作为约束条件,得出满足约束条件的设计参数范围,在此基础上将电枢磁场谐波作为优化目标对其进行优化,实现永磁-电枢双谐波的角度协同优化设计,确定出在高合成永磁磁场工作谐波基础上的最优电枢磁场谐波的电机设计方案。与目前单独从永磁或电枢的角度优化功率因数的方法相比,采用本发明设计方法可进一步提高电机的转矩密度和功率因数。3) The present invention uses the synthetic permanent magnet magnetic field working harmonic as a constraint condition to obtain the design parameter range satisfying the constraint condition. On this basis, the armature magnetic field harmonic is used as an optimization target to optimize it to realize permanent magnet-armature The angle synergistic optimization design of double harmonics determines the motor design scheme for the optimal armature magnetic field harmonics based on the high synthetic permanent magnetic field working harmonics. Compared with the current method of optimizing power factor solely from the perspective of permanent magnet or armature, the design method of the invention can further improve the torque density and power factor of the motor.
附图说明Description of drawings
图1是本发明实施例中一种磁场调制永磁电机永磁-电枢双谐波协同优化设计方法流程图;Fig. 1 is a flow chart of a permanent magnet-armature double harmonic collaborative optimization design method for a magnetic field modulation permanent magnet motor in an embodiment of the present invention;
图2是本发明磁场调制永磁电机磁链向量图;Fig. 2 is the flux linkage vector diagram of the magnetic field modulation permanent magnet motor of the present invention;
图3是本发明磁场调制永磁电机拓扑结构及参数分布图;Fig. 3 is the topological structure and parameter distribution diagram of the magnetic field modulation permanent magnet motor of the present invention;
图4是本发明磁场调制永磁电机设计参数关于合成永磁磁场工作谐波的敏感度分析结果;Fig. 4 is the sensitivity analysis result of the design parameters of the magnetic field modulation permanent magnet motor of the present invention about the working harmonic of the synthetic permanent magnet magnetic field;
图5(a)是本发明磁场调制永磁电机在w m、h m、变化下,合成永磁磁场工作谐波的克里格模型计算结果; Fig. 5 (a) is the Kriging model calculation result of the working harmonic of the synthetic permanent magnet magnetic field under the change of w m , h m , of the magnetic field modulation permanent magnet motor of the present invention;
图5(b)是本发明磁场调制永磁电机在w s、w p、变化下,合成永磁磁场工作谐波的克里格模型计算结果; Fig. 5(b) is the Kriging model calculation result of the working harmonics of the synthetic permanent magnet magnetic field under the change of ws , wp and the magnetic field modulation permanent magnet motor of the present invention;
图6(a)是本发明磁场调制永磁电机电枢磁场谐波敏感度分析图;Fig. 6 (a) is the analysis diagram of the harmonic sensitivity of the armature magnetic field of the magnetic field modulation permanent magnet motor of the present invention;
图6(b)是本发明磁场调制永磁电机电枢磁场谐波实验点分布图;Fig. 6 (b) is the experimental point distribution figure of magnetic field modulation permanent magnet motor armature magnetic field harmonic of the present invention;
图7(a)是本发明磁场调制永磁电机1次电枢磁场非工作谐波与9次电枢磁场非工作谐波之间的交互效应图;Fig. 7 (a) is the interaction effect diagram between the 1st armature magnetic field non-working harmonic and the 9th armature magnetic field non-working harmonic of the magnetic field modulation permanent magnet motor of the present invention;
图7(b)是本发明磁场调制永磁电机1次电枢磁场非工作谐波与11次电枢磁场非工作谐波之间的交互效应图;Fig. 7 (b) is the interaction effect diagram between the 1st non-working harmonic of the armature magnetic field and the 11 non-working harmonics of the armature magnetic field of the magnetic field modulation permanent magnet motor of the present invention;
图7(c)是本发明磁场调制永磁电机1次电枢磁场非工作谐波与29次电枢磁场非工作谐波之间的交互效应图;Fig. 7 (c) is the interaction effect diagram between the 1st non-working harmonic of the armature magnetic field and the 29 non-working harmonics of the armature magnetic field of the magnetic field modulation permanent magnet motor of the present invention;
图7(d)是本发明磁场调制永磁电机1次电枢磁场非工作谐波与31次电枢磁场工作谐波之间的交互效应图;Fig. 7 (d) is the interaction effect diagram between the 1st armature magnetic field non-working harmonic and the 31st armature magnetic field working harmonic of the magnetic field modulation permanent magnet motor of the present invention;
图8是本发明磁场调制永磁电机设计参数关于具备相对独立性的1次电枢磁场非工作谐波敏感度分析图;Fig. 8 is an analysis diagram of the design parameters of the magnetic field modulation permanent magnet motor of the present invention about the non-working harmonic sensitivity of the primary armature magnetic field with relative independence;
图9是本发明磁场调制永磁电机具备相对独立性的1次电枢磁场非工作谐波克里格模型;Fig. 9 is a non-working harmonic Kriging model of the primary armature magnetic field with relative independence of the magnetic field modulation permanent magnet motor of the present invention;
图10是本发明磁场调制永磁电机电枢谐波的帕累托前沿分布图;Fig. 10 is the distribution diagram of the Pareto frontier of the armature harmonics of the magnetic field modulation permanent magnet motor of the present invention;
图11是本发明磁场调制永磁电机优化前后转矩密度和功率因数结果对比图。Fig. 11 is a comparison chart of torque density and power factor before and after optimization of the magnetic field modulation permanent magnet motor of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the drawings in the embodiments of the present invention.
下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。The embodiments described below by referring to the figures are exemplary only for explaining the present invention and should not be construed as limiting the present invention.
图1为本发明实施例中一种磁场调制永磁电机永磁-电枢双谐波协同优化设计方法流程图。参阅图1,对本实施例中一种磁场调制永磁电机永磁-电枢双谐波协同优化设计方法进行详细说明。Fig. 1 is a flow chart of a permanent magnet-armature double harmonic collaborative optimization design method for a magnetic field modulation permanent magnet motor in an embodiment of the present invention. Referring to FIG. 1 , a design method for permanent magnet-armature double harmonic collaborative optimization of a magnetic field modulated permanent magnet motor in this embodiment will be described in detail.
本发明所述一种磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,具体实施方法如图1,包括以下步骤:A permanent magnet-armature double harmonic collaborative optimization design method for a magnetic field modulation permanent magnet motor according to the present invention, the specific implementation method is shown in Figure 1, including the following steps:
步骤1:根据永磁和电枢磁场的谐波特征,分别建立磁链与永磁和电枢磁场谐波之间的表达式;根据磁场调制永磁电机磁链向量图(图2),建立功率因数关于永磁和电枢磁链的表达式,根据向量图中转矩工作区域大小随永磁和电枢磁链变化情况,同时结合表达式分析永磁和电枢磁场谐波所对应磁链对转矩和功率因数的影响,建立出永磁-电枢双谐波协同优化设计思路,提升电机的转矩密度和功率因数。Step 1: According to the harmonic characteristics of the permanent magnet and armature magnetic fields, respectively establish the expressions between the flux linkage and the harmonics of the permanent magnet and armature magnetic fields; according to the magnetic field modulation permanent magnet motor flux vector diagram (Figure 2), establish The expression of the power factor on the permanent magnet and armature flux linkage, according to the change of the size of the torque working area with the permanent magnet and armature flux linkage in the vector diagram, combined with the expression to analyze the corresponding magnetic field of the permanent magnet and armature magnetic field harmonics Based on the influence of the chain on the torque and power factor, the permanent magnet-armature double harmonic collaborative optimization design idea is established to improve the torque density and power factor of the motor.
对于永磁磁场谐波而言,从永磁磁场谐波对于转矩是否有贡献的角度将永磁磁场谐波分为永磁磁场工作谐波和永磁磁场非工作谐波两类。对转矩有贡献的永磁磁场谐波为永磁磁场工作谐波,对转矩没有贡献的永磁磁场谐波为永磁磁场非工作谐波。对各次永磁磁场工作谐波进行加权可得到合成永磁磁场工作谐波,合成永磁磁场工作谐波所对应的永磁磁链与转矩和功率因数都成正比,因此提高合成永磁磁场工作谐波幅值可在提高功率因数的同时维持较高转矩密度;对于电枢磁场谐波而言,从电枢磁场谐波对于转矩是否有贡献的角度将电枢磁场谐波分为电枢磁场工作谐波和电枢磁场非工作谐波两类。对转矩有贡献的电枢磁场谐波为电枢磁场工作谐波,对转矩没有贡献的电枢磁场谐波为电枢磁场非工作谐波。根据磁链和功率因数表达式,电枢磁场非工作谐波所对应的磁链与功率因数成反比,同时根据磁链相量图,降低电枢磁场非工作谐波不会影响转矩所对应的工作区域大小,因此减少电枢磁场非工作谐波可以提升功率因数,同时不损失转矩密度。For the permanent magnetic field harmonics, from the perspective of whether the permanent magnetic field harmonics contribute to the torque, the permanent magnetic field harmonics are divided into two types: permanent magnetic field working harmonics and permanent magnetic field non-working harmonics. The harmonics of the permanent magnetic field that contribute to the torque are the working harmonics of the permanent magnetic field, and the harmonics of the permanent magnetic field that do not contribute to the torque are the non-working harmonics of the permanent magnetic field. The working harmonics of the synthetic permanent magnetic field can be obtained by weighting the working harmonics of the permanent magnetic field. The working harmonic amplitude of the magnetic field can maintain a high torque density while improving the power factor; for the armature magnetic field harmonics, the armature magnetic field There are two types of armature magnetic field working harmonics and armature magnetic field non-working harmonics. The armature magnetic field harmonics that contribute to the torque are the armature magnetic field working harmonics, and the armature magnetic field harmonics that do not contribute to the torque are the armature magnetic field non-working harmonics. According to the flux linkage and power factor expressions, the flux linkage corresponding to the non-working harmonic of the armature magnetic field is inversely proportional to the power factor, and according to the flux linkage phasor diagram, reducing the non-working harmonic of the armature magnetic field will not affect the torque corresponding to The size of the working area, so reducing the non-working harmonics of the armature field can improve the power factor without losing torque density.
根据针对电机向量图和功率因数表达式的分析,确定出在保证合成永磁磁场工作谐波处于较高水平的前提下,通过对电枢磁场非工作谐波最小值和电枢磁场工作谐波最大值的优化,即可实现对永磁-电枢磁场双谐波的协同优化设计,从而提升电机的转矩密度和功率因数。According to the analysis of the vector diagram of the motor and the power factor expression, it is determined that under the premise of ensuring that the working harmonic of the synthetic permanent magnet magnetic field is at a high level, the minimum value of the non-working harmonic of the armature magnetic field and the working harmonic of the armature magnetic field are determined. The optimization of the maximum value can realize the collaborative optimization design of the double harmonics of the permanent magnet-armature magnetic field, thereby improving the torque density and power factor of the motor.
其中,步骤1中永磁磁场谐波特征公式B m(θ,t)和电枢磁场谐波特征公式B a(θ,t)的表达式为: Among them, the expressions of the permanent magnet magnetic field harmonic characteristic formula B m (θ, t) and the armature magnetic field harmonic characteristic formula B a (θ, t) in step 1 are:
Figure PCTCN2022070687-appb-000021
Figure PCTCN2022070687-appb-000021
Figure PCTCN2022070687-appb-000022
Figure PCTCN2022070687-appb-000022
式中,C m是永磁磁动势傅里叶系数,D i和D j是电枢磁动势傅里叶系数,m是永磁磁动势阶次,k是磁导阶次,i和j为电枢磁动势阶次,P r是永磁体极对数,Ω r是电机机械转速,t为时间,Λ 0和Λ k是气隙磁导傅里叶系数,N s是电枢槽数。根据表达式可确定出永磁磁场谐波阶次为mP r,mP r±kN s,电枢磁场谐波阶次为i,j,i±kN s,j±kN s。F m(θ,t)为永磁磁动势表达式,Λ s(θ)为气隙磁导表达式,分别可以表示为: In the formula, C m is the Fourier coefficient of the permanent magnetomotive force, D i and D j are the Fourier coefficients of the armature magnetomotive force, m is the order of the permanent magnetomotive force, k is the order of the permeability, i and j are the order of the magnetomotive force of the armature, P r is the number of pole pairs of the permanent magnet, Ω r is the mechanical speed of the motor, t is the time, Λ 0 and Λ k are the Fourier coefficients of the air gap permeance, N s is the electrical Number of pivot slots. According to the expression, it can be determined that the harmonic order of the permanent magnet magnetic field is mP r , mP r ±kN s , and the harmonic order of the armature magnetic field is i, j, i±kN s , j±kN s . F m (θ,t) is the permanent magnet magnetomotive force expression, Λ s (θ) is the air gap permeance expression, which can be expressed as:
Figure PCTCN2022070687-appb-000023
Figure PCTCN2022070687-appb-000023
Figure PCTCN2022070687-appb-000024
Figure PCTCN2022070687-appb-000024
进一步,步骤1中永磁磁链关于永磁磁场谐波的表达式为:Further, the expression of the permanent magnet flux linkage on the harmonics of the permanent magnet field in step 1 is:
Figure PCTCN2022070687-appb-000025
Figure PCTCN2022070687-appb-000025
式中,r g为气隙半径,l ef为轴向长度,n c为绕组匝数。其中,基波永磁磁链幅值表达式为: In the formula, r g is the air gap radius, l ef is the axial length, n c is the number of turns of the winding. Among them, the amplitude expression of the fundamental permanent magnet flux linkage is:
Figure PCTCN2022070687-appb-000026
Figure PCTCN2022070687-appb-000026
式中,C 1是永磁磁动势基波傅里叶系数。 In the formula, C 1 is the Fourier coefficient of the fundamental wave of the permanent magnetomotive force.
进一步,步骤1中电枢磁链关于电枢磁场谐波的表达式为:Further, the expression of the armature flux linkage with respect to the harmonics of the armature magnetic field in step 1 is:
Figure PCTCN2022070687-appb-000027
Figure PCTCN2022070687-appb-000027
其中,基波电枢磁链幅值表达式为:Among them, the amplitude expression of the fundamental armature flux linkage is:
Figure PCTCN2022070687-appb-000028
Figure PCTCN2022070687-appb-000028
进一步,步骤1中功率因数关于永磁和电枢磁链的表达式为:Further, the expression of the power factor with respect to the permanent magnet and armature flux linkage in step 1 is:
Figure PCTCN2022070687-appb-000029
Figure PCTCN2022070687-appb-000029
式中,
Figure PCTCN2022070687-appb-000030
为漏磁,
Figure PCTCN2022070687-appb-000031
为永磁磁链,
Figure PCTCN2022070687-appb-000032
为电枢磁场工作谐波对应磁链,
Figure PCTCN2022070687-appb-000033
为电枢磁场非工作谐波对应磁链,U和ω r分别为相电压和频率,E 0为永磁反电势。
In the formula,
Figure PCTCN2022070687-appb-000030
For flux leakage,
Figure PCTCN2022070687-appb-000031
is the permanent magnet flux linkage,
Figure PCTCN2022070687-appb-000032
The corresponding flux linkage for the working harmonic of the armature magnetic field,
Figure PCTCN2022070687-appb-000033
is the flux linkage corresponding to the non-working harmonic of the armature magnetic field, U and ω r are the phase voltage and frequency respectively, and E 0 is the permanent magnet back EMF.
步骤2:从永磁磁场的角度对合成永磁磁场工作谐波幅值进行限制。将合成永磁磁场工作谐波幅值的最小值作为约束条件,通过敏感度分析选择出对永磁磁场工作谐波影响较大的设计参数,基于对永磁磁场工作谐波具有高敏感度的设计参数,建立合成永磁磁场工作谐波的优化模型,利用克里格模型表示合成永磁磁场工作谐波幅值与对永磁磁场工作谐波具有高敏感度的设计参数之间的关系,根据优化模型中设定的最小合成永磁磁场工作谐波幅值,基于所建立的克里格模型得出对应的对永磁磁场工作谐波具有高敏感度的设计参数的变化范围。Step 2: Limit the working harmonic amplitude of the synthetic permanent magnetic field from the perspective of the permanent magnetic field. Taking the minimum value of the working harmonic amplitude of the synthetic permanent magnetic field as a constraint condition, the design parameters that have a greater impact on the working harmonics of the permanent magnetic field are selected through sensitivity analysis. Based on the high sensitivity to the working harmonics of the permanent magnetic field Design parameters, establish an optimization model for the working harmonics of the synthetic permanent magnetic field, and use the Kriging model to express the relationship between the working harmonic amplitude of the synthetic permanent magnetic field and the design parameters with high sensitivity to the working harmonics of the permanent magnetic field, According to the minimum synthetic permanent magnet magnetic field working harmonic amplitude set in the optimization model, based on the established Kriging model, the corresponding variation range of the design parameters with high sensitivity to the permanent magnetic magnetic field working harmonic is obtained.
进一步,选择一种磁场调制永磁电机作为本优化设计的实施对象(图3),步骤2中敏感度计算公式表达式为:Furthermore, a magnetic field modulation permanent magnet motor is selected as the implementation object of this optimization design (Fig. 3), and the sensitivity calculation formula in step 2 is as follows:
Figure PCTCN2022070687-appb-000034
Figure PCTCN2022070687-appb-000034
式中,Y(x)表示的是不同设计参数下永磁和电枢磁场谐波幅值,N为采样数目,x为电机设计参数,包括h m、w m、w p、w s、w a、h b和w bIn the formula, Y(x) represents the harmonic amplitude of the permanent magnet and armature magnetic field under different design parameters, N is the sampling number, and x is the motor design parameters, including h m , w m , w p , w s , w a , h b , and w b .
步骤2中敏感度分析的具体步骤在于:首先,运用中心复合设计采样方法,对满足二阶回归旋转准则的二阶因子设计点,轴点以及零水平中心点进行采样,为进行敏感度计算提供数据支撑。然后,通过灵敏度公式计算不同设计参数对于永磁磁场工作谐波的敏感度,选择敏感度较大的设计参数作为限制合成永磁磁场工作谐波幅值的设计参数。The specific steps of the sensitivity analysis in step 2 are as follows: First, use the central composite design sampling method to sample the second-order factor design points, axis points and zero-level center points that meet the second-order regression rotation criterion, providing a basis for the sensitivity calculation. data support. Then, the sensitivity of different design parameters to the working harmonics of the permanent magnetic field is calculated through the sensitivity formula, and the design parameters with greater sensitivity are selected as the design parameters that limit the amplitude of the synthetic permanent magnetic field working harmonics.
进行敏感度分析后,根据永磁磁场工作谐波对应参数的敏感度分析结果(图4),选择敏感度高的设计参数建立的合成永磁磁场工作谐波的优化模型其表达式为:After the sensitivity analysis, according to the sensitivity analysis results of the parameters corresponding to the working harmonics of the permanent magnetic field (Fig. 4), the optimized model of the synthetic permanent magnetic field working harmonics established by selecting the design parameters with high sensitivity is expressed as:
Figure PCTCN2022070687-appb-000035
Figure PCTCN2022070687-appb-000035
Constraint:H sm(x 2)>g 1 Constraint:H sm (x 2 )>g 1
x 2∈{w s,w p,w m,h m} x 2 ∈{w s ,w p ,w m ,h m }
式中,H sm为电机合成永磁磁场工作谐波幅值,x 2表示的是对电机永磁磁场谐波具有高敏感度的设计参数,H mh为第h次永磁磁场工作谐波幅值,a h为对应永磁磁场工作谐波的加权系数,g 1为合成永磁磁场工作谐波的最小约束值。图5展示了本发明实施例磁场调制永磁电机的合成永磁磁场工作谐波的克里格模型。 In the formula, H sm is the harmonic amplitude of the synthetic permanent magnet magnetic field of the motor, x2 represents the design parameter with high sensitivity to the harmonic of the permanent magnet magnetic field of the motor, and H mh is the working harmonic amplitude of the hth permanent magnet magnetic field value, a h is the weighting coefficient corresponding to the working harmonic of the permanent magnetic field, g 1 is the minimum constraint value of the synthetic permanent magnetic field working harmonic. Fig. 5 shows the Kriging model of the working harmonics of the synthesized permanent magnet magnetic field of the magnetic field modulated permanent magnet motor according to the embodiment of the present invention.
步骤3:简化电枢磁场谐波优化目标。由于进行优化的电枢磁场谐波阶次较多,需要对电枢磁场谐波优化目标进行简化,降低优化目标个数。首先,分析电枢磁场谐波关于电机性能的敏感度,选择敏感度较大的电枢磁场谐波作为优化目标,进一步根据实验设计法计算出电枢磁场谐波的分布,通过电枢磁场谐波实验点分布图,选择出与电枢磁场工作谐波变化趋势不同的电枢磁场非工作谐波为需要优化的谐波,并采用加权的方式得到简化的电枢磁场谐波优化目标。Step 3: Simplify the armature field harmonics optimization objective. Since there are many harmonic orders of the armature magnetic field to be optimized, it is necessary to simplify the optimization objectives of the armature magnetic field harmonics and reduce the number of optimization objectives. First, analyze the sensitivity of the armature magnetic field harmonics to the performance of the motor, select the more sensitive armature magnetic field harmonics as the optimization target, and further calculate the distribution of the armature magnetic field harmonics according to the experimental design method. According to the distribution diagram of wave experiment points, the non-operating harmonics of the armature magnetic field, which are different from the changing trend of the working harmonics of the armature magnetic field, are selected as the harmonics to be optimized, and the simplified optimization target of the armature magnetic field harmonics is obtained by weighting.
进一步,步骤3中简化后电枢磁场谐波优化目标由敏感度分析和电枢磁场谐波实验点分布得出。首先,运用中心复合设计采样方法进行采样,为敏感度计算提供数据支撑。然后,通过灵敏度公式计算不同电枢磁场谐波对于功率因数的敏感度(图6(a)),选择敏感度较大的磁场谐波作为电枢谐波优化目标。进一步根据复合设计采样方法中的实验点,绘制电枢磁场谐波实验点分布图(图6(b)),如图6(b)所示,29次电枢磁场非工作谐波与31次电枢磁场工作谐波具有相同的变化趋势,因此29次电枢磁场非工作谐波不作为优化目标,同时选择其余电枢磁场非工作谐波为需要优化的谐波,然后将这些谐波进行线性加权作为优化目标,减少电枢磁场谐波优化目标的个数,实现对电枢磁场谐波优化目标的简化。Further, the simplified optimization target of armature magnetic field harmonics in step 3 is obtained from the sensitivity analysis and the distribution of experimental point distribution of armature magnetic field harmonics. First, the central composite design sampling method is used for sampling to provide data support for the sensitivity calculation. Then, the sensitivity of different armature magnetic field harmonics to power factor is calculated by the sensitivity formula (Fig. 6(a)), and the magnetic field harmonic with greater sensitivity is selected as the target of armature harmonic optimization. Further, according to the experimental points in the composite design sampling method, draw the distribution diagram of the experimental points of the armature magnetic field harmonics (Fig. 6(b)), as shown in Fig. 6(b), the 29th non-working harmonic of the armature magnetic The working harmonics of the armature magnetic field have the same change trend, so the 29th non-working harmonics of the armature magnetic field are not used as the optimization target, and the rest of the non-working harmonics of the armature magnetic field are selected as the harmonics that need to be optimized, and then these harmonics are optimized Linear weighting is used as the optimization objective to reduce the number of harmonic optimization objectives of the armature magnetic field, and realize the simplification of the harmonic optimization objectives of the armature magnetic field.
步骤4:电枢磁场非工作谐波的独立性判断。通过计算电枢磁场谐波之间的交互效应,分析判断电枢磁场谐波中是否存在具有相对独立性的非工作谐波。如果存在相对独立性的 电枢磁场非工作谐波,进入步骤5.1和5.2;如果不存在相对独立性的电枢磁场非工作谐波,进入步骤5.3。图7分析了磁场调制永磁电机1次电枢磁场非工作波与其它电枢磁场谐波之间的交互效应。如图所示1次电枢磁场非工作谐波与其它谐波之间相互独立,因此1次电枢磁场非工作波是具有相对独立性的电枢磁场非工作谐波。Step 4: Independent judgment of non-working harmonics of armature magnetic field. By calculating the interaction effect between the armature magnetic field harmonics, it is analyzed and judged whether there are relatively independent non-working harmonics in the armature magnetic field harmonics. If there are relatively independent non-working harmonics of the armature magnetic field, go to steps 5.1 and 5.2; if there are no relatively independent non-working harmonics of the armature magnetic field, go to step 5.3. Figure 7 analyzes the interaction effect between the non-operating wave of the primary armature magnetic field and other harmonics of the armature magnetic field of the magnetic field modulated permanent magnet motor. As shown in the figure, the 1st non-working harmonic of the armature magnetic field is independent of other harmonics, so the 1st non-working wave of the armature magnetic field is a relatively independent non-working harmonic of the armature magnetic field.
步骤5.1:因为存在具有相对独立性的1次电枢磁场非工作谐波,将其与其余电枢磁场谐波分开进行优化,以减小设计参数和优化目标维度,提高电枢磁场谐波优化结果的准确性。首先,对具有相对独立性的1次电枢磁场非工作谐波进行优化,利用敏感度分析选择对具有相对独立性的电枢磁场1次非工作谐波敏感度较大的设计参数,如图8所示,设计参数h b和w b具有较大的敏感度,因此基于h b和w b建立设计参数关于具有相对独立性的1次电枢磁场非工作谐波的克里格模型。如图9所示,根据所建立的克里格模型选择出具有相对独立性的电枢磁场非工作谐波的最优设计点。 Step 5.1: Because there is a relatively independent armature magnetic field non-working harmonic, optimize it separately from the rest of the armature magnetic field harmonics to reduce the design parameters and optimize the target dimension, and improve the armature magnetic field harmonic optimization the accuracy of the results. Firstly, optimize the first-order non-working harmonic of the armature magnetic field with relative independence, and use the sensitivity analysis to select the design parameters that are more sensitive to the first-order non-operating harmonic of the armature magnetic field with relative independence, as shown in Fig. As shown in 8, the design parameters h b and w b have greater sensitivity, so based on h b and w b , the Kriging model of the design parameters with respect to the 1st order non-working harmonic of the armature magnetic field with relative independence is established. As shown in Figure 9, according to the established Kriging model, the optimal design point of non-working harmonics of the armature magnetic field with relative independence is selected.
其中,步骤5.1中具有相对独立性的电枢磁场非工作谐波的优化模型表达式为:Among them, the optimized model expression of the non-working harmonics of the armature magnetic field with relative independence in step 5.1 is:
Objectives:Min[λ 1H a1(x 3)] Objectives:Min[λ 1 H a1 (x 3 )]
式中,x 3为经过敏感度分析后对具有相对独立性的电枢磁场非工作谐波影响较大的设计参数h b和w b,H a1为第1次具有相对独立性的电枢磁场非工作谐波,λ 1为对应具有相对独立性的电枢磁场非工作谐波的加权系数,优化模型的最小值为优化目标。 In the formula, x 3 is the design parameters h b and w b that have a great influence on the non-operating harmonics of the relatively independent armature magnetic field after the sensitivity analysis, and H a1 is the first relatively independent armature magnetic field Non-working harmonics, λ 1 is the weighting coefficient corresponding to the relatively independent non-working harmonics of the armature magnetic field, and the minimum value of the optimization model is the optimization target.
步骤5.2:对于电枢磁场工作谐波与其余不具备相对独立性的电枢磁场的非工作谐波进行优化。将经过合成永磁磁场工作谐波幅值限制后的设计参数范围作为约束条件,并将电枢磁场工作谐波与经过简化后的电枢磁场非工作谐波作为优化目标,采用多目标遗传算法对电机电枢磁场工作谐波和非工作谐波进行优化,得出电枢谐波最优帕累托前沿图(图10),最终确定出在高合成永磁工作谐波幅值基础上具有最优电枢谐波的电机设计方案,实现对永磁-电枢磁场双谐波的协同优化设计,从而提升电机的转矩密度和功率因数。如图11所示,优化后的磁场调制永磁电机的转矩密度从26.2Nm/L提升至30.9Nm/L,功率因数从0.49提升至0.7,,验证了所提出优化设计方法的有效性。Step 5.2: Optimize the working harmonics of the armature magnetic field and other non-working harmonics of the armature magnetic field that do not have relative independence. Taking the design parameter range limited by the working harmonic amplitude of the synthetic permanent magnet magnetic field as the constraint condition, and taking the working harmonic of the armature magnetic field and the simplified non-working harmonic of the armature magnetic field as the optimization objectives, a multi-objective genetic algorithm is adopted The working harmonics and non-working harmonics of the motor armature magnetic field are optimized, and the optimal Pareto frontier diagram of the armature harmonics is obtained (Figure 10), and finally it is determined that there is The motor design scheme with optimal armature harmonics realizes the collaborative optimization design of permanent magnet-armature magnetic field double harmonics, thereby improving the torque density and power factor of the motor. As shown in Figure 11, the torque density of the optimized magnetic field modulation permanent magnet motor increases from 26.2Nm/L to 30.9Nm/L, and the power factor increases from 0.49 to 0.7, which verifies the effectiveness of the proposed optimal design method.
其中,步骤5.2中电枢磁场工作谐波与简化后的其余不具备相对独立性的电枢磁场非工作谐波的优化模型表达式为:Among them, the optimized model expression of the working harmonics of the armature magnetic field in step 5.2 and the simplified non-working harmonics of the armature magnetic field that do not have relative independence is:
Objectives:Max[H a31(x 4)],Min[μ 9H a9(x 4)+μ 11H a11(x 4)] Objectives:Max[H a31 (x 4 )],Min[μ 9 H a9 (x 4 )+μ 11 H a11 (x 4 )]
Constraint:f 1(x 2)>0 Constraint: f 1 (x 2 )>0
式中,H aPr为电机P r次电枢磁场工作谐波,其最大值被设定为优化目标,从而保证优化后电机具有较高转矩。μ 9H a9(x 4)+μ 11H a11(x 4)为经过简化后的电机合成电枢磁场非工作谐波,其 最小值被设定为优化目标。x 4表示的是在总设计参数x 1中排除掉设计参数x 3后所剩余的设计参数w s、w p、w m、w a、h m,H a9和H a11为第9、11次电枢磁场非工作谐波,μ 9和μ 11为对应第9、11次电枢磁场非工作谐波的加权系数。H a31为第31次电枢磁场工作谐波,μ 31为对应第31次电枢磁场工作谐波的加权系数。f 1(x 2)为经过合成永磁磁场工作谐波限制后的设计参数取值范围函数。 In the formula, H aPr is the working harmonic of the armature magnetic field of the motor P r , and its maximum value is set as the optimization target, so as to ensure that the optimized motor has a higher torque. μ 9 H a9 (x 4 )+μ 11 H a11 (x 4 ) is the simplified non-working harmonic of the synthesized armature magnetic field of the motor, and its minimum value is set as the optimization target. x 4 represents the remaining design parameters w s , w p , w m , w a , h m after excluding the design parameter x 3 from the total design parameters x 1, H a9 and H a11 are the 9th and 11th times The non-working harmonics of the armature magnetic field, μ 9 and μ 11 are the weighting coefficients corresponding to the 9th and 11th non-working harmonics of the armature magnetic field. H a31 is the 31st working harmonic of the armature magnetic field, and μ 31 is the weighting coefficient corresponding to the 31st working harmonic of the armature magnetic field. f 1 (x 2 ) is a function of the value range of the design parameters after the working harmonics of the synthetic permanent magnet magnetic field are limited.
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。In the description of this specification, references to the terms "one embodiment," "some embodiments," "exemplary embodiments," "example," "specific examples," or "some examples" are intended to mean that the implementation A specific feature, structure, material, or characteristic described by an embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。Although the embodiments of the present invention have been shown and described, those skilled in the art can understand that various changes, modifications, substitutions and modifications can be made to these embodiments without departing from the principle and spirit of the present invention. The scope of the invention is defined by the claims and their equivalents.

Claims (9)

  1. 一种磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,其特征在于,包括以下步骤:A permanent magnet-armature dual-harmonic collaborative optimization design method for a magnetic field modulated permanent magnet motor, characterized in that it comprises the following steps:
    步骤1:根据永磁磁场谐波特征公式和电枢磁场谐波特征公式,建立永磁磁链关于永磁磁场谐波之间的表达式,以及建立电枢磁链关于电枢磁场谐波之间的表达式;根据磁场调制永磁电机磁链向量图,建立功率因数关于永磁和电枢磁链的表达式,获得功率因数关于永磁和电枢磁场谐波的表达式,根据向量图中转矩工作区域大小随永磁和电枢磁链变化情况,同时结合以上表达式分析永磁和电枢磁场谐波所对应磁链对转矩和功率因数的影响,建立出永磁-电枢双谐波协同优化设计思路,提升电机的转矩密度和功率因数;Step 1: According to the permanent magnet magnetic field harmonic characteristic formula and the armature magnetic field harmonic characteristic formula, establish the expression between the permanent magnet flux linkage with respect to the permanent magnet magnetic field harmonics, and establish the relationship between the armature flux linkage with respect to the armature magnetic field harmonics The expression between; according to the magnetic field modulation permanent magnet motor flux vector diagram, establish the expression of the power factor about the permanent magnet and armature flux linkage, and obtain the expression of the power factor about the permanent magnet and armature magnetic field harmonics, according to the vector diagram The size of the medium torque working area varies with the flux linkage of the permanent magnet and the armature. At the same time, combining the above expressions, the influence of the flux linkage corresponding to the harmonics of the permanent magnet and armature magnetic field on the torque and power factor is analyzed, and the permanent magnet-electric The design idea of double-harmonic synergistic optimization improves the torque density and power factor of the motor;
    对于永磁磁场谐波而言,从永磁磁场谐波对于转矩是否有贡献的角度将永磁磁场谐波分为永磁磁场工作谐波和永磁磁场非工作谐波两类;对转矩有贡献的永磁磁场谐波为永磁磁场工作谐波,对转矩没有贡献的永磁磁场谐波为永磁磁场非工作谐波。对各次永磁磁场工作谐波进行加权可得到合成永磁磁场工作谐波,合成永磁磁场工作谐波所对应的永磁磁链与转矩和功率因数都成正比,因此提高合成永磁磁场工作谐波幅值可在提高功率因数的同时维持较高转矩密度;对于电枢磁场谐波而言,从电枢磁场谐波对于转矩是否有贡献的角度将电枢磁场谐波分为电枢磁场工作谐波和电枢磁场非工作谐波两类;对转矩有贡献的电枢磁场谐波为电枢磁场工作谐波,对转矩没有贡献的电枢磁场谐波为电枢磁场非工作谐波。根据磁链和功率因数表达式,电枢磁场非工作谐波所对应的磁链与功率因数成反比,同时根据磁链相量图,降低电枢磁场非工作谐波不会影响转矩所对应的工作区域大小,因此减少电枢磁场非工作谐波可以提升功率因数,同时不损失转矩密度;For the permanent magnetic field harmonics, from the perspective of whether the permanent magnetic field harmonics contribute to the torque, the permanent magnetic field harmonics are divided into two types: permanent magnetic field working harmonics and permanent magnetic field non-working harmonics; The harmonics of the permanent magnetic field that contribute to the torque are the working harmonics of the permanent magnetic field, and the harmonics of the permanent magnetic field that do not contribute to the torque are the non-working harmonics of the permanent magnetic field. The working harmonics of the synthetic permanent magnetic field can be obtained by weighting the working harmonics of the permanent magnetic field. The working harmonic amplitude of the magnetic field can maintain a high torque density while improving the power factor; for the armature magnetic field harmonics, the armature magnetic field There are two types of armature magnetic field working harmonics and armature magnetic field non-working harmonics; the armature magnetic field harmonics that contribute to the torque are the armature magnetic field working harmonics, and the armature magnetic field harmonics that do not contribute to the torque are the armature magnetic field harmonics. The pivot magnetic field is non-working harmonic. According to the flux linkage and power factor expressions, the flux linkage corresponding to the non-working harmonic of the armature magnetic field is inversely proportional to the power factor, and according to the flux linkage phasor diagram, reducing the non-working harmonic of the armature magnetic field will not affect the torque corresponding to The size of the working area, so reducing the non-working harmonics of the armature magnetic field can improve the power factor without losing the torque density;
    根据针对电机向量图和功率因数表达式的分析,确定出在保证合成永磁磁场工作谐波处于较高水平的前提下,通过对电枢磁场非工作谐波最小值和电枢磁场工作谐波最大值的优化,即可实现对永磁-电枢磁场双谐波的协同优化设计,从而提升电机的转矩密度和功率因数;According to the analysis of the vector diagram of the motor and the power factor expression, it is determined that under the premise of ensuring that the working harmonic of the synthetic permanent magnet magnetic field is at a high level, the minimum value of the non-working harmonic of the armature magnetic field and the working harmonic of the armature magnetic field are determined. The optimization of the maximum value can realize the collaborative optimization design of the double harmonics of the permanent magnet-armature magnetic field, thereby improving the torque density and power factor of the motor;
    步骤2:从永磁磁场的角度对合成永磁磁场工作谐波幅值进行限制,将合成永磁磁场工作谐波幅值的最小值作为约束条件,通过敏感度分析选择出对永磁磁场工作谐波影响较大的设计参数,基于对永磁磁场工作谐波具有高敏感度的设计参数,建立合成永磁磁场工作谐波的优化模型,利用克里格模型表示合成永磁磁场工作谐波幅值与对永磁磁场工作谐波具有高敏感度的设计参数之间的关系,根据优化模型中设定的最小合成永磁磁场工作谐波幅值,基于所建立的克里格模型得出对应的对永磁磁场工作谐波具有高敏感度的设计参数的变化范围;Step 2: Limit the working harmonic amplitude of the synthetic permanent magnetic field from the perspective of the permanent magnetic field, take the minimum value of the working harmonic amplitude of the synthetic permanent magnetic field as a constraint condition, and select the working harmonic amplitude of the permanent magnetic field through sensitivity analysis Harmonics have a greater influence on the design parameters, based on the design parameters with high sensitivity to the permanent magnetic field working harmonics, an optimization model for the synthetic permanent magnetic field working harmonics is established, and the Kriging model is used to represent the synthetic permanent magnetic field working harmonics The relationship between the amplitude and the design parameters with high sensitivity to the working harmonics of the permanent magnetic field, according to the minimum synthetic permanent magnetic field working harmonic amplitude set in the optimization model, is obtained based on the established Kriging model Corresponding to the variation range of design parameters with high sensitivity to the working harmonics of the permanent magnetic field;
    步骤3:简化电枢磁场谐波优化目标,由于进行优化的电枢磁场谐波阶次较多,需要对电枢磁场谐波优化目标进行简化,降低优化目标个数;首先,分析电枢磁场谐波关于电机性能的敏感度,选择敏感度较大的电枢磁场谐波作为优化目标,进一步根据实验设计法计算出电枢磁场谐波的分布,通过电枢磁场谐波实验点分布图,选择出电枢磁场工作谐波变化趋势不同的电枢磁场非工作谐波为需要优化的谐波,并采用加权的方式得到简化的电枢磁场谐波优化目标;Step 3: Simplify the optimization target of the armature magnetic field harmonics. Since there are many harmonic orders of the armature magnetic field to be optimized, it is necessary to simplify the optimization target of the armature magnetic field harmonics and reduce the number of optimization targets; first, analyze the armature magnetic field The sensitivity of harmonics to motor performance, select the more sensitive armature magnetic field harmonics as the optimization target, and further calculate the distribution of armature magnetic field harmonics according to the experimental design method, through the armature magnetic field harmonic experimental point distribution map, The non-working harmonics of the armature magnetic field with different changing trends of the working harmonics of the armature magnetic field are selected as the harmonics to be optimized, and the simplified optimization target of the armature magnetic field harmonics is obtained by weighting;
    步骤4:电枢磁场非工作谐波的独立性判断。通过计算电枢磁场谐波之间的交互效应,分析判断电枢磁场谐波中是否存在具有相对独立性的非工作谐波。如果存在相对独立性的电枢磁场非工作谐波,进入步骤5.1和5.2;如果不存在相对独立性的电枢磁场非工作谐波,进入步骤5.3Step 4: Independent judgment of non-working harmonics of armature magnetic field. By calculating the interaction effect between the armature magnetic field harmonics, it is analyzed and judged whether there are relatively independent non-working harmonics in the armature magnetic field harmonics. If there are relatively independent non-working harmonics of the armature magnetic field, go to steps 5.1 and 5.2; if there are no relatively independent non-working harmonics of the armature magnetic field, go to step 5.3
    步骤5.1:如果存在具有相对独立性的电枢磁场非工作谐波,将具有相对独立性的电枢磁场非工作谐波与其余电枢磁场谐波分开进行优化,以减小设计参数和优化目标维度,提高电枢磁场谐波优化结果的准确性;首先,对具有相对独立性的电枢磁场非工作谐波进行优化,利用敏感度分析选择对具有相对独立性的电枢磁场非工作谐波敏感度较大的设计参数,建立设计参数关于具有相对独立性的电枢磁场非工作谐波的克里格模型,根据所建立的克里格模型选择出具有相对独立性的电枢磁场非工作谐波的最优设计点;Step 5.1: If there are non-operating harmonics of the armature field with relative independence, optimize the non-operating harmonics of the armature field with relative independence separately from the rest of the harmonics of the armature field to reduce the design parameters and optimization objectives dimension to improve the accuracy of the optimization results of armature magnetic field harmonics; firstly, optimize the relatively independent armature magnetic field non-working For design parameters with high sensitivity, establish the Kriging model of the design parameters about the non-working harmonic of the armature magnetic field with relative independence, and select the non-working armature magnetic field with relative independence according to the established Kriging model The optimal design point for harmonics;
    步骤5.2:对于电枢磁场工作谐波与其余不具备相对独立性的电枢磁场的非工作谐波进行优化,将经过合成永磁磁场工作谐波幅值限制后的设计参数范围作为约束条件,并将电枢磁场工作谐波与经过简化后的电枢磁场非工作谐波作为优化目标,采用多目标遗传算法对电机电枢磁场工作谐波和非工作谐波进行优化,最终确定出在高合成永磁工作谐波幅值基础上具有最优电枢谐波的电机设计方案,实现对永磁-电枢磁场双谐波的协同优化设计,从而提升电机的转矩密度和功率因数;Step 5.2: optimize the working harmonics of the armature magnetic field and other non-working harmonics of the armature magnetic field that do not have relative independence, and use the range of design parameters limited by the amplitude of the working harmonics of the synthesized permanent magnetic field as a constraint condition, Taking the working harmonics of the armature magnetic field and the simplified non-working harmonics of the armature magnetic field as optimization targets, the multi-objective genetic algorithm is used to optimize the working harmonics and non-working harmonics of the armature magnetic field of the motor, and finally determine the The motor design scheme with the optimal armature harmonics based on the synthetic permanent magnet working harmonic amplitude, realizes the collaborative optimization design of the permanent magnet-armature magnetic field double harmonics, thereby improving the torque density and power factor of the motor;
    步骤5.3:对于电枢磁场的工作谐波及非工作谐波进行优化。将经过合成永磁磁场工作谐波幅值限制后的参数范围作为约束条件,将电枢磁场工作谐波与经过简化后的电枢磁场非工作谐波作为优化目标,采用多目标遗传算法对电机电枢磁场工作谐波和非工作谐波进行优化,最终确定出在高合成永磁工作谐波幅值基础上具有最优电枢谐波的电机设计方案,实现对永磁-电枢磁场双谐波的协同优化设计,从而提升电机的转矩密度和功率因数。Step 5.3: Optimize the working harmonics and non-working harmonics of the armature magnetic field. Taking the parameter range limited by the working harmonic amplitude of the synthesized permanent magnet magnetic field as the constraint condition, taking the working harmonic of the armature magnetic field and the simplified non-working harmonic of the armature magnetic field as the optimization target, the multi-objective genetic algorithm is used to optimize the The working harmonics and non-working harmonics of the armature magnetic field of the machine are optimized, and finally the design scheme of the motor with the optimal armature harmonics is determined on the basis of the high synthetic permanent magnet working harmonic amplitude. The coordinated optimization design of harmonics can improve the torque density and power factor of the motor.
  2. 根据权利要求1所述磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,其特征在于,步骤1中永磁磁场谐波特征公式B m(θ,t)和电枢磁场谐波特征公式B a(θ,t)为: According to the magnetic field modulation permanent magnet motor permanent magnet-armature double harmonic collaborative optimization design method of claim 1, it is characterized in that, in step 1, the permanent magnet magnetic field harmonic characteristic formula B m (θ, t) and the armature magnetic field harmonic The wave characteristic formula B a (θ,t) is:
    Figure PCTCN2022070687-appb-100001
    Figure PCTCN2022070687-appb-100001
    Figure PCTCN2022070687-appb-100002
    Figure PCTCN2022070687-appb-100002
    式中,C m是永磁磁动势傅里叶系数,D i和D j是电枢磁动势傅里叶系数,m是永磁磁动势阶次,k是磁导阶次,i和j为电枢磁动势阶次,P r是永磁体极对数,Ω r是电机机械转速,t为时间,Λ 0和Λ k是气隙磁导傅里叶系数,N s是电枢槽数,根据表达式可确定出永磁磁场谐波阶次为mP r,mP r±kN s,电枢磁场谐波阶次为i,j,i±kN s,j±kN s,F a(θ,t)为电枢磁动势;F m(θ,t)为永磁磁动势表达式,Λ s(θ)为气隙磁导表达式,分别可以表示为: In the formula, C m is the Fourier coefficient of the permanent magnetomotive force, D i and D j are the Fourier coefficients of the armature magnetomotive force, m is the order of the permanent magnetomotive force, k is the order of the permeability, i and j are the order of the magnetomotive force of the armature, P r is the number of pole pairs of the permanent magnet, Ω r is the mechanical speed of the motor, t is the time, Λ 0 and Λ k are the Fourier coefficients of the air gap permeance, N s is the electrical According to the number of pivot slots, the harmonic order of the permanent magnet magnetic field can be determined as mP r , mP r ±kN s , and the harmonic order of the armature magnetic field is i, j, i±kN s , j±kN s , F a (θ,t) is the armature magnetomotive force; F m (θ,t) is the permanent magnet magnetomotive force expression, Λ s (θ) is the air gap permeance expression, which can be expressed as:
    Figure PCTCN2022070687-appb-100003
    Figure PCTCN2022070687-appb-100003
    Figure PCTCN2022070687-appb-100004
    Figure PCTCN2022070687-appb-100004
  3. 根据权利要求1所述磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,其特征在于,步骤1中永磁磁链关于永磁磁场谐波的表达式为:According to claim 1, the permanent magnet-armature double harmonic collaborative optimization design method for magnetic field modulation permanent magnet motor is characterized in that, in step 1, the expression of the permanent magnet flux linkage with respect to the permanent magnet magnetic field harmonic is:
    Figure PCTCN2022070687-appb-100005
    Figure PCTCN2022070687-appb-100005
    式中,r g为气隙半径,l ef为轴向长度,n c为绕组匝数;其中,基波永磁磁链幅值表达式为: In the formula, r g is the radius of the air gap, l ef is the axial length, n c is the number of turns of the winding; where, the expression of the amplitude of the fundamental permanent magnet flux linkage is:
    Figure PCTCN2022070687-appb-100006
    Figure PCTCN2022070687-appb-100006
    式中,C 1是永磁磁动势基波傅里叶系数。 In the formula, C 1 is the Fourier coefficient of the fundamental wave of the permanent magnetomotive force.
  4. 根据权利要求1所述磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,其特征在于,步骤1中电枢磁链关于电枢磁场谐波的表达式为:According to claim 1, the permanent magnet-armature double harmonic collaborative optimization design method for magnetic field modulation permanent magnet motor is characterized in that the expression of armature flux linkage in step 1 with respect to armature magnetic field harmonics is:
    Figure PCTCN2022070687-appb-100007
    Figure PCTCN2022070687-appb-100007
    其中,基波电枢磁链幅值表达式为:Among them, the amplitude expression of the fundamental armature flux linkage is:
    Figure PCTCN2022070687-appb-100008
    Figure PCTCN2022070687-appb-100008
  5. 根据权利要求1所述磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,其特征在于,步骤1中功率因数关于永磁和电枢磁链的表达式为:According to claim 1, the permanent magnet-armature dual-harmonic collaborative optimization design method for magnetic field modulation permanent magnet motor is characterized in that, in step 1, the expression of the power factor with respect to the permanent magnet and armature flux linkage is:
    Figure PCTCN2022070687-appb-100009
    Figure PCTCN2022070687-appb-100009
    式中,
    Figure PCTCN2022070687-appb-100010
    为漏磁,
    Figure PCTCN2022070687-appb-100011
    为永磁磁链,
    Figure PCTCN2022070687-appb-100012
    为电枢磁场工作谐波对应磁链,
    Figure PCTCN2022070687-appb-100013
    为电枢磁场非工作谐波对应磁链,U和ω r分别为相电压和频率,E 0为永磁反电势。
    In the formula,
    Figure PCTCN2022070687-appb-100010
    For flux leakage,
    Figure PCTCN2022070687-appb-100011
    is the permanent magnet flux linkage,
    Figure PCTCN2022070687-appb-100012
    The corresponding flux linkage for the working harmonic of the armature magnetic field,
    Figure PCTCN2022070687-appb-100013
    is the flux linkage corresponding to the non-working harmonic of the armature magnetic field, U and ω r are the phase voltage and frequency respectively, and E 0 is the back EMF of the permanent magnet.
  6. 根据权利要求1所述磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,其特征在于,步骤2中敏感度计算公式表达式为:According to claim 1, the permanent magnet-armature dual-harmonic collaborative optimization design method for magnetic field modulation permanent magnet motor is characterized in that, the sensitivity calculation formula in step 2 is expressed as:
    Figure PCTCN2022070687-appb-100014
    Figure PCTCN2022070687-appb-100014
    式中,Y(x)表示的是不同设计参数下永磁和电枢磁场谐波幅值,N为采样数目,x为电机设计参数;In the formula, Y(x) represents the harmonic amplitude of the permanent magnet and armature magnetic field under different design parameters, N is the sampling number, and x is the motor design parameter;
    步骤2中敏感度分析的具体步骤在于:首先,运用中心复合设计采样方法,对满足二阶回归旋转准则的二阶因子设计点,轴点以及零水平中心点进行采样,为进行敏感度计算提供数据支撑;然后,通过灵敏度公式计算不同设计参数对于永磁磁场工作谐波的敏感度,选择敏 感度较大的设计参数作为限制合成永磁磁场工作谐波幅值的设计参数;The specific steps of the sensitivity analysis in step 2 are as follows: First, use the central composite design sampling method to sample the second-order factor design points, axis points and zero-level center points that meet the second-order regression rotation criterion, providing a basis for the sensitivity calculation. Data support; then, the sensitivity of different design parameters to the working harmonics of the permanent magnetic field is calculated through the sensitivity formula, and the design parameters with greater sensitivity are selected as the design parameters that limit the amplitude of the synthetic permanent magnetic field working harmonics;
    进行敏感度分析后建立的合成永磁磁场工作谐波的优化模型的表达式为:The expression of the optimization model of the working harmonic of the synthetic permanent magnetic field established after the sensitivity analysis is:
    Figure PCTCN2022070687-appb-100015
    Figure PCTCN2022070687-appb-100015
    Constraint:H sm(x 2)>g 1 Constraint:H sm (x 2 )>g 1
    式中,H sm为电机合成永磁磁场工作谐波幅值,x 2表示的是对电机永磁磁场谐波具有高敏感度的设计参数,H mh为第h次永磁磁场工作谐波幅值,a h为对应永磁磁场工作谐波的加权系数,g 1为合成永磁磁场工作谐波的最小约束值。 In the formula, H sm is the harmonic amplitude of the synthetic permanent magnet magnetic field of the motor, x2 represents the design parameter with high sensitivity to the harmonic of the permanent magnet magnetic field of the motor, and H mh is the working harmonic amplitude of the hth permanent magnet magnetic field value, a h is the weighting coefficient corresponding to the working harmonic of the permanent magnetic field, g 1 is the minimum constraint value of the synthetic permanent magnetic field working harmonic.
  7. 根据权利要求1所述磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,其特征在于,步骤3中,简化后电枢磁场谐波优化目标由敏感度分析和电枢磁场谐波实验点分布得出;首先,运用中心复合设计采样方法进行采样,为敏感度计算提供数据支撑;然后,通过灵敏度公式计算不同电枢磁场谐波对于功率因数的敏感度,选择敏感度较大的磁场谐波作为电枢谐波优化目标;进一步根据复合设计采样方法中的实验点,绘制电枢磁场谐波实验点分布图,根据图中不同电枢磁场谐波的变化趋势,选择出与电枢磁场工作谐波变化趋势不同的电枢磁场非工作谐波为需要优化的谐波,而与电枢磁场工作谐波变化趋势相同的电枢磁场非工作谐波不作为需要优化的谐波;然后将这些需要优化的谐波进行线性加权作为优化目标,减少电枢磁场谐波优化目标的个数,实现对电枢磁场谐波优化目标的简化。According to the magnetic field modulation permanent magnet motor permanent magnet-armature double harmonic collaborative optimization design method of claim 1, it is characterized in that, in step 3, the armature magnetic field harmonic optimization target after simplification is determined by sensitivity analysis and armature magnetic field harmonic First, use the central composite design sampling method to sample to provide data support for the sensitivity calculation; then, use the sensitivity formula to calculate the sensitivity of different armature magnetic field harmonics to the power factor, and choose a higher sensitivity The magnetic field harmonics of the armature are used as the optimization target of the armature harmonics; further, according to the experimental points in the composite design sampling method, the distribution map of the experimental point distribution of the armature magnetic field harmonics is drawn, and according to the change trend of different armature magnetic field harmonics in the figure, the The non-working harmonics of the armature magnetic field with different changing trends of the working harmonics of the armature magnetic field are the harmonics that need to be optimized, while the non-working harmonics of the armature magnetic field that have the same changing trend as the working harmonics of the armature magnetic field are not considered as the harmonics that need to be optimized ; Then these harmonics that need to be optimized are linearly weighted as optimization targets, reducing the number of optimization targets for armature magnetic field harmonics, and realizing the simplification of the optimization targets for armature magnetic field harmonics.
  8. 根据权利要求1所述磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,其特征在于,步骤5.1中具有相对独立性的电枢磁场非工作谐波的优化模型表达式为:According to claim 1, the permanent magnet-armature dual-harmonic collaborative optimization design method of magnetic field modulation permanent magnet motor is characterized in that, in step 5.1, the optimization model expression of the non-working harmonic of the armature magnetic field having relative independence is:
    Figure PCTCN2022070687-appb-100016
    Figure PCTCN2022070687-appb-100016
    式中,x 3为经过敏感度分析后对具有相对独立性的电枢磁场非工作谐波影响较大的设计参数,H al为第l次具有相对独立性的电枢磁场非工作谐波,λ l为对应具有相对独立性的电枢磁场非工作谐波的加权系数,优化模型的最小值为优化目标; In the formula, x3 is a design parameter that has a great influence on the non-operating harmonic of the armature magnetic field with relative independence after sensitivity analysis, and Hal is the non-operating harmonic of the armature magnetic field with relative independence for the lth time, λ l is the weighting coefficient corresponding to the non-working harmonic of the armature magnetic field with relative independence, and the minimum value of the optimization model is the optimization goal;
    步骤5.2中电枢磁场工作谐波与简化后的其余不具备相对独立性的电枢磁场非工作谐波的优化模型表达式为:In step 5.2, the optimized model expression of the working harmonics of the armature magnetic field and the simplified non-working harmonics of the armature magnetic field that do not have relative independence is:
    Figure PCTCN2022070687-appb-100017
    Figure PCTCN2022070687-appb-100017
    Constraint:f 1(x 2)>0 Constraint: f 1 (x 2 )>0
    式中,H aPr为电机P r次电枢磁场工作谐波,其最大值被设定为优化目标,从而保证优化后电机具有较高转矩,
    Figure PCTCN2022070687-appb-100018
    为经过简化后的电机合成电枢磁场非工作谐波,其最小值被设定为优化目标,x 4表示的是在总设计参数x 1中排除掉设计参数x 3后所剩余的设计参数,H as为第s 次电枢磁场非工作谐波,μ s为对应s次电枢磁场非工作谐波的加权系数,f 1(x 2)为经过合成永磁磁场工作谐波限制后的设计参数取值范围函数。
    In the formula, H aPr is the working harmonic of the armature magnetic field of the motor P r , and its maximum value is set as the optimization target, so as to ensure that the optimized motor has a higher torque,
    Figure PCTCN2022070687-appb-100018
    The non-working harmonic of the armature magnetic field is synthesized for the simplified motor, and its minimum value is set as the optimization target. x 4 represents the remaining design parameters after excluding the design parameter x 3 from the total design parameters x 1 . H as is the sth non-working harmonic of the armature magnetic field, μ s is the weighting coefficient corresponding to the s-th non-working armature magnetic field harmonic, and f 1 (x 2 ) is the design after limiting the working harmonic of the synthetic permanent magnet magnetic field Parameter value range function.
  9. 根据权利要求1所述磁场调制永磁电机永磁-电枢双谐波协同优化设计方法,其特征在于,步骤5.3中,电枢磁场工作谐波与经过简化后的电枢磁场非工作谐波的优化模型表达式为:According to claim 1, the permanent magnet-armature double harmonic collaborative optimization design method for magnetic field modulation permanent magnet motor is characterized in that, in step 5.3, the working harmonic of the armature magnetic field and the non-working harmonic of the simplified armature magnetic field The optimization model expression of is:
    Figure PCTCN2022070687-appb-100019
    Figure PCTCN2022070687-appb-100019
    Constraint:f 1(x 2)>0 Constraint: f 1 (x 2 )>0
    式中,H aPr为电机P r次电枢磁场工作谐波,其最大值被设定为优化目标,从而保证优化后电机具有较高转矩;
    Figure PCTCN2022070687-appb-100020
    为经过简化后的电机合成电枢磁场非工作谐波,其最小值被设定为优化目标;x 1为电机总的设计参数,x 2表示的是对电机永磁磁场谐波具有高敏感度的设计参数,H as为第s次电枢磁场非工作谐波,μ s为对应s次电枢磁场非工作谐波的加权系数;f 1(x 2)为经过合成永磁磁场工作谐波限制后的设计参数取值范围函数。
    In the formula, H aPr is the working harmonic of the armature magnetic field of the motor P r , and its maximum value is set as the optimization target, so as to ensure that the optimized motor has a higher torque;
    Figure PCTCN2022070687-appb-100020
    The non-working harmonics of the armature magnetic field are synthesized for the simplified motor, and its minimum value is set as the optimization target; x 1 is the overall design parameter of the motor, and x 2 indicates that it is highly sensitive to the harmonics of the permanent magnet magnetic field of the motor H as is the sth non-working harmonic of the armature magnetic field, μ s is the weighting coefficient corresponding to the s-th non-working armature magnetic field harmonic; f 1 (x 2 ) is the working harmonic of the synthesized permanent magnet magnetic field Restricted design parameter value range function.
PCT/CN2022/070687 2021-12-07 2022-01-07 Method for collaborative optimization design of permanent magnet-armature double harmonics of magnetic field-modulated permanent magnet motor WO2023103138A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
GB2218908.8A GB2612711A (en) 2021-12-07 2022-01-07 Permanent magnet-armature double-harmonic collaborative optimization design method for magnetic field-modulated permanent magnet motor

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202111485800.3A CN114204711B (en) 2021-12-07 2021-12-07 Permanent magnet-armature double-harmonic collaborative optimization design method for magnetic field modulation permanent magnet motor
CN202111485800.3 2021-12-07

Publications (1)

Publication Number Publication Date
WO2023103138A1 true WO2023103138A1 (en) 2023-06-15

Family

ID=80650996

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/070687 WO2023103138A1 (en) 2021-12-07 2022-01-07 Method for collaborative optimization design of permanent magnet-armature double harmonics of magnetic field-modulated permanent magnet motor

Country Status (2)

Country Link
CN (1) CN114204711B (en)
WO (1) WO2023103138A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117610356A (en) * 2023-11-24 2024-02-27 哈尔滨工业大学 Permanent magnet motor cogging torque approximate modeling method based on improved vector superposition

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2612711A (en) * 2021-12-07 2023-05-10 Univ Jiangsu Permanent magnet-armature double-harmonic collaborative optimization design method for magnetic field-modulated permanent magnet motor
CN115276335B (en) * 2022-07-29 2023-11-10 江苏大学 Magnetic field modulation permanent magnet motor optimization method considering magnetic leakage and magnetic field harmonic wave
CN116861720B (en) * 2023-06-01 2024-04-05 青岛大学 Multi-objective layered optimization method for field modulation motor based on power factor expansion

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120091850A1 (en) * 2010-10-14 2012-04-19 Hitachi, Ltd. Rotating electric machine
CN105391265A (en) * 2015-12-21 2016-03-09 东南大学 Hybrid excitation fault-tolerant motor system based on brushless harmonic excitation
CN112366912A (en) * 2020-10-19 2021-02-12 江苏大学 Double-winding low-harmonic design method based on fractional-slot concentrated winding permanent magnet motor
CN112953046A (en) * 2021-03-17 2021-06-11 江苏大学 Stator modularization magnetic field modulation motor with low harmonic loss

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112421919B (en) * 2020-11-09 2021-10-08 华中科技大学 Permanent magnet motor topology construction method based on working magnetic field harmonic wave orientation and motor
CN112436706B (en) * 2020-11-24 2022-03-22 江苏大学 Loss analysis and suppression method for magnetic field modulation permanent magnet motor
CN113555986B (en) * 2021-06-23 2022-09-16 江苏大学 Magnetic field modulation type radial permanent magnet motor and multi-harmonic optimization design method thereof

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120091850A1 (en) * 2010-10-14 2012-04-19 Hitachi, Ltd. Rotating electric machine
CN105391265A (en) * 2015-12-21 2016-03-09 东南大学 Hybrid excitation fault-tolerant motor system based on brushless harmonic excitation
CN112366912A (en) * 2020-10-19 2021-02-12 江苏大学 Double-winding low-harmonic design method based on fractional-slot concentrated winding permanent magnet motor
CN112953046A (en) * 2021-03-17 2021-06-11 江苏大学 Stator modularization magnetic field modulation motor with low harmonic loss

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117610356A (en) * 2023-11-24 2024-02-27 哈尔滨工业大学 Permanent magnet motor cogging torque approximate modeling method based on improved vector superposition
CN117610356B (en) * 2023-11-24 2024-05-28 哈尔滨工业大学 Permanent magnet motor cogging torque approximate modeling method based on improved vector superposition

Also Published As

Publication number Publication date
CN114204711B (en) 2023-09-26
CN114204711A (en) 2022-03-18

Similar Documents

Publication Publication Date Title
WO2023103138A1 (en) Method for collaborative optimization design of permanent magnet-armature double harmonics of magnetic field-modulated permanent magnet motor
Gao et al. Improved genetic optimization algorithm with subdomain model for multi-objective optimal design of SPMSM
CN113094911B (en) High-power factor design method for magnetic field modulation permanent magnet fault-tolerant motor
CN109378918A (en) A kind of double permanent magnetism vernier motors of DC bias current type rotor
CN112737160B (en) Method for improving power factor of concentrated winding outer rotor vernier motor
WO2018214236A1 (en) Method for adjusting high efficiency region of permanent magnet synchronous motor
CN115276335B (en) Magnetic field modulation permanent magnet motor optimization method considering magnetic leakage and magnetic field harmonic wave
CN106295004B (en) Permanent magnet motor optimization design method considering disturbance type design variable interval
Tavana et al. Analysis and design of magnetic pole shape in linear permanent-magnet machine
Nazir et al. Investigation of sinusoidal shaped rotor to reduce torque ripple in axial flux permanent magnet machine
Nobahari et al. Modified Halbach design for torque profile improvement in a permanent magnet synchronous motor dedicated to an electric city bus
Hoang et al. Flux control range broadening and torque ripple minimization of a double excitation synchronous motor
Baek et al. Optimal design and comparison of stator winding configurations in permanent magnet assisted synchronous reluctance generator
Rahmani Fard et al. Dynamic performance of the novel axial flux-switching permanent magnet motor
Sardar et al. Design and Numerical Optimization for High-Performance IPM Synchronous Motor Using PSO Search Algorithm
Gao et al. A memetic-inspired harmony search method in optimal wind generator design
Mandasari et al. Design and optimization of brushless dc motor for electric boat thruster
CN116317233B (en) Surface-embedded permanent magnet type double-stator hybrid excitation motor, design analysis method thereof and performance optimization method of air gap field harmonic wave
Yang et al. Design and analysis of a new bearingless switched reluctance motor
Lee et al. Improvement of PMa-SynRM output characteristics by optimizing the rib structure
Song et al. Cogging torque optimization and analysis of hybrid stator permanent magnet motor
Li et al. Analysis and optimization of a novel dual PM vernier machine with halbach array
CN112398253B (en) Five-phase unequal tooth shoe double-harmonic permanent magnet synchronous motor and tooth shoe radian optimization method
Wi et al. Surrogate Assisted Contour Algorithm for Optimal Design of Interior Permanent Magnet Synchronous Motor for Electric Vehicles
Fuzhen et al. A novel spoke-type IPM machine with surface-mounted magnets to improve torque characteristics

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 202218908

Country of ref document: GB

Kind code of ref document: A

Free format text: PCT FILING DATE = 20220107

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22902583

Country of ref document: EP

Kind code of ref document: A1