CN110112979B - Permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit - Google Patents
Permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit Download PDFInfo
- Publication number
- CN110112979B CN110112979B CN201910399476.XA CN201910399476A CN110112979B CN 110112979 B CN110112979 B CN 110112979B CN 201910399476 A CN201910399476 A CN 201910399476A CN 110112979 B CN110112979 B CN 110112979B
- Authority
- CN
- China
- Prior art keywords
- permanent magnet
- magnet synchronous
- torque
- synchronous motor
- flux linkage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 93
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000004907 flux Effects 0.000 claims abstract description 83
- 239000013598 vector Substances 0.000 claims abstract description 53
- 238000005070 sampling Methods 0.000 claims abstract description 7
- 230000002457 bidirectional effect Effects 0.000 claims description 33
- 238000013178 mathematical model Methods 0.000 claims description 10
- 230000003068 static effect Effects 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 3
- OFHCOWSQAMBJIW-AVJTYSNKSA-N alfacalcidol Chemical compound C1(/[C@@H]2CC[C@@H]([C@]2(CCC1)C)[C@H](C)CCCC(C)C)=C\C=C1\C[C@@H](O)C[C@H](O)C1=C OFHCOWSQAMBJIW-AVJTYSNKSA-N 0.000 claims description 3
- 230000009466 transformation Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 9
- 238000013461 design Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000011217 control strategy Methods 0.000 description 3
- 230000006698 induction Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000002902 bimodal effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/20—Estimation of torque
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
- H02P25/02—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
- H02P25/022—Synchronous motors
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
- H02P6/14—Electronic commutators
- H02P6/16—Circuit arrangements for detecting position
- H02P6/18—Circuit arrangements for detecting position without separate position detecting elements
- H02P6/182—Circuit arrangements for detecting position without separate position detecting elements using back-emf in windings
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
The invention provides a permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit, which comprises the following steps: sampling three-phase current, angular speed and rotor position angle of the permanent magnet synchronous motor at the moment k, and obtaining output current and equivalent counter electromotive force through coordinate transformation; calculating an output voltage vector and a stator voltage according to the switching state of the voltage source inverter, and predicting an output current, a stator flux linkage amplitude and a torque at the moment of k +1 according to the stator voltage; calculating a torque target function, a stator flux linkage target function and respective corresponding standard deviations thereof by using the stator flux linkage amplitude and the torque, and further obtaining a new target function; and finally, using the voltage vector corresponding to the minimum objective function to control the permanent magnet synchronous motor. The method realizes the prediction torque control without the weight coefficient of the permanent magnet synchronous motor by performing standard deviation per unit on the torque and flux linkage target functions, simplifies the complexity of the system, reduces the torque ripple and improves the torque control precision.
Description
Technical Field
The invention relates to the field of power electronics, in particular to a permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit.
Background
In recent years, new energy electric vehicle technology has been developed vigorously to deal with energy crisis. Compared with an asynchronous Motor, a Permanent Magnet Synchronous Motor (PMSM) is widely applied to the field of electric vehicles due to the advantages of high efficiency, high power density and the like. In order to improve the torque dynamic control characteristic of the permanent magnet synchronous motor, the model prediction control mechanism is widely applied to the drive control of the permanent magnet synchronous motor. However, the conventional method for controlling the predicted torque of the permanent magnet synchronous motor has the disadvantage that the weight coefficient needs to be designed. Although the literature researches a permanent magnet synchronous motor non-weight coefficient prediction torque control method, no better weight coefficient theoretical design method exists so far.
The invention provides a rapid non-weight coefficient model predictive control method, which takes a target voltage vector as a target function to avoid using a weight coefficient, thereby realizing rapid non-weight coefficient control of a multilevel converter. However, this method cannot be used to achieve predicted torque control of a permanent magnet synchronous motor. The application number is 201810724498.4, the name of the invention is a linear induction motor weightless coefficient model prediction thrust control method, and the linear induction motor weightless coefficient model prediction thrust control method taking thrust and conjugate thrust as objective functions is provided. However, the invention is directed to a linear induction motor and cannot be directly applied to a permanent magnet synchronous motor. The invention discloses a non-weight prediction torque control method of a permanent magnet motor system based on discrete duty ratio, which is 201811426304.9 and is characterized in that a reference voltage command is calculated according to a flux linkage and a torque command, and an objective function is established by using the reference voltage, so that a weight coefficient is eliminated. Although the method can realize the non-weight prediction torque control of the permanent magnet motor system, the complicated calculation is required to obtain the target voltage vector. The document [ Xuyanping, Liyuyuan, Zhoujin. permanent magnet synchronous motor double-model prediction torque control strategy [ J ]. power electronic technology, 2018,52(06):37-39] provides permanent magnet synchronous motor double-model prediction torque control, effectively reduces torque pulsation and reduces the calculation amount of an algorithm. However, this method requires the design of weighting factors, which is difficult to design.
Disclosure of Invention
Aiming at the technical problem of complex calculated amount of the existing control method for predicting the torque of the permanent magnet synchronous motor, the invention provides a permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit, two objective functions of torque and flux linkage are established, and a new objective function without a weight coefficient is obtained through per unit of standard deviation, so that the non-weight coefficient prediction torque control of the permanent magnet synchronous motor is finally realized, the algorithm complexity is simplified, and the torque ripple is reduced.
The technical scheme of the invention is realized as follows:
a permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit comprises the following steps:
step one, sampling three-phase current i of permanent magnet synchronous motor at moment ka、ib、icAnd angular velocity omega of permanent magnet synchronous motorrRotor position angle thetarAnd apply three-phase current ia、ib、icObtaining the output current i under a static αβ coordinate system through coordinate transformationαAnd iβ;
Step two, outputting the current iα、iβUsing rotor position angle thetarDetermination of d-axis current i by coordinate transformationdAnd q-axis current iqThen the angular velocity omega in the step one is usedrRotor position angle thetarAnd d-axis current idSubstituting into the equivalent back electromotive force mathematical model of the permanent magnet synchronous motor to calculate the equivalent back electromotive force e of the permanent magnet synchronous motorαAnd eβ;
Step three, according to the switching state S of the voltage source invertera、Sb、ScObtaining the voltage vector V output by the voltage source inverteri(SaSbSc) Where i is 0,1,2,3,4,5,6,7, switch state Sa、Sb、ScIs equal to 0 or 1;
step four, according to the direct current voltage U of the voltage source inverterdcCalculating the voltage vector V in step threei(SaSbSc) The output voltage of the corresponding inverter, i.e. the stator voltage u of the PMSMαiAnd uβi;
Step five, the output current i obtained in the step oneα、iβAnd the equivalent back electromotive force e obtained in the step twoα、eβAnd the stator voltage u obtained in the step fourαi、uβiPredicting output current i at k +1 moment by substituting discrete current prediction model of permanent magnet synchronous motorαi(k +1) and iβi(k+1);
Step six, the output current i obtained in the step five is usedαi(k+1)、iβi(k +1) and the equivalent back electromotive force e obtained in the second stepα、eβSubstituting the stator flux linkage prediction model of the permanent magnet synchronous motor to predict the stator flux linkage psi at the k +1 momentαi(k+1)、ψβi(k +1) and calculating the stator flux linkage amplitude psisi(k+1);
Step seven, the output current i obtained in the step five is usedαi(k+1)、iβi(k +1) and the stator flux linkage psi obtained in step sixαi(k+1)、ψβiSubstituting (k +1) into the torque prediction model of the permanent magnet synchronous motor to calculate the torque T of the permanent magnet synchronous motorei(k+1);
Step eight, setting reference motor torque according to load requirementsAnd according to a reference motor torqueCalculation referenceStator flux linkage
Step nine, calculating reference motor torqueTorque T obtained by subtracting step seveneiThe absolute value of (k +1) yields the torque target function gTiCalculating the reference stator flux linkageSubtracting the stator flux linkage amplitude psi obtained in the step sixsiObtaining a stator flux linkage objective function g by the absolute value of (k +1)ψiThen, for eight torque target functions gTiAnd eight stator flux linkage objective functions gψiAveraging to obtain the average value of the torque target functionAnd stator flux linkage objective function mean
Step ten, obtaining a torque target function g according to the step nineTiAnd stator flux linkage objective function gψiAnd torque target function meanAnd stator flux linkage objective function meanCalculating to obtain a standard deviation g of a torque target functionTσAnd standard deviation g of flux linkage objective functionψσ;
Step eleven, obtaining a torque target function g according to the step nineTiStator flux linkage objective function gψiMean value of torque objective functionMean value of stator flux linkage objective functionAnd the standard deviation g of the torque target function obtained in the step tenTσStandard deviation g of magnetic linkage objective functionψσSubstituting the mathematical model of the target function based on per unit to calculate a new target function Gi;
Twelfth, comparing the target function GiValue of (d), the smallest objective function GiCorresponding voltage vector Vi(SaSbSc) And the optimal vector is used for controlling the permanent magnet synchronous motor.
Preferably, the three-phase current i in the first stepa、ib、icObtaining the output current i under a static αβ coordinate system through coordinate transformationαAnd iβComprises the following steps:
preferably, the d-axis and q-axis currents i in the second stepdAnd current iqThe obtaining method comprises the following steps:
the equivalent back electromotive force mathematical model of the permanent magnet synchronous motor is as follows:
wherein, ω isrAngular velocity, psi, of a permanent magnet synchronous machinefIs a permanent magnet flux linkage, L, of a permanent magnet synchronous motordAnd LqAre all stator inductances of the permanent magnet synchronous motor.
Preferably, the voltage vector V in step threei(SaSbSc) The obtaining method comprises the following steps:
Satable 1 (the attached drawings)The upper pipe of the a-phase bridge arm of the bidirectional AC-DC converter is connected, and the lower pipe is disconnected;
Sawhen the upper tube of the a-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube of the a-phase bridge arm of the bidirectional alternating current-direct current converter is turned on, the upper tube of the a-phase bridge arm of the bidirectional alternating current-direct current converter is;
Sbwhen the upper tube of the b-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube of the b-phase bridge arm of the bidirectional alternating current-direct current converter is turned on, the upper tube of the b-phase bridge arm of the bidirectional alternating current-direct current converter is;
Scwhen the upper tube of the c-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube is turned on, and when the upper tube of the c-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube of the c-phase bridge arm of;
if Sa=0,Sb=0,ScVoltage vector is denoted as V when it is 00(000);
If Sa=1,Sb=0,ScVoltage vector is denoted as V when it is 01(100);
If Sa=1,Sb=1,ScVoltage vector is denoted as V when it is 02(110);
If Sa=0,Sb=1,ScVoltage vector is denoted as V when it is 03(010);
If Sa=0,Sb=1,ScVoltage vector is denoted as V ═ 14(011);
If Sa=0,Sb=0,ScVoltage vector is denoted as V ═ 15(001);
If Sa=1,Sb=0,ScVoltage vector is denoted as V ═ 16(101);
If Sa=1,Sb=1,ScVoltage vector is denoted as V ═ 17(111)。
Preferably, the stator voltage u of the permanent magnet synchronous motor in step fourαiAnd uβiThe obtaining method comprises the following steps:
wherein i is 0,1,2,3,4,5,6,7, SaiIs a voltage vector Vi(SaSbSc) Corresponding switch state Sa,SbiIs a voltage vector Vi(SaSbSc) Corresponding switch state Sb,SciIs a voltage vector Vi(SaSbSc) Corresponding switch state Sc。
Preferably, the discrete current prediction model of the permanent magnet synchronous motor in the step five is as follows:
wherein e isαAnd eβAre all equivalent back electromotive forces, T, of PMSMsFor a sampling period, RsIs the stator resistance, L, of a permanent magnet synchronous machineqIs the stator inductance of the permanent magnet synchronous motor;
the stator flux linkage prediction model of the permanent magnet synchronous motor in the sixth step is as follows:
the stator flux linkage amplitude psi in the sixth stepsiThe method for obtaining (k +1) is as follows:
the torque prediction model of the permanent magnet synchronous motor in the step seven comprises the following steps:
Preferably, said step eight uses a reference motor torque Te *Calculating the reference flux linkage psis *The method comprises the following steps:
Preferably, the torque target function g in the step nineTiAnd stator flux linkage objective function gψiThe obtaining method comprises the following steps:
the torque target function mean valueAnd stator flux linkage objective function meanThe obtaining method comprises the following steps:
preferably, the standard deviation g of the torque target function in the step ten isTσAnd standard deviation g of flux linkage objective functionψσThe obtaining method comprises the following steps:
preferably, the mathematical model of the per-unit-based objective function in the step eleven is as follows:
the beneficial effect that this technical scheme can produce: by adopting two objective functions of torque and flux linkage, the control of the torque and flux linkage of the permanent magnet synchronous motor can be realized without a weight coefficient, the complexity of system design and debugging is reduced, torque ripples are reduced, and the torque control precision is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a block diagram of the overall structure of the present invention.
FIG. 2 is a graph of simulation results of a document [ Xuyanping, Liyuan, Zhoujin ] permanent magnet synchronous motor double model predicted torque control strategy [ J ]. Power electronics, 2018,52(06):37-39 ]; (a) the simulation result diagram of the torque error and the torque is shown, and the simulation result diagram of the flux linkage error and the flux linkage is shown.
FIG. 3 is a graph of simulation results for the present invention; (a) the simulation result diagram of the torque error and the torque is shown, and the simulation result diagram of the flux linkage error and the flux linkage is shown.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive effort based on the embodiments of the present invention, are within the scope of the present invention.
As shown in fig. 1, a method for predicting torque based on per unit without weighting factor of a permanent magnet synchronous motor includes the following steps:
step one, sampling three-phase current i of permanent magnet synchronous motor at moment ka、ib、icAnd angular velocity omega of permanent magnet synchronous motorrRotor position angle thetarAnd the three-phase current i is converted by using the formula (1)a、ib、icObtaining the output current i under a static αβ coordinate system through coordinate transformationαAnd iβ:
Step two, outputting the current i according to a formula (2)α、iβAnd rotor position angle thetarSeparately determining d-axis and q-axis currents i by coordinate transformationdAnd current iq:
Then the angular velocity omega is adjustedrRotor position angle thetarAnd d-axis current idSubstituting into the equivalent back electromotive force mathematical model of the permanent magnet synchronous motor to calculate the equivalent back electromotive force e of the permanent magnet synchronous motorαAnd eβAs shown in equation (3):
wherein psifIs a permanent magnet flux linkage, L, of a permanent magnet synchronous motordAnd LqAre all stator inductances of the permanent magnet synchronous motor.
Step three, according to the switching state S of the voltage source invertera、Sb、ScObtaining the voltage vector V output by the voltage source inverteri(SaSbSc) Where i is 0,1,2,3,4,5,6,7, switch state Sa、Sb、ScIs equal to 0 or 1:
Sawhen the upper tube of the a-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube of the a-phase bridge arm of the bidirectional alternating current-direct current converter is turned on, the upper tube of the a-phase bridge arm of the bidirectional alternating current-direct current converter is;
Sbwhen the upper tube of the b-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube of the b-phase bridge arm of the bidirectional alternating current-direct current converter is turned on, the upper tube of the b-phase bridge arm of the bidirectional alternating current-direct current converter is;
Scwhen the upper tube of the c-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube is turned on, and when the upper tube of the c-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube of the c-phase bridge arm of;
if Sa=0,Sb=0,ScVoltage vector is denoted as V when it is 00(000);
If Sa=1,Sb=0,ScVoltage vector is denoted as V when it is 01(100);
If Sa=1,Sb=1,ScVoltage vector is denoted as V when it is 02(110);
If Sa=0,Sb=1,ScVoltage vector is denoted as V when it is 03(010);
If Sa=0,Sb=1,ScVoltage vector is denoted as V ═ 14(011);
If Sa=0,Sb=0,ScVoltage vector is denoted as V ═ 15(001);
If Sa=1,Sb=0,ScVoltage vector is denoted as V ═ 16(101);
If Sa=1,Sb=1,ScVoltage vector is denoted as V ═ 17(111)。
Therefore, eight voltage vectors output by the voltage source inverter are respectively marked as V0(000)、V1(100)、V2(110)、V3(010)、V4(011)、V5(001)、V6(101) And V7(111)。
Step four, according to the direct current voltage U of the voltage source inverterdcCalculating the voltage vector V in step threei(SaSbSc) The output voltage of the corresponding inverter, i.e. the stator voltage u of the PMSMαiAnd uβiAs shown in equation (4):
wherein i is 0,1,2,3,4,5,6,7, SaiIs equal to the voltage vector Vi(SaSbSc) Corresponding switch state Sa,SbiIs equal to the voltage vector Vi(SaSbSc) Corresponding switch state Sb,SciIs equal to the voltage vector Vi(SaSbSc) Corresponding switch state Sc。
Step five, the output current i obtained in the step oneα、iβAnd the equivalent back electromotive force e obtained in the step twoα、eβAnd the stator voltage u obtained in the step fourαi、uβiPredicting output current i at k +1 moment by substituting discrete current prediction model of permanent magnet synchronous motorαi(k +1) and iβi(k +1) as shown in formula (5):
wherein i is 0,1,2,3,4,5,6,7, TsFor a sampling period, RsIs the stator resistance, L, of a permanent magnet synchronous machineqIs the stator inductance of the permanent magnet synchronous motor.
Step six, the output current i obtained in the step five is usedαi(k+1)、iβi(k +1) and the equivalent back electromotive force e obtained in the second stepα、eβSubstituting the stator flux linkage prediction model of the permanent magnet synchronous motor to predict the stator flux linkage psi at the k +1 momentαi(k+1)、ψβi(k +1) and calculating the stator flux linkage amplitude psisi(k +1) as shown in formula (6) and formula (7):
wherein i is 0,1,2,3,4,5,6,7, eαAnd eβAre all equivalent back electromotive forces, omega, of permanent magnet synchronous motorsrIs the angular velocity of the permanent magnet synchronous motor.
Step seven, obtaining the output current i according to the step fiveαi(k+1)、iβi(k +1), stator flux linkage psi obtained in step sixαi(k+1)、ψβi(k +1) and a torque prediction model of the permanent magnet synchronous motor to calculate the torque T of the permanent magnet synchronous motorei(k +1) as shown in formula (8):
wherein i is 0,1,2,3,4,5,6,7, npThe number of pole pairs of the permanent magnet synchronous motor is shown.
Step eight, setting reference motor torque according to load requirementsAnd according to a reference motor torqueObtaining the reference stator flux linkage through the calculation of the formula (9)
Wherein n ispIs the pole pair number psi of the permanent magnet synchronous motorfIs a permanent magnet flux linkage, L, of a permanent magnet synchronous motorqIs the stator inductance of the permanent magnet synchronous motor.
Step nine, calculating a reference stator flux linkageSubtracting the stator flux linkage amplitude psi obtained in the step sixsiObtaining a stator flux linkage objective function g by the absolute value of (k +1)ψiCalculating a reference motor torqueTorque T obtained by subtracting step seveneiThe absolute value of (k +1) yields the torque target function gTiAs shown in equation (10):
then, the eight torque target functions g are respectively processed according to the formula (11)TiAnd eight stator flux linkage objective functions gψiAveraging to obtain the average value of the torque target functionAnd stator flux linkage objective function mean
Step ten, obtaining a torque target function g according to the step nineTiAnd stator flux linkage objective function gψiAnd torque target function meanAnd stator flux linkage objective function meanCalculating to obtain a standard deviation g of a torque target functionTσAnd standard deviation g of flux linkage objective functionψσAs shown in equation (12):
wherein i is 0,1,2,3,4,5,6, 7.
Step eleven, obtaining a torque target function g according to the step nineTiStator flux linkage objective function gψiTorque target function flatMean valueMean value of stator flux linkage objective functionStep ten, obtaining a standard deviation g of a torque target functionTσStandard deviation g of magnetic linkage objective functionψσAnd calculating to obtain a new target function G based on the mathematical model of the per-unit target functioniAs shown in equation (13):
wherein i is 0,1,2,3,4,5,6, 7.
Twelfth, comparing the target function GiValue of (d), the smallest objective function GiCorresponding voltage vector Vi(SaSbSc) The most vector is used for controlling the permanent magnet synchronous motor.
In order to verify the effectiveness of the invention, simulation verification is carried out on the scheme of the control of the predicted torque of the traditional permanent magnet synchronous motor. During simulation, the direct-current side voltage U of the grid-connected inverterdcIs 600V, and the stator resistance R of the motors0.0154 omega, and a permanent magnet flux linkage psifIs 1.5Wb, d-axis inductance LdIs 4mH, q-axis inductance Lq9mH, the number of pole pairs n of the motorpTo 3, reference is made to the motor torqueIs 300 Nm. FIG. 2 shows a bimodal predicted torque control strategy for PMSM [ Xuyanping, Liyuan, Zhoujin ]]Power electronics, 2018,52(06):37-39]Fig. 3 shows the simulation result of the present invention. As shown in fig. 2, in the conventional permanent magnet synchronous motor predicted torque control scheme, because there is no mature weight factor design theory at present, it is difficult to obtain optimal stator flux linkage and torque control effects, that is, torque error and torque waveform, and flux linkage error and flux linkage waveform wave simultaneouslyThe dynamic phase is relatively large; as shown in FIG. 3, two objective functions of torque and flux linkage are established, and a new objective function without weight coefficients is obtained through standard deviation per unit, so that the torque and flux linkage control of the permanent magnet synchronous motor is realized, and the torque error, the torque waveform, the flux linkage error and the flux linkage waveform fluctuation range are reduced. The invention does not need a weight coefficient, thereby not only reducing the complexity of system design and debugging, but also being beneficial to reducing torque ripple and improving the torque control precision.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (10)
1. A permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit is characterized by comprising the following steps:
step one, sampling three-phase current i of permanent magnet synchronous motor at moment ka、ib、icAnd angular velocity omega of permanent magnet synchronous motorrRotor position angle thetarAnd apply three-phase current ia、ib、icObtaining the output current i under a static αβ coordinate system through coordinate transformationαAnd iβ;
Step two, outputting the current iα、iβUsing rotor position angle thetarDetermination of d-axis current i by coordinate transformationdAnd q-axis current iqThen the angular velocity omega in the step one is usedrRotor position angle thetarAnd d-axis current idSubstituting into the equivalent back electromotive force mathematical model of the permanent magnet synchronous motor to calculate the equivalent back electromotive force e of the permanent magnet synchronous motorαAnd eβ;
Step three, according to the switching state S of the voltage source invertera、Sb、ScObtaining the voltage vector V output by the voltage source inverteri(SaSbSc) Wherein i is 0,1,2,3,4,5,6,7, KOff state Sa、Sb、ScIs equal to 0 or 1;
step four, according to the direct current voltage U of the voltage source inverterdcCalculating the voltage vector V in step threei(SaSbSc) The output voltage of the corresponding inverter, i.e. the stator voltage u of the PMSMαiAnd uβi;
Step five, the output current i obtained in the step oneα、iβAnd the equivalent back electromotive force e obtained in the step twoα、eβAnd the stator voltage u obtained in the step fourαi、uβiPredicting output current i at k +1 moment by substituting discrete current prediction model of permanent magnet synchronous motorαi(k +1) and iβi(k+1);
Step six, the output current i obtained in the step five is usedαi(k+1)、iβi(k +1) and the equivalent back electromotive force e obtained in the second stepα、eβSubstituting the stator flux linkage prediction model of the permanent magnet synchronous motor to predict the stator flux linkage psi at the k +1 momentαi(k+1)、ψβi(k +1) and calculating the stator flux linkage amplitude psisi(k+1);
Step seven, the output current i obtained in the step five is usedαi(k+1)、iβi(k +1) and the stator flux linkage psi obtained in step sixαi(k+1)、ψβiSubstituting (k +1) into the torque prediction model of the permanent magnet synchronous motor to calculate the torque T of the permanent magnet synchronous motorei(k+1);
Step eight, setting reference motor torque according to load requirementsAnd according to a reference motor torqueCalculating reference stator flux linkage
Step nine, calculating a reference motorTorque momentTorque T obtained by subtracting step seveneiThe absolute value of (k +1) yields the torque target function gTiCalculating the reference stator flux linkageSubtracting the stator flux linkage amplitude psi obtained in the step sixsiObtaining a stator flux linkage objective function g by the absolute value of (k +1)ψiThen, for eight torque target functions gTiAnd eight stator flux linkage objective functions gψiAveraging to obtain the average value of the torque target functionAnd stator flux linkage objective function mean
Step ten, obtaining a torque target function g according to the step nineTiAnd stator flux linkage objective function gψiAnd torque target function meanAnd stator flux linkage objective function meanCalculating to obtain a standard deviation g of a torque target functionTσAnd standard deviation g of flux linkage objective functionψσ;
Step eleven, obtaining a torque target function g according to the step nineTiStator flux linkage objective function gψiMean value of torque objective functionMean value of stator flux linkage objective functionAnd the standard deviation g of the torque target function obtained in the step tenTσStandard deviation g of magnetic linkage objective functionψσSubstituting the mathematical model of the target function based on per unit to calculate a new target function Gi;
Twelfth, comparing the target function GiValue of (d), the smallest objective function GiCorresponding voltage vector Vi(SaSbSc) And the optimal vector is used for controlling the permanent magnet synchronous motor.
2. The per-unit-based permanent magnet synchronous motor non-weight coefficient prediction torque control method according to claim 1, wherein the three-phase current i in the first stepa、ib、icObtaining the output current i under a static αβ coordinate system through coordinate transformationαAnd iβComprises the following steps:
3. the method according to claim 1, wherein the d-axis and q-axis currents i in the second step are currents idAnd current iqThe obtaining method comprises the following steps:
the equivalent back electromotive force mathematical model of the permanent magnet synchronous motor is as follows:
4. The method according to claim 1, wherein the voltage vector V in the third step is a voltage vector Vi(SaSbSc) The obtaining method comprises the following steps:
Sa1 represents that an upper tube of a phase bridge arm of a bidirectional alternating current-direct current converter is conducted and a lower tube is turned off;
Sawhen the upper tube of the a-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube of the a-phase bridge arm of the bidirectional alternating current-direct current converter is turned on, the upper tube of the a-phase bridge arm of the bidirectional alternating current-direct current converter is;
Sb1 represents that the upper tube of a b-phase bridge arm of the bidirectional AC-DC converter is conducted and the lower tube is turned off;
Sbwhen the upper tube of the b-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube of the b-phase bridge arm of the bidirectional alternating current-direct current converter is turned on, the upper tube of the b-phase bridge arm of the bidirectional alternating current-direct current converter is;
Sc1 represents that the upper tube of the c-phase bridge arm of the bidirectional AC-DC converter is conducted and the lower tube is turned off;
Scwhen the upper tube of the c-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube is turned on, and when the upper tube of the c-phase bridge arm of the bidirectional alternating current-direct current converter is turned off, the lower tube of the c-phase bridge arm of;
if Sa=0,Sb=0,ScVoltage vector is denoted as V when it is 00(000);
If Sa=1,Sb=0,ScVoltage vector is denoted as V when it is 01(100);
If Sa=1,Sb=1,ScVoltage vector is denoted as V when it is 02(110);
If Sa=0,Sb=1,ScVoltage vector is denoted as V when it is 03(010);
If Sa=0,Sb=1,ScVoltage vector is denoted as V ═ 14(011);
If Sa=0,Sb=0,ScVoltage vector is denoted as V ═ 15(001);
If Sa=1,Sb=0,ScVoltage vector is denoted as V ═ 16(101);
If Sa=1,Sb=1,ScVoltage vector is denoted as V ═ 17(111)。
5. The method according to claim 1 or 4, wherein the stator voltage u of the PMSM in step four is the stator voltage u of the PMSMαiAnd uβiThe obtaining method comprises the following steps:
6. The method for controlling the torque without the weight coefficient prediction of the per unit permanent magnet synchronous motor according to claim 5, wherein the discrete current prediction model of the permanent magnet synchronous motor in the step five is:
wherein e isαAnd eβAre all equivalent back electromotive forces, T, of PMSMsFor a sampling period, RsIs the stator resistance, L, of a permanent magnet synchronous machineqIs the stator inductance of the permanent magnet synchronous motor;
the stator flux linkage prediction model of the permanent magnet synchronous motor in the sixth step is as follows:
the stator flux linkage amplitude psi in the sixth stepsiThe method for obtaining (k +1) is as follows:
the torque prediction model of the permanent magnet synchronous motor in the step seven comprises the following steps:
7. The PMSM non-weight-coefficient predicted torque control method according to claim 6, wherein reference motor torque is used in step eightCalculating reference flux linkageThe method comprises the following steps:
8. The method according to claim 7, wherein the torque objective function g in the step nine is a torque objective function gTiAnd stator flux linkage objective function gψiThe obtaining method comprises the following steps:
the torque target function mean valueAnd stator flux linkage objective function meanThe obtaining method comprises the following steps:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910399476.XA CN110112979B (en) | 2019-05-14 | 2019-05-14 | Permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910399476.XA CN110112979B (en) | 2019-05-14 | 2019-05-14 | Permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110112979A CN110112979A (en) | 2019-08-09 |
CN110112979B true CN110112979B (en) | 2020-09-04 |
Family
ID=67489869
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910399476.XA Expired - Fee Related CN110112979B (en) | 2019-05-14 | 2019-05-14 | Permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110112979B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111224604A (en) * | 2019-11-21 | 2020-06-02 | 西安理工大学 | Simplified model prediction control method for asynchronous motor |
CN110995072B (en) * | 2019-12-19 | 2021-09-10 | 华中科技大学 | Motor rotor position estimation method |
CN112688587B (en) * | 2020-12-28 | 2022-02-15 | 珠海创芯科技有限公司 | Robust prediction control method of impedance source inverter |
CN112821832B (en) * | 2021-02-10 | 2022-09-30 | 北方工业大学 | Control method and device of permanent magnet synchronous motor and motor controller |
CN112865654B (en) * | 2021-04-13 | 2022-09-09 | 山东大学 | Torque maximum utilization control system and method for permanent magnet magnetic concentration type synchronous reluctance motor |
CN113411028B (en) * | 2021-06-28 | 2022-12-30 | 珠海创芯科技有限公司 | Permanent magnet synchronous motor model prediction control method and device, electronic equipment and medium |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3457557B1 (en) * | 2017-09-15 | 2021-03-17 | Hamilton Sundstrand Corporation | Active stabilization of dc link in motor drive systems |
CN108736787B (en) * | 2018-05-30 | 2019-11-05 | 中车青岛四方车辆研究所有限公司 | Permanent-magnet synchronous traction drive control system and method |
-
2019
- 2019-05-14 CN CN201910399476.XA patent/CN110112979B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN110112979A (en) | 2019-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110112979B (en) | Permanent magnet synchronous motor non-weight coefficient prediction torque control method based on per unit | |
CN110120763B (en) | Permanent magnet synchronous motor non-weight coefficient prediction torque control method | |
Pan et al. | Voltage balancing control of diode-clamped multilevel rectifier/inverter systems | |
CN108631672B (en) | Permanent magnet synchronous motor prediction flux linkage control method considering optimal duty ratio modulation | |
WO2022110277A1 (en) | Model predictive decomposition control method and apparatus for open-winding five-phase permanent magnet synchronous motor | |
CN110995076A (en) | Permanent magnet synchronous motor model prediction current control method | |
CN109412482B (en) | Unified predictive control method for quasi-Z-source inverter-permanent magnet synchronous motor system | |
CN107196571B (en) | Double-motor series prediction type direct torque control method | |
CN110829922B (en) | Semi-controlled open winding PMSG dual-vector model prediction flux linkage control method | |
CN108736778B (en) | Dual-vector prediction flux linkage control method for permanent magnet synchronous motor | |
CN111726046B (en) | Asymmetric six-phase PMSM model prediction flux linkage control method considering duty ratio optimization | |
JP2010068653A (en) | Inverter control apparatus and motor drive system | |
CN110504889B (en) | Fault-tolerant direct torque control method for five-phase permanent magnet synchronous motor | |
CN110011588B (en) | Semi-control open winding permanent magnet synchronous generator model prediction flux linkage control method | |
CN109600065B (en) | Discontinuous pulse width modulation method for three-level converter | |
Guo et al. | Hybrid multi-vector modulated model predictive control strategy for voltage source inverters based on a new visualization analysis method | |
CN111800050A (en) | Permanent magnet synchronous motor three-vector model predicted torque control method based on voltage vector screening and optimization | |
Shen et al. | A mixed SVPWM technique for three-phase current reconstruction with single DC negative rail current sensor | |
Lee et al. | Speed-sensorless DTC-SVM for matrix converter drives with simple nonlinearity compensation | |
CN109347392B (en) | Instantaneous power decoupling control method for open-winding permanent magnet synchronous motor | |
CN105871293B (en) | Low-cost double-PWM power converter model prediction control method | |
CN116667732A (en) | Model prediction current control method for three-level inverter permanent magnet synchronous motor | |
Vargas et al. | Predictive strategy to reduce common-mode voltages on power converters | |
Reimers et al. | Switched reluctance motor drive with three-phase integrated battery charger for electric vehicle applications | |
Siami et al. | Improved predictive torque control of a permanent magnet synchronous motor fed by a matrix converter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200904 |
|
CF01 | Termination of patent right due to non-payment of annual fee |