CN115347841A - Dead-beat prediction current loop control method for permanent magnet synchronous motor - Google Patents

Dead-beat prediction current loop control method for permanent magnet synchronous motor Download PDF

Info

Publication number
CN115347841A
CN115347841A CN202211159454.4A CN202211159454A CN115347841A CN 115347841 A CN115347841 A CN 115347841A CN 202211159454 A CN202211159454 A CN 202211159454A CN 115347841 A CN115347841 A CN 115347841A
Authority
CN
China
Prior art keywords
current
speed
modulation period
observed value
permanent magnet
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.)
Pending
Application number
CN202211159454.4A
Other languages
Chinese (zh)
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.)
State Grid Heilongjiang Electric Power Co Ltd Electric Power Research Institute
Harbin Institute of Technology
State Grid Corp of China SGCC
Original Assignee
State Grid Heilongjiang Electric Power Co Ltd Electric Power Research Institute
Harbin Institute of Technology
State Grid Corp of China SGCC
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 State Grid Heilongjiang Electric Power Co Ltd Electric Power Research Institute, Harbin Institute of Technology, State Grid Corp of China SGCC filed Critical State Grid Heilongjiang Electric Power Co Ltd Electric Power Research Institute
Priority to CN202211159454.4A priority Critical patent/CN115347841A/en
Publication of CN115347841A publication Critical patent/CN115347841A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/22Current control, e.g. using a current control loop
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • H02P21/0007Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control using sliding mode control
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/13Observer control, e.g. using Luenberger observers or Kalman filters
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18Estimation of position or speed
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements 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/022Synchronous motors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05Synchronous machines, e.g. with permanent magnets or DC excitation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

A dead-beat prediction current loop control method for a permanent magnet synchronous motor belongs to the technical field of motor control. The invention solves the problems of parameter mismatching in the current loop and strong interference of parameters and loads in the speed loop in the existing compensation process. The method comprises the steps of establishing a mathematical model of the permanent magnet synchronous motor, and obtaining a current loop supercoiled sliding-mode observer by using the mathematical model; establishing a motor motion equation with disturbance, and designing a speed ring high-order sliding mode observer; observing the speed ring disturbance of the motor by using a speed ring supercoiled sliding-mode observer to obtain a disturbance observation result; improving a sliding mode approximation law by using the actual rotating speed and the given value of the rotating speed of the motor to obtain a speed control law; and inputting the observed value and the speed control law of the current loop supercoiled sliding-mode observer into a dead-beat current controller to obtain a current loop control signal. The invention is suitable for dead-beat prediction current control of the permanent magnet synchronous motor.

Description

Dead-beat prediction current loop control method for permanent magnet synchronous motor
Technical Field
The invention belongs to the technical field of motor control.
Background
In a permanent magnet synchronous motor driving system, the traditional control method mainly comprises vector control and direct torque control. However, for the rapidly developing industrial field, the traditional control strategy cannot meet the continuously developing industrial demand, so that some new modern control methods are proposed by the scholars one after another. Model predictive control has been successful in complex industrial processes since the advent, and has evolved from the original heuristic control algorithm to a new branch of discipline in the industrial field. Dead-beat predictive current control is a typical model-based method, and the mismatch of system models can seriously degrade the performance of the controller. The method is an effective method for estimating disturbance and uncertainty variables by designing an observer and then compensating the estimated disturbance into a predictive controller, but the existing compensation process has the problems of parameter mismatch in a current loop and strong disturbance of parameters and loads in a speed loop.
Disclosure of Invention
The invention aims to solve the problems of parameter mismatch in a current loop and strong interference of parameters and loads in a speed loop in the conventional compensation process, and provides a dead-beat prediction current control method for a permanent magnet synchronous motor.
The invention discloses a dead-beat prediction current control method for a permanent magnet synchronous motor, which comprises the following steps of:
step one, establishing a mathematical model of a permanent magnet synchronous motor, and obtaining a current loop supercoiled sliding-mode observer by using the mathematical model;
establishing a motor motion equation with disturbance, and designing a speed ring high-order sliding mode observer;
observing the speed ring disturbance of the motor by using a speed ring supercoiled sliding-mode observer to obtain a disturbance observation result;
improving a sliding mode approach law by utilizing the actual rotating speed and a rotating speed given value of the motor, and compensating the disturbance observation result into the improved sliding mode approach law to obtain a speed control law;
and fourthly, observing the current loop by using the current loop supercoiled sliding-mode observer, inputting the observed value and the speed control law of the current loop supercoiled sliding-mode observer to the deadbeat current controller, and acquiring a current loop control signal.
Further, in the present invention, the method for establishing the mathematical model of the permanent magnet synchronous motor comprises:
establishing a stator current equation of the permanent magnet synchronous motor under a dq coordinate system:
Figure BDA0003858952520000021
wherein: u. of d ,u q D-axis component and q-axis component of stator voltage under dq coordinate system respectively; i.e. i d ,i q D-axis component and q-axis component of stator current under dq coordinate system respectively; r is a stator resistor; l is a radical of an alcohol d ,L q Stator inductances of d and q axes respectively; omega e Is the angular velocity of the motor; psi f The amplitude of the flux linkage of the permanent magnet of the rotor;
discretizing the stator current equation by adopting a forward Euler method, and controlling an actual current vector i (k + 1) to reach a reference current value i after one modulation period * (k+1),i(k+1)=i * (k + 1), obtaining a discrete current model of the permanent magnet synchronous motor:
Figure BDA0003858952520000022
wherein u is d (k) Is the d-axis component of the kth modulation period of the stator voltage in the dq coordinate system; u. of q (k) Is the q-axis component of the kth modulation period of the stator voltage in the dq coordinate system; t is s Is a sampling period, i d (k) Is the d-axis component, i, of the k-th modulation period of the stator current in dq coordinate system d (k) Is the q-axis component of the k-th modulation period of the stator current in the dq coordinate system,
Figure BDA0003858952520000023
is a given value of a d-axis component of a (k + 1) th modulation period of the stator current under a dq coordinate system,
Figure BDA0003858952520000024
is a given value of a q-axis component of a k +1 th modulation period of the stator current in a dq coordinate system.
Further, in the present invention, in the first step, the method for obtaining the current-loop supercoiled sliding-mode observer by using the mathematical model includes:
establishing a voltage equation containing parameter uncertainty by using the mathematical model, and acquiring an expression of the current loop superspiral sliding-mode observer according to the voltage equation containing parameter uncertainty, wherein the expression is as follows:
Figure BDA0003858952520000025
Figure BDA0003858952520000031
wherein the content of the first and second substances,
Figure BDA0003858952520000032
is i d Is detected by the measured values of (a) and (b),
Figure BDA0003858952520000033
is i q Observed value of (a), k 1 Is a first control parameter, k, of a d-axis current observer 2 Is a second control parameter, k, of the d-axis current observer 3 A third control parameter that is a d-axis current observer; k is a radical of 4 Is a first control parameter, k, of a q-axis current observer 5 Is a second control parameter, k, of the q-axis current observer 6 A third control parameter, F, for the q-axis current observer d ,F q Are respectively f d ,f q The derivative of (a) of (b),
Figure BDA0003858952520000034
are respectively as
Figure BDA0003858952520000035
The derivative of (a) of (b),
Figure BDA0003858952520000036
is F q Is detected by the measured values of (a) and (b),
Figure BDA0003858952520000037
is F d Is detected by the measured values of (a) and (b),
Figure BDA0003858952520000038
is F q The derivative of (a) of (b),
Figure BDA0003858952520000039
is composed of
Figure BDA00038589525200000310
Sgn () is a sign function,
Figure BDA00038589525200000311
is f d Is detected by the measured values of (a) and (b),
Figure BDA00038589525200000312
is f q Observed value of f d 、f q Perturbation of the d-axis and q-axis parameters, respectively.
Further, the invention also includes a step of discretizing the current loop supercoiled sliding mode observer, specifically:
by utilizing a forward Euler method, the errors of observed current and disturbance are converged to zero in a limited time, and a discrete time equation of a current loop supercoiled sliding-mode observer is as follows:
Figure BDA00038589525200000313
Figure BDA00038589525200000314
wherein k is the number of modulation cycles, T s Which is indicative of the period of the modulation,
Figure BDA00038589525200000315
is i d The observed value at the k +1 modulation period,
Figure BDA00038589525200000316
is i q The observed value at the k +1 th modulation period,
Figure BDA00038589525200000317
is f d The observed value at the k-th modulation period,
Figure BDA00038589525200000318
is f d The observed value at the k +1 modulation period,
Figure BDA00038589525200000319
is F d The observed value at the k +1 th modulation period,
Figure BDA00038589525200000320
is F d The observed value at the k-th modulation period,
Figure BDA00038589525200000321
is f q The observed value at the k-th modulation period,
Figure BDA00038589525200000322
is f q The observed value at the k +1 th modulation period,
Figure BDA0003858952520000041
is F q The observed value at the k +1 th modulation period,
Figure BDA0003858952520000042
is F q The observed value at the k-th modulation period,
Figure BDA0003858952520000043
is i d Observed value at k modulation period.
Further, in the present invention, in the first step, the motor motion equation with disturbance is:
Figure BDA0003858952520000044
wherein, P n Is the number of magnetic pole pairs; omega m The mechanical angular speed of the motor; b is a friction coefficient; j is the moment of inertia, and ρ is the parametric perturbation and the total perturbation of the external load.
Further, in the present invention, in the first step, the speed loop high-order sliding mode observer is:
Figure BDA0003858952520000045
in the formula (I), the compound is shown in the specification,
Figure BDA0003858952520000046
in order to observe the error in the speed,
Figure BDA0003858952520000047
is the observed value of p, D is the derivative of p,
Figure BDA0003858952520000048
is an observed value of D, k w1 Is a first control parameter, k, of a high-order sliding-mode observer of the speed loop w2 Is a second control parameter, k, of a speed-loop high-order sliding-mode observer w3 Is a third control parameter of the speed loop high-order sliding mode observer,
Figure BDA0003858952520000049
for mechanical angular velocity omega of the motor m The observed value of (1).
Further, in the present invention, in the first step, a step of discretizing a high-order sliding mode observer of the velocity ring is further included, specifically:
Figure BDA00038589525200000410
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA00038589525200000411
is the observed value, omega, of the motor mechanical angular velocity at the k +1 th modulation period m (k) Is the observed value of the k modulation period of the mechanical angular speed of the motor,
Figure BDA00038589525200000412
is observed value at the k modulation period.
Further, in the third step, the method for improving the sliding mode approach law by using the actual rotating speed and the given value of the rotating speed of the motor comprises the following steps:
defining a slip form surface:
Figure BDA00038589525200000413
wherein e is w Is the error of the rotating speed; c is an integral coefficient of an integral sliding mode surface; s is a slip form surface;
Figure BDA00038589525200000414
in the formula (I), the compound is shown in the specification,
Figure BDA0003858952520000051
setting a rotating speed value;
the improved sliding mode approach law is as follows:
Figure BDA0003858952520000052
in the formula (I), the compound is shown in the specification,
Figure BDA0003858952520000053
the derivative of slip form surface s, k s For switching the gain, ε is a variable term ε + (1- ε) e -δ|s| Gain of (k) t In order to obtain a linear gain, the gain is, δ is an exponential term e -δ|s| And k is s >0,k t Greater than 0, delta greater than 0,0 < epsilon < 1, when the system is far away from the sliding mode surface, the coefficient in the front of the sign function approaches k s |e|/ε,
Figure BDA0003858952520000054
A value of greater than k s The convergence speed is accelerated, and when the system reaches the sliding mode, the coefficient in front of the sign function approaches to k s |e w |,
Figure BDA0003858952520000055
Value of (A)Less than k s And buffeting is limited.
Further, in the third step, the method for obtaining the speed control law includes:
compensating the observed parameter perturbation and the total disturbance rho of the external load into a control law, wherein the obtained speed control law is as follows:
Figure BDA0003858952520000056
Figure BDA0003858952520000057
the set value of the mechanical angular speed of the motor,
Figure BDA0003858952520000058
is composed of
Figure BDA0003858952520000059
The derivative of (c).
The high-order (third order, more than second order defined in the field is high order) supercoiled sliding mode observer designed by the invention can estimate the concentrated disturbance caused by the mismatch of future current values and parameters, effectively eliminate the influence of the mismatch of the parameters and improve the anti-interference performance of a current loop. In addition, a speed controller is designed by adopting an improved sliding mode approach law in a speed ring, and a third-order supercoiled sliding mode observer is designed to estimate parameter perturbation and compensate total disturbance of an external load into the speed controller, so that a stable current reference value is provided for a current ring.
Drawings
FIG. 1 is a schematic block diagram of dead-beat prediction current control of a permanent magnet synchronous motor based on a high-order sliding-mode observer according to the present invention;
FIG. 2 is a graph of PI control for the speed loop, 10 times the nominal resistance parameter for the current loop, and conventional dead-beat predictive current control for the current loop d ,i q A current effect graph;
FIG. 3 shows the speed loop using PI control, the current loop using conventional dead-beat predictive current control with 2 times the nominal inductance parameter, i d ,i q A current effect graph;
FIG. 4 is a graph of PI control for the speed loop, 4 times the flux linkage parameter to the nominal parameter for the current loop using conventional dead-beat predictive current control d ,i q A current effect graph;
FIG. 5 shows that the speed loop adopts PI control, the current loop adopts dead-beat prediction current control combined with a high-order sliding-mode observer, and the resistance parameter is changed to 10 times of the nominal parameter, i d ,i q A current effect graph;
FIG. 6 shows that the speed loop adopts PI control, the current loop adopts deadbeat prediction current control combined with a high-order sliding-mode observer, and inductance parameters are changed to 2 times, i, of nominal parameters d ,i q A current effect graph;
FIG. 7 shows that the speed loop adopts PI control, the current loop adopts dead-beat prediction current control combined with a high-order sliding-mode observer, and flux linkage parameters are changed to be 4 times, i, of nominal parameters d ,i q A current effect graph;
FIG. 8 is a rotation speed diagram of an improved sliding mode approach law speed control method adopted;
FIG. 9 is a rotation speed diagram of a control method combining an improved sliding mode approach law with a high-order sliding mode observer.
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 efforts based on the embodiments of the present invention, shall fall within the scope of protection of the present invention.
It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict.
The first embodiment is as follows: the present embodiment is described below with reference to fig. 1, and a method for dead-beat prediction current control of a permanent magnet synchronous motor according to the present embodiment includes:
establishing a mathematical model of a permanent magnet synchronous motor, and acquiring a current loop supercoiled sliding-mode observer by using the mathematical model;
establishing a motor motion equation with disturbance, and designing a speed ring high-order sliding mode observer;
observing the speed ring disturbance of the motor by using a speed ring supercoiled sliding-mode observer to obtain a disturbance observation result;
improving a sliding mode approach law by using the actual rotating speed and the given value of the rotating speed of the motor, and compensating the disturbance observation result to the improved sliding mode approach law to obtain a speed control law;
and fourthly, observing the current loop by using the current loop supercoil sliding mode observer, inputting the observed value and the speed control law of the current loop supercoil sliding mode observer into the dead-beat current controller, and acquiring a current loop control signal.
The invention improves the speed and current anti-interference performance and tracking precision of the permanent magnet synchronous motor. The invention regards parameter change and external load change as lumped disturbance, establishes a mathematical model of the built-in permanent magnet synchronous motor, further constructs two third-order supercoiled sliding-mode observers to respectively estimate the lumped disturbance in a speed loop and a current loop, and feeds the estimated disturbance forward to a corresponding controller for compensation, so as to improve the robustness and the tracking precision of a system current loop. In addition, in order to obtain a stable current loop set value, an improved sliding mode approach law is adopted in the rotating speed control, so that buffeting in a control signal is reduced, the time required by a system state to reach a sliding mode surface is shortened, and the anti-interference performance is further improved. The experimental results prove the effectiveness of the invention. In the invention, a current loop control signal is obtained and input to a current controller, and the current controller generates a voltage control signal
Figure BDA0003858952520000071
Is transformed into by coordinate
Figure BDA0003858952520000072
Entering SVPMW modulation module, generating control pulse to inverter module to driveThe permanent magnet synchronous motor rotates.
Further, in the present invention, the method for establishing the mathematical model of the permanent magnet synchronous motor comprises:
establishing a stator current equation of the permanent magnet synchronous motor under a dq coordinate system:
Figure BDA0003858952520000073
wherein: u. of d ,u q Respectively a d-axis component and a q-axis component of the stator voltage under the dq coordinate system; i all right angle d ,i q D-axis component and q-axis component of stator current under dq coordinate system respectively; r is a stator resistor; l is a radical of an alcohol d ,L q Stator inductances of d and q axes respectively; omega e Is the angular velocity of the motor; psi f The amplitude of the flux linkage of the permanent magnet of the rotor;
discretizing the stator current equation by adopting a forward Euler method, and controlling an actual current vector i (k + 1) to reach a reference current value i after a modulation period * (k+1),i(k+1)=i * (k + 1), obtaining a discrete current model of the permanent magnet synchronous motor:
Figure BDA0003858952520000074
wherein u is d (k) Is the d-axis component of the kth modulation period of the stator voltage in the dq coordinate system; u. u q (k) Is the q-axis component of the kth modulation period of the stator voltage under the dq coordinate system; t is a unit of s Is a sampling period, i d (k) Is the d-axis component, i, of the k-th modulation period of the stator current in dq coordinate system d (k) Is the q-axis component of the k-th modulation period of the stator current in the dq coordinate system,
Figure BDA0003858952520000075
is a given value of a d-axis component of a (k + 1) th modulation period of the stator current under a dq coordinate system,
Figure BDA0003858952520000076
is the dq coordinate systemThe given value of the q-axis component of the (k + 1) th modulation period of the lower stator current.
Further, in the present invention, in the first step, the method for obtaining the current-loop supercoiled sliding-mode observer by using the mathematical model includes:
establishing a voltage equation containing parameter uncertainty by using the mathematical model, and obtaining an expression of the current loop supercoiled sliding mode observer according to the voltage equation containing parameter uncertainty, wherein the expression is as follows:
Figure BDA0003858952520000081
Figure BDA0003858952520000082
wherein the content of the first and second substances,
Figure BDA0003858952520000083
is i d Is measured in a time-domain manner by a time-domain,
Figure BDA0003858952520000084
is i q Observed value of (a), k 1 Is a first control parameter, k, of a d-axis current observer 2 Is a second control parameter, k, of the d-axis current observer 3 A third control parameter that is a d-axis current observer; k is a radical of formula 4 Is a first control parameter, k, of a q-axis current observer 5 As a second control parameter, k, of the q-axis current observer 6 A third control parameter, F, of the q-axis current observer d ,F q Are respectively f d ,f q The derivative of (a) is determined,
Figure BDA0003858952520000085
are respectively as
Figure BDA0003858952520000086
The derivative of (a) of (b),
Figure BDA0003858952520000087
is F q Is detected by the measured values of (a) and (b),
Figure BDA0003858952520000088
is F d Is detected by the measured values of (a) and (b),
Figure BDA0003858952520000089
is F q The derivative of (a) of (b),
Figure BDA00038589525200000810
is composed of
Figure BDA00038589525200000811
Sgn () is a sign function,
Figure BDA00038589525200000812
is f d Is measured in a time-domain manner by a time-domain,
Figure BDA00038589525200000813
is f q Observed value of (a), f d 、f q Perturbation of the d-axis and q-axis parameters, respectively.
Further, the invention also includes a step of discretizing the current loop supercoiled sliding mode observer, specifically:
by utilizing a forward Euler method, the errors of observed current and disturbance are converged to zero in a limited time, and a discrete time equation of the current loop superspiral sliding-mode observer is as follows:
Figure BDA00038589525200000814
Figure BDA0003858952520000091
wherein k is the number of modulation periods, T s Which is indicative of the period of the modulation,
Figure BDA0003858952520000092
is i d The observed value at the k +1 th modulation period,
Figure BDA0003858952520000093
is i q The observed value at the k +1 th modulation period,
Figure BDA0003858952520000094
is f d The observed value at the k-th modulation period,
Figure BDA0003858952520000095
is f d The observed value at the k +1 th modulation period,
Figure BDA0003858952520000096
is F d The observed value at the k +1 modulation period,
Figure BDA0003858952520000097
is F d The observed value at the k-th modulation period,
Figure BDA0003858952520000098
is f q The observed value at the k-th modulation period,
Figure BDA0003858952520000099
is f q The observed value at the k +1 th modulation period,
Figure BDA00038589525200000910
the observed value at the k +1 th modulation period,
Figure BDA00038589525200000911
is F q The observed value at the k-th modulation period,
Figure BDA00038589525200000912
is i d Observed value at k modulation period.
Further, in the present invention, in the first step, the motor motion equation with disturbance is:
Figure BDA00038589525200000913
wherein, P n The number of magnetic pole pairs is shown; omega m The mechanical angular speed of the motor; b is a friction coefficient; j is the moment of inertia, and ρ is the parameter perturbation and the total disturbance of the external load.
Further, in the present invention, in the first step, the speed loop high-order sliding mode observer is:
Figure BDA00038589525200000914
in the formula (I), the compound is shown in the specification,
Figure BDA00038589525200000915
in order to observe the error in the speed,
Figure BDA00038589525200000916
is the observed value of p, D is the derivative of p,
Figure BDA00038589525200000917
is an observed value of D, k w1 Is a first control parameter, k, of a high-order sliding mode observer of a speed ring w2 Is a second control parameter, k, of a speed-loop high-order sliding-mode observer w3 Is a third control parameter of the speed loop high-order sliding mode observer,
Figure BDA00038589525200000918
for mechanical angular velocity omega of the motor m The observed value of (a).
Further, in the present invention, in the first step, a step of discretizing the speed ring high-order sliding-mode observer is further included, specifically:
Figure BDA0003858952520000101
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003858952520000102
is the observed value, omega, of the motor mechanical angular velocity at the k +1 th modulation period m (k) Is the observed value of the k modulation period of the mechanical angular speed of the motor,
Figure BDA0003858952520000103
is observed value at the k modulation period.
Further, in the third step, the method for improving the sliding mode approach law by using the actual rotating speed and the given value of the rotating speed of the motor comprises the following steps:
defining a slip form surface:
Figure BDA0003858952520000104
wherein e is w Is the error of the rotating speed; c is an integral coefficient of an integral sliding mode surface; s is a slip form surface;
Figure BDA0003858952520000105
in the formula (I), the compound is shown in the specification,
Figure BDA0003858952520000106
setting a rotating speed value;
the improved sliding mode approach law is as follows:
Figure BDA0003858952520000107
in the formula (I), the compound is shown in the specification,
Figure BDA0003858952520000108
the derivative of the slip-form surface s, k s For switched gain, ε is the variable term ε + (1- ε) e -δ|s| Gain of (k) t For linear gain, δ is an exponential term e -δ|s| And k is s >0,k t Greater than 0, delta greater than 0,0 < epsilon < 1, when the system is far away from the sliding mode surface, the coefficient in front of the sign function approaches k s |e|/ε,
Figure BDA0003858952520000109
A value of greater than k s The convergence speed is accelerated, and when the system reaches the sliding mode, the coefficient in front of the sign function approaches to k s |e w |,
Figure BDA00038589525200001010
Value of less than k s And buffeting limitation is realized.
Further, in the third step, the method for obtaining the speed control law includes:
compensating the observed parameter perturbation and the total disturbance rho of the external load into a control law, wherein the obtained speed control law is as follows:
Figure BDA00038589525200001011
Figure BDA00038589525200001012
the set value of the mechanical angular speed of the motor,
Figure BDA00038589525200001013
is composed of
Figure BDA00038589525200001014
The derivative of (c).
In the invention, a stator current equation of the built-in permanent magnet synchronous motor in a dq coordinate system is adopted in the first step as follows:
Figure BDA0003858952520000111
electromagnetic torque equation:
Figure BDA0003858952520000112
equation of motion:
Figure BDA0003858952520000113
wherein: u. u d ,u q ;i d ,i q Stator voltage and stator current under dq coordinate system respectively; r is a stator resistor; l is d ,L q Dq-axis inductances, respectively; omega e Is the electrical angular velocity; psi f The amplitude of the flux linkage of the permanent magnet of the rotor; t is a unit of e Is the electromagnetic torque; p is n The number of magnetic pole pairs is shown; omega m Is the mechanical angular velocity; t is L Is the load torque; b is a friction coefficient; j is moment of inertia.
The forward euler method is applied to discretize the model shown in (1), and the discrete current model of the permanent magnet synchronous motor can be represented as follows:
Figure BDA0003858952520000114
wherein, T s Is the sampling period. According to a discrete current model of a discrete permanent magnet synchronous motor, after one modulation period, an actual current vector reaches a reference current, namely i (k + 1) = i * (k + 1), the stator voltage expression at this time is as follows:
Figure BDA0003858952520000115
the parameters set in the controller may be different from actual as the motor operates or external conditions change. The voltage equation containing the parameter uncertainty is expressed as:
Figure BDA0003858952520000116
wherein f is d 、f q Is a parametric perturbation.
The effectiveness of the invention is demonstrated by experimental results, the system parameters of the invention are shown in Table 1, and the overall control block diagram of the invention is shown in Table 1As shown in fig. 1. The motor is started without load, the given rotating speed is 1500r/min, and the load torque is suddenly changed into 12 N.m in the running process. As can be seen from a comparison between fig. 2 and fig. 7, the change in resistance does not significantly affect the current of the system, but when the inductance and flux linkage change, the current ripple of the conventional control method becomes large, and the tracking accuracy decreases. And most obvious when the magnetic linkage changes, current i q Fully-tracking not-on reference current i qref . The system designed by the invention shows better dynamic and steady-state performances when the resistance, the inductance and the flux linkage change. As can be seen from fig. 8 and 9, the system designed by the present invention has a shorter regulation time of the rotation speed and a smaller rotation speed drop when the load changes. Experiments show that the system designed by the invention has better dynamic and steady-state performance.
TABLE 1
Figure BDA0003858952520000121
Although the invention herein has been described with reference to particular embodiments, it is to be understood that these embodiments are merely illustrative of the principles and applications of the present invention. It is therefore to be understood that numerous modifications may be made to the illustrative embodiments and that other arrangements may be devised without departing from the spirit and scope of the present invention as defined by the appended claims. It should be understood that features described in different dependent claims and herein may be combined in ways different from those described in the original claims. It is also to be understood that features described in connection with individual embodiments may be used in other described embodiments.

Claims (9)

1. A dead-beat prediction current loop control method for a permanent magnet synchronous motor is characterized by comprising the following steps:
establishing a mathematical model of a permanent magnet synchronous motor, and acquiring a current loop supercoiled sliding-mode observer by using the mathematical model;
establishing a motor motion equation with disturbance, and designing a speed ring high-order sliding mode observer;
observing the speed ring disturbance of the motor by using a speed ring supercoiled sliding-mode observer to obtain a disturbance observation result;
improving a sliding mode approach law by using the actual rotating speed and the given value of the rotating speed of the motor, and compensating the disturbance observation result to the improved sliding mode approach law to obtain a speed control law;
and fourthly, observing the current loop by using the current loop supercoiled sliding-mode observer, inputting the observed value and the speed control law of the current loop supercoiled sliding-mode observer to the deadbeat current controller, and acquiring a current loop control signal.
2. The dead-beat prediction current loop control method of the permanent magnet synchronous motor according to claim 1, wherein the method for establishing the mathematical model of the permanent magnet synchronous motor comprises the following steps:
establishing a stator current equation of the permanent magnet synchronous motor in a dq coordinate system:
Figure FDA0003858952510000011
wherein: u. of d ,u q Respectively a d-axis component and a q-axis component of the stator voltage under the dq coordinate system; i all right angle d ,i q D-axis component and q-axis component of stator current under dq coordinate system respectively; r is a stator resistor; l is d ,L q Stator inductances of d and q axes respectively; omega e Is the angular velocity of the motor; psi f The amplitude of the flux linkage of the permanent magnet of the rotor;
discretizing the stator current equation by adopting a forward Euler method, and controlling an actual current vector i (k + 1) to reach a reference current value i after a modulation period * (k+1),i(k+1)=i * (k + 1), acquiring a discrete current model of the permanent magnet synchronous motor:
Figure FDA0003858952510000012
wherein u is d (k) Is the d-axis component of the kth modulation period of the stator voltage in the dq coordinate system; u. u q (k) Is the q-axis component of the kth modulation period of the stator voltage under the dq coordinate system; t is a unit of s Is a sampling period, i d (k) Is the d-axis component, i, of the kth modulation period of the stator current in dq coordinate system d (k) Is the q-axis component of the k-th modulation period of the stator current in the dq coordinate system,
Figure FDA0003858952510000013
is a given value of a d-axis component of a k +1 th modulation period of the stator current in a dq coordinate system,
Figure FDA0003858952510000014
is a given value of the q-axis component of the k +1 th modulation period of the stator current in the dq coordinate system.
3. The method for controlling the dead-beat prediction current loop of the permanent magnet synchronous motor according to claim 2, wherein in the first step, the method for obtaining the current loop supercoiled sliding-mode observer by using the mathematical model comprises the following steps:
establishing a voltage equation containing parameter uncertainty by using the mathematical model, and obtaining an expression of the current loop supercoiled sliding mode observer according to the voltage equation containing parameter uncertainty, wherein the expression is as follows:
Figure FDA0003858952510000021
Figure FDA0003858952510000022
wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003858952510000023
Figure FDA0003858952510000024
is i d Is detected by the measured values of (a) and (b),
Figure FDA0003858952510000025
is i q Observed value of (a), k 1 Is a first control parameter, k, of a d-axis current observer 2 Second control parameter, k, of d-axis current observer 3 A third control parameter of the d-axis current observer; k is a radical of 4 Is a first control parameter, k, of a q-axis current observer 5 Is a second control parameter, k, of the q-axis current observer 6 A third control parameter, F, for the q-axis current observer d ,F q Are respectively f d ,f q The derivative of (a) of (b),
Figure FDA0003858952510000026
are respectively as
Figure FDA0003858952510000027
The derivative of (a) of (b),
Figure FDA0003858952510000028
is F q Is detected by the measured values of (a) and (b),
Figure FDA0003858952510000029
is F d Is detected by the measured values of (a) and (b),
Figure FDA00038589525100000210
is F q The derivative of (a) of (b),
Figure FDA00038589525100000211
is composed of
Figure FDA00038589525100000212
Sgn () is a sign function,
Figure FDA00038589525100000213
is f d Is measured in a time-domain manner by a time-domain,
Figure FDA00038589525100000214
is f q Observed value of (a), f d 、f q Perturbation of the d-axis and q-axis parameters, respectively.
4. The dead-beat prediction current loop control method of the permanent magnet synchronous motor according to claim 3, further comprising a step of discretizing a current loop supercoiled sliding-mode observer, specifically:
by utilizing a forward Euler method, the errors of observed current and disturbance are converged to zero in a limited time, and a discrete time equation of a current loop supercoiled sliding-mode observer is as follows:
Figure FDA00038589525100000215
Figure FDA0003858952510000031
wherein k is the number of modulation periods, T s Which is indicative of the period of the modulation,
Figure FDA0003858952510000032
is i d The observed value at the k +1 th modulation period,
Figure FDA0003858952510000033
is i q The observed value at the k +1 th modulation period,
Figure FDA0003858952510000034
is f d The observed value at the k-th modulation period,
Figure FDA0003858952510000035
is f d The observed value at the k +1 th modulation period,
Figure FDA0003858952510000036
is F d The observed value at the k +1 th modulation period,
Figure FDA0003858952510000037
is F d The observed value at the k-th modulation period,
Figure FDA0003858952510000038
is f q The observed value at the k-th modulation period,
Figure FDA0003858952510000039
is f q The observed value at the k +1 th modulation period,
Figure FDA00038589525100000310
is F q The observed value at the k +1 th modulation period,
Figure FDA00038589525100000311
is F q The observed value at the k-th modulation period,
Figure FDA00038589525100000312
is i d Observed value at k modulation period.
5. The dead-beat prediction current loop control method for the permanent magnet synchronous motor according to claim 3 or 4, characterized in that a motor motion equation with disturbance is as follows:
Figure FDA00038589525100000313
wherein, P n Is the number of magnetic pole pairs; omega m The mechanical angular speed of the motor; b is a friction coefficient; j is the moment of inertia, and ρ is the parameter perturbation and the total disturbance of the external load.
6. The method for controlling the dead-beat prediction current loop of the permanent magnet synchronous motor according to claim 5, wherein in the first step, the speed loop high-order sliding mode observer is as follows:
Figure FDA00038589525100000314
in the formula (I), the compound is shown in the specification,
Figure FDA00038589525100000315
in order to observe the error in the speed,
Figure FDA00038589525100000316
as an observed value of p,
Figure FDA00038589525100000317
is the observed value of D, which is the derivative of p,
Figure FDA00038589525100000318
is composed of
Figure FDA00038589525100000319
Derivative of (k), k w1 Is a first control parameter, k, of a high-order sliding mode observer of a speed ring w2 Is a second control parameter, k, of a speed-loop high-order sliding-mode observer w3 Is a third control parameter of the speed loop high-order sliding mode observer,
Figure FDA00038589525100000320
for mechanical angular velocity omega of the motor m Is measured in a time-domain manner by a time-domain,
Figure FDA00038589525100000321
is composed of
Figure FDA00038589525100000322
The derivative of (c).
7. The method for controlling the dead-beat prediction current loop of the permanent magnet synchronous motor according to claim 6, wherein the first step further comprises a step of discretizing a high-order sliding mode observer of the speed loop, and specifically comprises the following steps:
Figure FDA0003858952510000041
wherein the content of the first and second substances,
Figure FDA0003858952510000042
is the observed value of the mechanical angular speed of the motor in the k +1 modulation period,
Figure FDA0003858952510000043
is the observed value of the k modulation period of the mechanical angular speed of the motor,
Figure FDA0003858952510000044
is the observed value of p at the k-th modulation period,
Figure FDA0003858952510000045
is the observed value of p at the k +1 modulation period,
Figure FDA0003858952510000046
is the observed value of D at the k-th modulation period,
Figure FDA0003858952510000047
is the observed value of D at the k +1 modulation period.
8. The method for controlling the dead-beat prediction current loop of the permanent magnet synchronous motor according to claim 7, wherein in the third step, the method for improving the sliding mode approach law by using the actual rotating speed and the given value of the rotating speed of the motor comprises the following steps:
defining a slip form surface:
Figure FDA0003858952510000048
wherein e is w Is the error of the rotating speed; c is an integral coefficient of an integral sliding mode surface; s is a slip form surface;
Figure FDA0003858952510000049
in the formula (I), the compound is shown in the specification,
Figure FDA00038589525100000410
setting a rotating speed value;
the improved sliding mode approach law is as follows:
Figure FDA00038589525100000411
in the formula (I), the compound is shown in the specification,
Figure FDA00038589525100000412
the derivative of the slip-form surface s, k s For switching the gain, ε is a variable term ε + (1- ε) e -δ|s| Gain of (k) t For linear gain, δ is an exponential term e -δ|s| And k is s >0,k t More than 0, delta more than 0, more than 0 and less than 1, when the system is far away from the sliding mode surface,
Figure FDA00038589525100000413
a value of greater than k s The convergence speed is accelerated, when the system reaches the sliding mode,
Figure FDA00038589525100000414
value of less than k s And buffeting is limited.
9. The method for controlling the dead-beat predicted current loop of the permanent magnet synchronous motor according to claim 8, wherein in the third step, the method for obtaining the speed control law comprises the following steps:
compensating the observed parameter perturbation and the total disturbance rho of the external load into a control law, wherein the obtained speed control law is as follows:
Figure FDA00038589525100000415
Figure FDA00038589525100000416
the set value of the mechanical angular speed of the motor,
Figure FDA00038589525100000417
is composed of
Figure FDA00038589525100000418
The derivative of (c).
CN202211159454.4A 2022-09-22 2022-09-22 Dead-beat prediction current loop control method for permanent magnet synchronous motor Pending CN115347841A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211159454.4A CN115347841A (en) 2022-09-22 2022-09-22 Dead-beat prediction current loop control method for permanent magnet synchronous motor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211159454.4A CN115347841A (en) 2022-09-22 2022-09-22 Dead-beat prediction current loop control method for permanent magnet synchronous motor

Publications (1)

Publication Number Publication Date
CN115347841A true CN115347841A (en) 2022-11-15

Family

ID=83956494

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211159454.4A Pending CN115347841A (en) 2022-09-22 2022-09-22 Dead-beat prediction current loop control method for permanent magnet synchronous motor

Country Status (1)

Country Link
CN (1) CN115347841A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117811445A (en) * 2024-02-28 2024-04-02 华侨大学 Novel ultra-spiral sliding mode robust load observation method for permanent magnet synchronous motor

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117811445A (en) * 2024-02-28 2024-04-02 华侨大学 Novel ultra-spiral sliding mode robust load observation method for permanent magnet synchronous motor
CN117811445B (en) * 2024-02-28 2024-05-21 华侨大学 Ultra-spiral sliding mode robust load observation method for permanent magnet synchronous motor

Similar Documents

Publication Publication Date Title
Jiang et al. Improved deadbeat predictive current control combined sliding mode strategy for PMSM drive system
WO2021017237A1 (en) Deadbeat control system and method for permanent magnet synchronous motor under low carrier ratio
CN111342720B (en) Permanent magnet synchronous motor self-adaptive continuous sliding mode control method based on torque observation
CN110764418B (en) Active disturbance rejection controller based on limited time convergence extended state observer
CN112422004B (en) Disturbance suppression method for permanent magnet synchronous motor in weak magnetic control mode
CN110190795B (en) Permanent magnet synchronous motor cascade type robust prediction current control method
Kuang et al. Speed control of load torque feedforward compensation based on linear active disturbance rejection for five-phase PMSM
CN111600518B (en) Design method of permanent magnet synchronous current controller based on extended state observer
Liu et al. Improved backstepping control with nonlinear disturbance observer for the speed control of permanent magnet synchronous motor
CN108377117B (en) Permanent magnet synchronous motor composite current control system and method based on predictive control
CN111211717B (en) IPMSM (intelligent power management system) position-sensorless motor closed-loop control method of nonsingular sliding mode structure
CN110995102A (en) Direct torque control method and system for permanent magnet synchronous motor
CN115250085A (en) Complementary terminal sliding mode speed control method, system, equipment and storage medium
CN112953335A (en) Finite time self-adaptive composite control method and system for permanent magnet synchronous motor
CN115347841A (en) Dead-beat prediction current loop control method for permanent magnet synchronous motor
CN110112980B (en) Control method and system for closed-loop iterative learning based on fractional calculus
CN111313773A (en) Permanent magnet synchronous motor parameter identification method based on variable step length LMS algorithm
CN112422014B (en) Permanent magnet synchronous motor rotating speed prediction method based on high-order sliding mode compensation
Lin et al. Adaptive backstepping PI sliding-mode control for interior permanent magnet synchronous motor drive systems
Sun et al. Improved deadbeat predictive current control with extended state observer for dual three-phase PMSMs
CN113364375A (en) Sensorless control method for PMSM (permanent magnet synchronous motor) driving system of variable-structure current regulator
Wang et al. An improved deadbeat predictive current control based on parameter identification for PMSM
CN115102444B (en) Permanent magnet synchronous motor self-adaptive integral sliding mode prediction control method
Yang et al. Sliding-mode observer based rotor resistance updating method for indirect vector controlled induction motor
CN114499312A (en) Model-free current prediction control method and system for double three-phase permanent magnet synchronous motor

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