CN110784145B - Permanent magnet synchronous motor full-speed domain model prediction flux linkage control method - Google Patents

Permanent magnet synchronous motor full-speed domain model prediction flux linkage control method Download PDF

Info

Publication number
CN110784145B
CN110784145B CN201911124614.XA CN201911124614A CN110784145B CN 110784145 B CN110784145 B CN 110784145B CN 201911124614 A CN201911124614 A CN 201911124614A CN 110784145 B CN110784145 B CN 110784145B
Authority
CN
China
Prior art keywords
flux linkage
speed
formula
moment
obtaining
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.)
Active
Application number
CN201911124614.XA
Other languages
Chinese (zh)
Other versions
CN110784145A (en
Inventor
於锋
霍闯
茅靖峰
李凯凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jieyang Shengjia Motor Co.,Ltd.
Shenzhen Lizhuan Technology Transfer Center Co ltd
Original Assignee
Nantong University
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 Nantong University filed Critical Nantong University
Priority to CN201911124614.XA priority Critical patent/CN110784145B/en
Publication of CN110784145A publication Critical patent/CN110784145A/en
Application granted granted Critical
Publication of CN110784145B publication Critical patent/CN110784145B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/141—Flux estimation
    • 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/10—Arrangements for controlling torque ripple, e.g. providing reduced torque ripple
    • 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/34—Modelling or simulation for control purposes

Abstract

The invention discloses a permanent magnet synchronous motorA flux linkage control method for speed domain model prediction includes the steps that firstly, a reference torque T is obtained through a rotating speed outer ring PI controllere ref(ii) a Then obtaining the electrical angle theta of the permanent magnet synchronous motorrAnd electrical angular velocity ωrObtaining three-phase stator current at the moment k, and obtaining d-q axis components of the stator current at the moment k after coordinate transformation; then, calculating flux linkage by combining a flux linkage equation and a load angle to obtain a flux linkage predicted value and a flux linkage reference value at the (k +1) moment; further, a full-speed domain cost function is constructed by using the predicted value of the flux linkage at the (k +1) moment and the reference value, and the optimal voltage vector of the inverter is obtained by minimizing the full-speed domain cost function; and finally, calculating the duty ratio according to the q-axis flux linkage dead-beat concept, and distributing the time of the optimal voltage vector and the zero vector acting on the inverter. The invention can effectively reduce current and torque pulsation, has good dynamic and steady performance, and is applicable to a constant torque area and a constant power area.

Description

Permanent magnet synchronous motor full-speed domain model prediction flux linkage control method
Technical Field
The invention relates to a permanent magnet synchronous motor full-speed domain model prediction flux linkage control method, and belongs to the field of motor driving and control.
Background
A Permanent Magnet Synchronous Motor (PMSM) has the advantages of small loss, low temperature rise, high power factor, high starting torque, short starting time, high overload capability and the like. The conventional control method for the permanent magnet synchronous motor mainly includes Vector Control (VC) and Direct Torque Control (DTC). VC is a linear control method for independently controlling flux linkage and torque by controlling exciting current and torque current. However, coordinate transformation is complex, a current regulator is usually PI control, and the PI control has inertia link and hysteresis, so that the dynamic performance of the system is poor. The DTC adopts a hysteresis controller and an offline switching meter, simplifies the system structure, but has large switching frequency change, and causes poor performance when the motor runs at low speed. As microprocessor performance has improved, Model Predictive Control (MPC) has gained widespread attention. MPC uses mathematical model to predict the change trend of system variable at next time to select the corresponding optimal control action at current time for rolling optimization, thereby accurately controlling the control target. Therefore, a novel model prediction control method derived from the DTC is provided, and the principle of the method is similar to that of the DTC, namely, the torque and the flux linkage are directly predicted through a prediction model of a system, and an optimal voltage vector is selected through a corresponding cost function. The algorithm has the advantages of simple structure, quick dynamic response and the like, but the weight coefficient design process is complex.
Aiming at the problem of complex weight coefficient design, the model prediction flux linkage control converts the control on the stator flux linkage amplitude and the electromagnetic torque into the control on the equivalent stator flux linkage vector by deeply deducing the analytic relation between the flux linkage and the torque, thereby eliminating the complicated weight coefficient design in the traditional method, and the algorithm is simple and is easy to realize. Since the MPC strategy has a good application prospect in the field of PMSM driving, in recent years, many scholars are dedicated to research and improvement of MPC in a full-speed motor domain operation state. When the motor runs below the basic speed, namely in a constant torque area, the Maximum torque current ratio (MTPA) is adopted for control, so that the reluctance torque can be effectively utilized, and the system efficiency of the motor is improved. When the motor operates above the basic speed, namely in a constant power area, Flux Weakening Control (FWC) is adopted, so that the operating speed range of the PMSM can be widened. In addition, the current and torque ripple can be effectively reduced by adopting the double-vector model predictive control of duty ratio control.
The invention patent 201810417636.4, entitled model-based torque prediction control method for a permanent magnet synchronous motor considering switching frequency optimization, discloses a method and device for model-based torque prediction control when the permanent magnet synchronous motor operates in a constant torque region and a constant power region. The method utilizes a prediction model to predict the d-q component i of the stator current at the (k +1) moment on lined k+1，iq k+1And electromagnetic torque Te k+1Constructing a full-speed domain cost function by combining a plurality of control targets; and finally, the optimal voltage vector of the inverter is obtained by minimizing the full-speed domain cost function, so that the switching frequency of the inverter is reduced, and a better system control effect is obtained. However, the patent does not consider the problem of complicated weight coefficient design in the full-speed domain cost function and does not consider the problem in model prediction torque controlThe current and the torque ripple are large.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the prior art, the permanent magnet synchronous motor full-speed domain model prediction flux linkage control method is provided, good dynamic and steady-state performance can be obtained while weight coefficients are eliminated, current and torque pulsation are effectively reduced, and the method is suitable for a constant torque area and a constant power area.
The technical scheme is as follows: a permanent magnet synchronous motor full-speed domain model prediction flux linkage control method comprises the following steps: firstly, a reference torque T is obtained through a rotating speed outer ring PI controllere ref(ii) a Then obtaining the electrical angle theta of the permanent magnet synchronous motor from the motor encoderrAnd electrical angular velocity ωrAnd acquiring three-phase stator current i at the time k by using a current sensora，ibAnd icObtaining d-q axis component i of stator current at the moment k after coordinate transformationdAnd iq(ii) a Then, flux linkage calculation is carried out by combining a flux linkage equation and a load angle to obtain a predicted value psi of flux linkage at the (k +1) momentsd(k+1)、ψsq(k +1) and flux linkage reference value psisd ref(k+1)、ψsqref (k + 1); further, a full-speed domain cost function is constructed by using the flux linkage predicted value and the flux linkage reference value at the (k +1) moment, and an optimal voltage vector of the inverter is obtained by minimizing the full-speed domain cost function; and finally, calculating the duty ratio according to the q-axis flux linkage dead-beat concept, and distributing the time of the optimal voltage vector and the zero vector acting on the inverter.
Further, the reference torque Te refThe acquisition method comprises the following steps: the difference e between the reference speed and the actual speed of the motor is calculatednInputting a rotating speed PI controller, and obtaining a reference torque T according to a formula (1)e ref；
Figure GDA0003170828780000021
In the formula, kpAnd kiProportional gain and integral gain of the rotating speed PI controller are respectively shown, and s represents a complex variable.
Further, the electrical angle θrElectrical angular velocity omegarAnd d-q axis component i of stator current at time kd，iqThe acquisition method comprises the following steps: obtaining the electrical angle theta of the permanent magnet synchronous motor from the encoderrThen, the electrical angle theta is obtained through the formula (2)rWith respect to the differentiation of time, an electrical angular velocity ω is obtainedr(ii) a Measuring k-time three-phase stator current i of permanent magnet synchronous motor by using current sensora，ibAnd icObtaining d-q axis component i of stator current at the moment k after coordinate transformationdAnd iq；
Figure GDA0003170828780000022
Further, a predicted value psi of flux linkage at the time (k +1) is calculatedsd(k+1)、ψsqReference value ψ of flux linkage at time (k +1) and (k +1)sd ref(k+1)、ψsq refThe method of (k +1) is: the obtained d-q axis current component idAnd iqAngular velocity ω of rotorrAnd rotor electrical angle thetarInputting a model prediction flux linkage control module, obtaining a prediction current model at the moment (k +1) according to a formula (3), and then obtaining a predicted flux linkage value psi at the moment (k +1) according to a formula (4)sd(k+1)、ψsq(k + 1); obtaining the load angle delta and the electromagnetic torque T according to the formula (5)eAnd the load angle delta is derived according to the formula (6) to obtain the load angle increment delta of the formula (7), and the reference value delta of the load angle delta at the moment of (k +1) is obtained according to the formula (8)refThe reference value ψ of the flux linkage at the time (k +1) is obtained according to the formula (9)sd ref(k+1)、ψsq ref(k+1)；
Figure GDA0003170828780000031
Figure GDA0003170828780000032
Figure GDA0003170828780000033
Figure GDA0003170828780000034
Figure GDA0003170828780000035
δref＝Δδ+δ (8)
Figure GDA0003170828780000036
In the formula id k+1、iq k+1The predicted value of the current at the moment (k + 1); rsIs a stator phase resistance; l isd、LqThe inductor is a direct axis inductor and a quadrature axis inductor; t is the sampling period of the system; u. ofd k、uq kThe voltage of the stator voltage on the d-q axis component at the moment k; psifIs a rotor permanent magnet flux linkage; n ispIs the number of pole pairs; psisFor the resultant flux linkage psi on the d-q axis component at time ks(k) The amplitude of (d); psisd ref(k+1)、ψsq ref(k +1) is a flux linkage reference value of the stator flux linkage on the d-q axis component at the moment (k + 1); delta TeIs the electromagnetic torque increment; deltarefIs the reference value of the load angle at the moment (k + 1).
Further, constructing a full-speed domain cost function comprises constructing a cost function of the low-medium speed region control target and a cost function of the high-speed region control target in a full-speed domain cost function module;
the method for constructing the low and medium speed region control target comprises the following steps: obtaining a flux linkage error function g of the low and medium speed region at the time of (k +1) according to the formula (10)F(ii) a When the permanent magnet synchronous motor operates in the maximum torque current ratio mode, the low-medium speed region at the moment (k +1) is obtained according to the formula (11)Regional convergence function gMTPA(ii) a Obtaining a flux linkage limiting condition function g of the low and medium speed region at the moment (k +1) according to the formula (12)FMAX(ii) a Obtaining the direction selection function g of the low and medium speed regions at the moment (k +1) according to the formula (13)dir；
Figure GDA0003170828780000041
Figure GDA0003170828780000042
Figure GDA0003170828780000043
Figure GDA0003170828780000044
The method for constructing the high-speed area control target comprises the following steps: calculating a flux linkage error function g 'of the high-speed region at the moment (k + 1)'FFlux linkage error function g 'of high speed region'FFlux linkage error function g of low and medium speed regionFThe consistency is achieved; neglecting the stator resistance voltage drop when the permanent magnet synchronous motor operates stably at high speed to obtain a formula (14), and obtaining a high-speed area convergence function g at the moment of (k +1) according to the formula (15)FW(ii) a Calculating a high-speed region flux linkage limitation condition function g 'at the moment of (k + 1)'FMAXHigh speed region flux linkage constraint function g'FMAXFlux linkage limiting condition function g of low and medium speed regionsFMAXThe consistency is achieved; when the permanent magnet synchronous motor operates above the basic speed, the voltage limiting condition function g of the high-speed area at the moment of (k +1) is obtained according to the formula (16) under the constraint of the maximum output voltage of the inverterumax(ii) a Obtaining a high-speed region stable operation function g at the moment (k +1) according to the formula (17)stab；
Figure GDA0003170828780000045
Figure GDA0003170828780000046
Figure GDA0003170828780000047
Figure GDA0003170828780000051
In the formula usIs the stator voltage; u. ofsmaxIs the maximum output voltage of the inverter; vdcIs a dc bus voltage; lambda [ alpha ]mIs a voltage coefficient; eta is a voltage limiting condition intermediate variable; zeta is the intermediate variable of the motor high speed stable running condition.
Further, the step of obtaining the optimal voltage vector of the inverter by minimizing the full-speed domain cost function comprises the following steps: obtaining a full-speed domain cost function according to the formula (18), respectively substituting eight basic voltage vectors into the cost function, and outputting a switch state S which enables the cost function to be minimumabcFeeding the inverter;
g(min)＝gF+gc+gL (18)
definition of ωcFor the corresponding electrical angular velocity when the permanent magnet synchronous motor is operated at the base speed, when omegar<ωcWhen g isc＝gMTPAAnd g isL＝gFMAX+gdir(ii) a When ω isr>ωcWhen g isc＝gFWAnd g isL＝gFMAX+gumax+gstab。
Further, the duty ratio calculation method comprises: the q-axis flux linkage obtained according to the formula (19) reaches the given value psi at the moment of (k +1) under the combined action of the optimal voltage vector and the zero voltage vectorsq ref(ii) a Combining a stator flux linkage equation and a voltage equation, and obtaining the slope S of the q-axis flux linkage when the zero vector acts according to the formula (20)0The slope S of the q-axis flux linkage at the time of the optimum voltage vector action is obtained from the equation (21)optObtaining the optimum voltage vector action time t according to the equation (22)opt；
Figure GDA0003170828780000052
Figure GDA0003170828780000053
Figure GDA0003170828780000054
Figure GDA0003170828780000055
In the formula, #sqRepresenting flux linkage, ψ, in the q-axis componentsq(k) Is the flux linkage on the q-axis component at time k; s0Is the slope of the q-axis flux linkage when zero vector is applied; soptIs the slope of the q-axis flux linkage when the optimal voltage vector acts; t is toptIs the optimal voltage vector action time;
Figure GDA0003170828780000056
representing the q-axis voltage under the optimal vector at time k.
Has the advantages that: compared with the prior art, the method is based on the model prediction flux linkage control principle, double-vector model prediction control based on duty ratio control is introduced into the MPFC, a full-speed domain cost function of a plurality of control targets including flux linkage control, MTPA optimization, flux linkage limitation, voltage limitation and the like is constructed, the optimal voltage vector acting on the inverter is obtained through the cost function, and after zero vector action is added, not only can good dynamic and steady-state performance be obtained, but also current and torque ripple are effectively reduced, and the method is suitable for a constant torque area and a constant power area.
Drawings
FIG. 1 is a schematic diagram of a permanent magnet synchronous motor full-speed domain model prediction flux linkage control method provided by the invention;
FIG. 2 is a flowchart of a permanent magnet synchronous motor full-speed domain model predictive flux linkage control method provided by the invention;
FIG. 3 is a simulation result of the permanent magnet synchronous motor full-speed domain model predicting the working condition below the flux linkage control base speed; fig. 3 (a) is a simulation result of a steady state in which a single vector model predicts a flux linkage control base speed or lower, and fig. 3 (b) is a simulation result of a steady state in which a double vector model predicts a flux linkage control base speed or lower;
FIG. 4 is a simulation result of the condition below the basic speed of the permanent magnet synchronous motor full-speed domain model prediction flux linkage control method; fig. 4 (a) is a simulation result of a sudden change in rotation speed when the single vector model predicts that the flux linkage control base speed is lower, and fig. 4 (b) is a simulation result of a sudden change in rotation speed when the double vector model predicts that the flux linkage control base speed is lower;
FIG. 5 shows a simulation result of the full-speed domain condition of model predictive flux linkage control.
Detailed Description
The invention is further explained below with reference to the drawings.
A schematic diagram of a full-speed domain model prediction flux linkage control method of a permanent magnet synchronous motor is shown in fig. 1, and the full-speed domain model prediction flux linkage control method comprises a rotating speed PI controller 1, a model prediction flux linkage control module 2, a full-speed domain cost function module 3, a duty ratio calculation module 4, an inverter 5, a coordinate transformation module 6, a permanent magnet synchronous motor 7 and an encoder 8.
Firstly, a reference torque T is obtained through a rotating speed outer ring PI controllere ref(ii) a Then obtaining the electrical angle theta of the permanent magnet synchronous motor from the motor encoderrAnd electrical angular velocity ωrAnd acquiring three-phase stator current i at the time k by using a current sensora，ibAnd icObtaining d-q axis component i of stator current at the moment k after coordinate transformationdAnd iq(ii) a Then, flux linkage calculation is carried out by combining a flux linkage equation and a load angle to obtain a predicted value psi of flux linkage at the (k +1) momentsd(k+1)、ψsq(k +1) and flux linkage reference value psisd ref(k+1)、ψsq ref(k + 1); and further, constructing a full-speed domain cost function by using the flux linkage predicted value and the flux linkage reference value at the (k +1) moment, and performing full-speed domain cost functionObtaining an optimal voltage vector of the inverter by minimizing a full-speed domain cost function; and finally, calculating the duty ratio according to the q-axis flux linkage dead-beat concept, and distributing the time of the optimal voltage vector and the zero vector acting on the inverter.
The method specifically comprises the following steps:
(1) calculating a given torque Te ref: the difference e between the reference speed and the actual speed of the motor is calculatednInputting the rotation speed PI controller 1, and obtaining the reference torque T according to the formula (1)e ref；
Figure GDA0003170828780000061
In the formula, kpAnd kiProportional gain and integral gain of the rotating speed PI controller are respectively shown, and s represents a complex variable.
(2) Calculating the electrical angle thetarElectrical angular velocity omegarAnd d-q component i of stator current at time kdAnd iq: obtaining the electrical angle theta of the motor from the encoderrThen, the electrical angle theta is obtained through the formula (2)rWith respect to the differentiation of time, an electrical angular velocity ω is obtainedr(ii) a Three-phase stator current i at motor k moment is measured againa，ibAnd icObtaining i through a coordinate transformation moduledAnd iq；
Figure GDA0003170828780000071
(3) Calculating a predicted value psi of flux linkage at time (k +1)sd(k+1)、ψsqReference value ψ of flux linkage at time (k +1) and (k +1)sd ref(k+1)、ψsq refThe method of (k +1) is: the obtained d-q axis current component idAnd iqAngular velocity ω of rotorrAnd rotor electrical angle thetarInputting a model prediction flux linkage control module 2, obtaining a prediction current model at the moment (k +1) according to a formula (3), and then obtaining a predicted flux linkage value psi at the moment (k +1) according to a formula (4)sd(k+1)、ψsq(k + 1). Root of herbaceous plantObtaining the load angle delta and the electromagnetic torque T according to the formula (5)eAnd the load angle delta is derived according to the formula (6) to obtain the load angle increment delta of the formula (7), and the reference value delta of the load angle delta at the moment of (k +1) is obtained according to the formula (8)refThe reference value ψ of the flux linkage at the time (k +1) is obtained according to the formula (9)sd ref(k+1)、ψsq ref(k+1)；
Figure GDA0003170828780000072
Figure GDA0003170828780000073
Figure GDA0003170828780000074
Figure GDA0003170828780000075
Figure GDA0003170828780000076
δref＝Δδ+δ (8)
Figure GDA0003170828780000077
In the formula id k+1、iq k+1The predicted value of the current at the moment (k + 1); rsIs a stator phase resistance; l isd、LqThe inductor is a direct axis inductor and a quadrature axis inductor; t is the sampling period of the system; u. ofd k、uq kThe voltage of the stator voltage on the d-q axis component at the moment k; psifIs a rotor permanent magnet flux linkage; n ispIs the number of pole pairs; psisFor the resultant flux linkage psi on the d-q axis component at time ks(k) The amplitude of (d); psisd ref(k+1)、ψsq ref(k +1) is a flux linkage reference value of the stator flux linkage on the d-q axis component at the moment (k + 1); delta TeIs the electromagnetic torque increment; deltarefIs the reference value of the load angle at the moment (k + 1).
(4) Full-speed domain cost functions are constructed in full-speed domain cost function module 3, including cost functions of low and medium speed region control targets and cost functions of high-speed region control targets.
The method for constructing the low and medium speed region control target comprises the following steps: obtaining a flux linkage error function g of the low and medium speed region at the time of (k +1) according to the formula (10)F(ii) a When the permanent magnet synchronous motor operates in the mode of maximum torque current ratio, a convergence function g of a low-medium speed region at the moment of (k +1) is obtained according to a formula (11)MTPA(ii) a Obtaining a flux linkage limiting condition function g of the low and medium speed region at the moment (k +1) according to the formula (12)FMAX(ii) a Obtaining the direction selection function g of the low and medium speed regions at the moment (k +1) according to the formula (13)dir；
Figure GDA0003170828780000081
Figure GDA0003170828780000082
Figure GDA0003170828780000083
Figure GDA0003170828780000084
The method for constructing the high-speed area control target comprises the following steps: calculating a flux linkage error function g 'of the high-speed region at the moment (k + 1)'FFlux linkage error function g 'of high speed region'FFlux linkage error function g of low and medium speed regionFThe consistency is achieved; neglecting the stator resistance voltage drop when the permanent magnet synchronous motor operates stably at high speed to obtain a formula (14), and obtaining the formula (1)5) Obtaining a convergence function g of the high-speed area at the time of (k +1)FW(ii) a Calculating a high-speed region flux linkage limitation condition function g 'at the moment of (k + 1)'FMAXHigh speed region flux linkage constraint function g'FMAXFlux linkage limiting condition function g of low and medium speed regionsFMAXThe consistency is achieved; when the permanent magnet synchronous motor operates above the basic speed, the voltage limiting condition function g of the high-speed area at the moment of (k +1) is obtained according to the formula (16) under the constraint of the maximum output voltage of the inverterumax(ii) a Obtaining a high-speed region stable operation function g at the moment (k +1) according to the formula (17)stab；
Figure GDA0003170828780000091
Figure GDA0003170828780000092
Figure GDA0003170828780000093
Figure GDA0003170828780000094
In the formula usIs the stator voltage; u. ofsmaxIs the maximum output voltage of the inverter; vdcIs a dc bus voltage; lambda [ alpha ]mIs a voltage coefficient; eta is a voltage limiting condition intermediate variable; zeta is the intermediate variable of the motor high speed stable running condition.
(5) Designing a value function: obtaining a value function of the full-speed-domain permanent magnet synchronous motor according to the formula (18), respectively substituting eight basic voltage vectors in the table 1 into the value function, and outputting a switching state S enabling the value function to be minimumabcFeeding the inverter;
g(min)＝gF+gc+gL (18)
TABLE 1 basic Voltage vector Table
Figure GDA0003170828780000095
Definition of ωcFor the corresponding electrical angular velocity when the permanent magnet synchronous motor is operated at the base speed, when omegar<ωcWhen g isc＝gMTPAAnd g isL＝gFMAX+gdir(ii) a When ω isr>ωcWhen g isc＝gFWAnd g isL＝gFMAX+gumax+gstab。
(6) Calculating the duty ratio: the q-axis flux linkage obtained according to the formula (19) reaches the given value psi at the moment of (k +1) under the combined action of the optimal voltage vector and the zero voltage vectorsq ref(ii) a Combining a stator flux linkage equation and a voltage equation, and obtaining the slope S of the q-axis flux linkage when the zero vector acts according to the formula (20)0The slope S of the q-axis flux linkage at the time of the optimum voltage vector action is obtained from the equation (21)optObtaining the optimum voltage vector action time t according to the equation (22)opt；
Figure GDA0003170828780000101
Figure GDA0003170828780000102
Figure GDA0003170828780000103
Figure GDA0003170828780000104
In the formula, #sqRepresenting flux linkage, ψ, in the q-axis componentsq(k) Is the flux linkage on the q-axis component at time k; s0Is the slope of the q-axis flux linkage when zero vector is applied; soptIs the slope of the q-axis flux linkage when the optimal voltage vector acts; t is toptIs the optimal voltage vector action time;
Figure GDA0003170828780000105
representing the q-axis voltage under the optimal vector at time k.
The permanent magnet synchronous motor full-speed domain model prediction flux linkage control method is shown in a flow chart of fig. 2, and firstly, a stator current d-q axis component i at the moment k is obtaineddAnd iqElectric angle of rotor thetarAngular velocity ω of rotorrAnd a reference torque Te ref(ii) a And then predicting a flux linkage predicted value and a flux linkage reference value at the (k +1) moment, and predicting and constructing a flux linkage error function g at the (k +1) moment by utilizing three control requirements of a maximum torque-current ratio control or weak magnetic control strategy based on model prediction flux linkage controlFZone convergence function gcAnd a constraint function gL(ii) a Then, according to the rotating speed, selecting a value function of an MTPA (maximum Transmission Power Amplifier) area or a weak magnetic area and carrying out online rolling optimization to obtain an optimal voltage vector of the inverter; and finally, calculating the duty ratio by using the q-axis flux linkage dead beat, and distributing the time of the optimal voltage vector and the zero vector acting on the inverter.
FIG. 3 is a comparison of steady-state simulation results of model prediction flux linkage control of the single-double vector permanent magnet synchronous motor below the base speed. The simulation working condition is set as follows: the motor is given with the rotating speed of 400r/min and the electromagnetic torque of 4 N.m. At this time, (a) of fig. 3 predicts a torque ripple of flux linkage control to be 33% by the single vector model, and (b) of fig. 3 predicts a torque ripple of flux linkage control to be only 20% by the double vector model; comparing the current THD of the two schemes, the single vector model predicted flux linkage control current THD was 17.49%, while the dual vector model predicted flux linkage control current THD was only 7.96%. FIG. 4 is a comparison of simulation results of single and double vector permanent magnet synchronous motor model prediction flux linkage control when the rotation speed is suddenly changed below the base speed. Given a load of 4N · m, the initial rotational speed is 200r/min, and the rotational speed is abruptly changed to 400r/min at t ═ 0.6 s. For sudden change of the rotating speed, the two schemes can both rapidly respond, the rotating speed overshoot does not exceed 3%, and the current waveform can also keep good sine degree. The simulation results of the two schemes are compared to see that the dynamic performances of the single-double vector MPFC are basically consistent.
FIG. 5 shows the full speed rangePredicting a flux linkage control simulation result by a permanent magnet synchronous motor model, wherein the simulation working condition is set as follows: the motor is started from no load to the basic speed of 600r/min, then the flux weakening control is adopted, the rotating speed of the motor reaches 1800r/min, and the rotating speed is suddenly reduced to 0r/min when t is 0.4 s. FIG. 5 (a) shows a stator current isAnd d-q axis component id、iqA waveform diagram; fig. 5 (b) is a motor rotation speed map; FIG. 5 (c) is a three-phase current waveform diagram when the motor is operated in the full speed domain; fig. 5 (d) is a waveform diagram of electromagnetic torque of the motor. In the starting stage of the motor, in order to enable the motor to quickly reach a rotating speed set value 1800r/min, the three-phase current reaches a maximum value 10A, the torque reaches a maximum value 9N · m, the motor reaches the set rotating speed 1800r/min after being started for 5ms, the three-phase current value is 4A, the torque is 0N · m, when the rotating speed suddenly drops to 0r/min when the set rotating speed is t being 0.4s, the rotating speed of the motor becomes 0r/min after 10ms, and at the moment, the three-phase current and the torque respond quickly and reach the maximum value again.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and decorations can be made without departing from the principle of the present invention, and these modifications and decorations should also be regarded as the protection scope of the present invention.

Claims (3)

1. A permanent magnet synchronous motor full-speed domain model prediction flux linkage control method is characterized by comprising the following steps: the system comprises a rotating speed outer-loop PI controller (1), a model prediction flux linkage control module (2), a full-speed domain value function module (3), a duty ratio calculation module (4), an inverter (5), a coordinate transformation module (6), a permanent magnet synchronous motor (7) and an encoder (8); firstly, a reference torque T is obtained through a rotating speed outer ring PI controllere ref(ii) a Then obtaining the electrical angle theta of the permanent magnet synchronous motor from the motor encoderrAnd electrical angular velocity ωrAnd acquiring three-phase stator current i at the time k by using a current sensora，ibAnd icObtaining d-q axis component i of stator current at the moment k after coordinate transformationdAnd iq(ii) a Then, flux linkage calculation is carried out by combining a flux linkage equation and a load angle to obtain a predicted value psi of flux linkage at the (k +1) momentsd(k+1)、ψsq(k +1) and flux linkageReference value psisd ref(k+1)、ψsq ref(k + 1); further, a full-speed domain cost function is constructed by using the flux linkage predicted value and the flux linkage reference value at the (k +1) moment, and an optimal voltage vector of the inverter is obtained by minimizing the full-speed domain cost function; finally, duty ratio is calculated according to the q-axis flux linkage dead beat concept, and time of acting an optimal voltage vector and a zero vector on the inverter is distributed;
calculating a predicted value psi of flux linkage at time (k +1)sd(k+1)、ψsqReference value ψ of flux linkage at time (k +1) and (k +1)sd ref(k+1)、ψsq refThe method of (k +1) is: the obtained d-q axis current component idAnd iqAngular velocity ω of rotorrAnd rotor electrical angle thetarInputting a model prediction flux linkage control module, obtaining a prediction current model at the moment (k +1) according to a formula (3), and then obtaining a predicted flux linkage value psi at the moment (k +1) according to a formula (4)sd(k+1)、ψsq(k + 1); obtaining the load angle delta and the electromagnetic torque T according to the formula (5)eAnd the load angle delta is derived according to the formula (6) to obtain the load angle increment delta of the formula (7), and the reference value delta of the load angle delta at the moment of (k +1) is obtained according to the formula (8)refThe reference value ψ of the flux linkage at the time (k +1) is obtained according to the formula (9)sd ref(k+1)、ψsq ref(k+1)；
Figure FDA0003170828770000014
Figure FDA0003170828770000011
Figure FDA0003170828770000012
Figure FDA0003170828770000013
Figure FDA0003170828770000021
δref＝Δδ+δ (8)
Figure FDA0003170828770000022
In the formula id k+1、iq k+1The predicted value of the current at the moment (k + 1); rsIs a stator phase resistance; l isd、LqThe inductor is a direct axis inductor and a quadrature axis inductor; t is the sampling period of the system; u. ofd k、uq kThe voltage of the stator voltage on the d-q axis component at the moment k; psifIs a rotor permanent magnet flux linkage; n ispIs the number of pole pairs; psisFor the resultant flux linkage psi on the d-q axis component at time ks(k) The amplitude of (d); psisd ref(k+1)、ψsq ref(k +1) is a flux linkage reference value of the stator flux linkage on the d-q axis component at the moment (k + 1); delta TeIs the electromagnetic torque increment; deltarefA reference value of the load angle at the time (k + 1);
constructing a full-speed domain cost function comprises constructing a cost function of a low-medium speed region control target and a cost function of a high-speed region control target in a full-speed domain cost function module (3);
the method for constructing the low and medium speed region control target comprises the following steps: obtaining a flux linkage error function g of the low and medium speed region at the time of (k +1) according to the formula (10)F(ii) a When the permanent magnet synchronous motor operates in the mode of maximum torque current ratio, a convergence function g of a low-medium speed region at the moment of (k +1) is obtained according to a formula (11)MTPA(ii) a Obtaining a flux linkage limiting condition function g of the low and medium speed region at the moment (k +1) according to the formula (12)FMAX(ii) a Obtaining the direction selection function g of the low and medium speed regions at the moment (k +1) according to the formula (13)dir；
Figure FDA0003170828770000023
Figure FDA0003170828770000024
Figure FDA0003170828770000025
Figure FDA0003170828770000026
The method for constructing the high-speed area control target comprises the following steps: calculating a flux linkage error function g 'of the high-speed region at the moment (k + 1)'FFlux linkage error function g 'of high speed region'FFlux linkage error function g of low and medium speed regionFThe consistency is achieved; neglecting the stator resistance voltage drop when the permanent magnet synchronous motor operates stably at high speed to obtain a formula (14), and obtaining a high-speed area convergence function g at the moment of (k +1) according to the formula (15)FW(ii) a Calculating a high-speed region flux linkage limitation condition function g 'at the moment of (k + 1)'FMAXHigh speed region flux linkage constraint function g'FMAXFlux linkage limiting condition function g of low and medium speed regionsFMAXThe consistency is achieved; when the permanent magnet synchronous motor operates above the basic speed, the voltage limiting condition function g of the high-speed area at the moment of (k +1) is obtained according to the formula (16) under the constraint of the maximum output voltage of the inverterumax(ii) a Obtaining a high-speed region stable operation function g at the moment (k +1) according to the formula (17)stab；
Figure FDA0003170828770000035
Figure FDA0003170828770000031
Figure FDA0003170828770000032
Figure FDA0003170828770000033
In the formula usIs the stator voltage; u. ofsmaxIs the maximum output voltage of the inverter; vdcIs a dc bus voltage; lambda [ alpha ]mIs a voltage coefficient; eta is a voltage limiting condition intermediate variable; zeta is the intermediate variable of the motor high speed stable operation condition;
the method for obtaining the optimal voltage vector of the inverter by minimizing the full-speed domain cost function comprises the following steps: obtaining a full-speed domain cost function according to the formula (18), respectively substituting eight basic voltage vectors into the cost function, and outputting a switch state S which enables the cost function to be minimumabcFeeding the inverter;
g(min)＝gF+gc+gL (18)
definition of ωcFor the corresponding electrical angular velocity when the permanent magnet synchronous motor is operated at the base speed, when omegar<ωcWhen g isc＝gMTPAAnd g isL＝gFMAX+gdir(ii) a When ω isr>ωcWhen g isc＝gFWAnd g isL＝gFMAX+gumax+gstab；
The duty ratio calculation method comprises the following steps: the q-axis flux linkage obtained according to the formula (19) reaches the given value psi at the moment of (k +1) under the combined action of the optimal voltage vector and the zero voltage vectorsq ref(ii) a Combining a stator flux linkage equation and a voltage equation, and obtaining the slope S of the q-axis flux linkage when the zero vector acts according to the formula (20)0The slope S of the q-axis flux linkage at the time of the optimum voltage vector action is obtained from the equation (21)optObtaining the optimum voltage vector action time t according to the equation (22)opt：
Figure FDA0003170828770000034
Figure FDA0003170828770000041
Figure FDA0003170828770000042
Figure FDA0003170828770000043
In the formula, #sqRepresenting flux linkage, ψ, in the q-axis componentsq(k) Is the flux linkage on the q-axis component at time k; s0Is the slope of the q-axis flux linkage when zero vector is applied; soptIs the slope of the q-axis flux linkage when the optimal voltage vector acts; t is toptIs the optimal voltage vector action time;
Figure FDA0003170828770000044
representing the q-axis voltage under the optimal vector at time k.
2. The permanent magnet synchronous motor full-speed domain model prediction flux linkage control method according to claim 1, characterized in that: the reference torque Te refThe acquisition method comprises the following steps: the difference e between the reference speed and the actual speed of the motor is calculatednInputting a rotating speed PI controller, and obtaining a reference torque T according to a formula (1)e ref；
Figure FDA0003170828770000045
In the formula, kpAnd kiProportional gain and integral gain of the rotating speed PI controller are respectively shown, and s represents a complex variable.
3. The permanent magnet synchronous motor full-speed domain model prediction flux linkage control method according to claim 1, characterized in that: the electrical angle thetarElectrical angular velocity omegarAnd d-q axis component i of stator current at time kd，iqThe acquisition method comprises the following steps: obtaining the electrical angle theta of the permanent magnet synchronous motor from the encoderrThen, the electrical angle theta is obtained through the formula (2)rWith respect to the differentiation of time, an electrical angular velocity ω is obtainedr(ii) a Measuring k-time three-phase stator current i of permanent magnet synchronous motor by using current sensora，ibAnd icObtaining d-q axis component i of stator current at the moment k after coordinate transformationdAnd iq；
Figure FDA0003170828770000046
CN201911124614.XA 2019-11-18 2019-11-18 Permanent magnet synchronous motor full-speed domain model prediction flux linkage control method Active CN110784145B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911124614.XA CN110784145B (en) 2019-11-18 2019-11-18 Permanent magnet synchronous motor full-speed domain model prediction flux linkage control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911124614.XA CN110784145B (en) 2019-11-18 2019-11-18 Permanent magnet synchronous motor full-speed domain model prediction flux linkage control method

Publications (2)

Publication Number Publication Date
CN110784145A CN110784145A (en) 2020-02-11
CN110784145B true CN110784145B (en) 2021-11-26

Family

ID=69391454

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911124614.XA Active CN110784145B (en) 2019-11-18 2019-11-18 Permanent magnet synchronous motor full-speed domain model prediction flux linkage control method

Country Status (1)

Country Link
CN (1) CN110784145B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111342726B (en) * 2020-03-26 2023-05-30 武汉科技大学 Direct torque control method and system for permanent magnet synchronous motor
CN111431453A (en) * 2020-05-06 2020-07-17 南通大学 Fault-tolerant control method for open-winding permanent magnet synchronous motor based on single-phase fault of inverter
CN112910359A (en) * 2020-12-21 2021-06-04 南通斯迈尔精密设备有限公司 Improved permanent magnet synchronous linear motor model prediction current control method
CN112953358A (en) * 2021-03-12 2021-06-11 南通大学 Fault-tolerant control method for open-winding flux switching motor based on single-phase fault of inverter
CN114123904B (en) * 2021-06-10 2023-08-08 浙江大学先进电气装备创新中心 Predictive current increment control method suitable for operation of permanent magnet synchronous motor in high-speed region
CN115347829A (en) * 2022-08-12 2022-11-15 华北电力大学 Induction motor control method, system, device, electronic equipment and storage medium

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5387614B2 (en) * 2011-05-16 2014-01-15 株式会社デンソー Rotating machine control device
CN108448961B (en) * 2018-05-04 2019-09-03 南通大学 Meter and the permanent magnet synchronous motor model prediction method for controlling torque of switching frequency optimization
CN110460281B (en) * 2019-03-28 2021-04-13 南通大学 Three-level permanent magnet synchronous motor double-vector model prediction flux linkage control method

Also Published As

Publication number Publication date
CN110784145A (en) 2020-02-11

Similar Documents

Publication Publication Date Title
CN110784145B (en) Permanent magnet synchronous motor full-speed domain model prediction flux linkage control method
JP5357232B2 (en) Synchronous machine controller
JP5035641B2 (en) Control device for motor drive device
CN110460281B (en) Three-level permanent magnet synchronous motor double-vector model prediction flux linkage control method
CN109861609B (en) Five-bridge arm two-permanent magnet motor system optimization model prediction control device and method
JP5595835B2 (en) Electric motor drive
JP2000032799A (en) Controller and control method for electric rotating machine
CN111277180B (en) Rotating speed control method of square wave permanent magnet synchronous motor under two-axis rotating coordinate system
CN109600085B (en) Permanent magnet synchronous motor direct prediction duty ratio control method based on variable control set
WO2015056541A1 (en) Drive device for electric motor
CN111800056A (en) Permanent magnet synchronous motor three-vector model predicted torque control method based on novel switch table
CN109995293A (en) The switching method of I/F starting and closed-loop control under permanent magnet synchronous motor senseless control
CN112886880A (en) Three-level permanent magnet synchronous motor position sensorless model prediction current control method
Dahmardeh et al. A novel combined DTC method and SFOC system for three-phase induction machine drives with PWM switching method
CN110474589B (en) Control method for optimizing efficiency of induction motor
CN110096077B (en) Nonsingular rapid terminal sliding mode rotating speed control method and system for switched reluctance motor
CN110529979A (en) A kind of motor control method, device and air conditioner
Dang et al. Model Predictive Direct Torque Control for PMSM Drives in M–T Frame
CN116526919A (en) Permanent magnet synchronous motor servo system and current prediction control method and device thereof
CN115189610A (en) Multi-mode current prediction control method for hybrid excitation axial magnetic field permanent magnet motor
CN112234894B (en) Dead-beat direct torque-flux linkage control system and method for variable flux memory motor
JP5326444B2 (en) Rotating machine control device
CN109660176B (en) Control method of synchronous reluctance motor for new energy vehicle based on alternating current stepping control
CN107482961B (en) Brushless double-fed motor torque ripple minimization control method
JP2003111490A (en) Method and apparatus of controlling inverter

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
TR01 Transfer of patent right

Effective date of registration: 20230801

Address after: 509 Kangrui Times Square, Keyuan Business Building, No. 39 Huarong Road, Gaofeng Community, Dalang Street, Longhua District, Shenzhen City, Guangdong Province, 518109

Patentee after: Shenzhen lizhuan Technology Transfer Center Co.,Ltd.

Address before: 226019 Jiangsu Province, Nantong City Chongchuan District sik Road No. 9

Patentee before: NANTONG University

Effective date of registration: 20230801

Address after: 522021 Road 3, Penglin Industrial Zone, Rongdong Office, Rongcheng District, Jieyang City, Guangdong Province

Patentee after: Jieyang Shengjia Motor Co.,Ltd.

Address before: 509 Kangrui Times Square, Keyuan Business Building, No. 39 Huarong Road, Gaofeng Community, Dalang Street, Longhua District, Shenzhen City, Guangdong Province, 518109

Patentee before: Shenzhen lizhuan Technology Transfer Center Co.,Ltd.

TR01 Transfer of patent right