CN111211718B - Automatic Disturbance Rejection Controller Parameter Automatic Adjustment System for Permanent Magnet Synchronous Motor Vector Control - Google Patents

Automatic Disturbance Rejection Controller Parameter Automatic Adjustment System for Permanent Magnet Synchronous Motor Vector Control Download PDF

Info

Publication number
CN111211718B
CN111211718B CN202010038985.2A CN202010038985A CN111211718B CN 111211718 B CN111211718 B CN 111211718B CN 202010038985 A CN202010038985 A CN 202010038985A CN 111211718 B CN111211718 B CN 111211718B
Authority
CN
China
Prior art keywords
controller
parameter
disturbance rejection
lion
axis current
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
CN202010038985.2A
Other languages
Chinese (zh)
Other versions
CN111211718A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202010038985.2A priority Critical patent/CN111211718B/en
Publication of CN111211718A publication Critical patent/CN111211718A/en
Application granted granted Critical
Publication of CN111211718B publication Critical patent/CN111211718B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Molecular Biology (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Power Engineering (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种用于永磁同步电机矢量控制的自抗扰控制器参数自动调节系统,包括:自抗扰控制器参数选取指标计算模块(COM)与自抗扰控制器参数寻优模块(OPT)。所述的COM负责对转速、d轴电流、q轴电流的误差以及它们的自抗扰控制器的输出信号进行采样并对参数选取指标H进行计算;所述的OPT以H为目标函数对3个控制器的参数进行寻优,最后输出优化后的转速、d轴电流、q轴电流自抗扰控制器的参数组给控制器。本发明基于狮群算法这种性能卓越的群智能算法,可以实现更加高效的控制器参数自整定;用自整定代替传统的人工调节,从而大大降低了操作人员的工作量。

Figure 202010038985

The invention discloses a parameter automatic adjustment system of an active disturbance rejection controller for vector control of a permanent magnet synchronous motor, comprising: an active disturbance rejection controller parameter selection index calculation module (COM) and an active disturbance rejection controller parameter optimization module (OPT). The COM is responsible for sampling the rotational speed, the error of the d-axis current, the q-axis current and the output signals of their active disturbance rejection controllers and calculates the parameter selection index H; the OPT takes H as the objective function to 3 The parameters of each controller are optimized, and finally the optimized parameters of the speed, d-axis current, and q-axis current of the ADR controller are output to the controller. Based on the swarm intelligence algorithm with excellent performance, the present invention can realize more efficient self-tuning of controller parameters; the traditional manual adjustment is replaced by self-tuning, thereby greatly reducing the workload of operators.

Figure 202010038985

Description

Automatic parameter adjusting system of active disturbance rejection controller for vector control of permanent magnet synchronous motor
Technical Field
The invention belongs to the field of motor control, and particularly relates to an auto-disturbance-rejection controller parameter automatic adjusting system for vector control of a permanent magnet synchronous motor.
Background
A Permanent Magnet Synchronous Motor (PMSM) is a Synchronous motor that generates a Synchronous rotating magnetic field by excitation of Permanent magnets, the Permanent magnets serve as a rotor to generate a rotating magnetic field, and a three-phase stator winding reacts through an armature under the action of the rotating magnetic field to induce three-phase symmetrical current.
Modeling the PMSM, assuming:
(1) no core saturation;
(2) no winding leakage inductance;
(3) no hysteresis and eddy current loss;
(4) no higher harmonics of the magnetic field;
(5) the motor rotor is not provided with a damping winding;
(6) the electrical conductivity of the permanent magnet is zero.
Therefore, a mathematical model of the PMSM in a three-phase coordinate system can be established, wherein the mathematical model comprises a voltage equation, a flux linkage equation, a torque equation and a motion equation.
The voltage equation is
Figure BDA0002367053820000011
The flux linkage equation is
Figure BDA0002367053820000021
The torque equation is
Figure BDA0002367053820000022
Equation of motion of
Figure BDA0002367053820000023
The three-phase windings of the PMSM are assumed to be in a neutral star connection, and if the three-phase windings are in a triangular connection, the three-phase windings can be equivalently in a star connection. The following constraints can be obtained
Figure BDA0002367053820000024
It can be seen that under the above conditions, only two phases of the three-phase variables are independent, so the three-phase mathematical model is not the simplest description of the physical object, and therefore it is necessary to replace the two-phase model.
And (3) after the model equation under the three-phase static coordinate system is subjected to Clarke and Park coordinate transformation, a mathematical model of the permanent magnet synchronous motor under a d-q coordinate system can be obtained. As the control system of the PMSM mostly adopts a vector control system, the basic idea of the vector control system is to orient according to the rotor flux linkage, namely, the d axis is coincided with the rotor flux linkage vector, so that a mathematical model of the PMSM under a d-q coordinate system is obtained.
The voltage equation is
Figure BDA0002367053820000025
The flux linkage equation is
Figure BDA0002367053820000026
The torque equation is
Te=pnfiq+(Ld-Lq)idiq] (8)
Equation of motion of
Figure BDA0002367053820000027
To simplify the control scheme, i is often useddControl strategy is 0. If the d-axis current is 0 at steady state, then the equation at PMSM steady state is
Figure BDA0002367053820000031
It can be seen that the torque of the motor is only related to the magnitude of the q-axis current, and the electromagnetic torque only has a permanent magnet torque component and does not contain a reluctance torque component. The armature reaction has no direct-axis demagnetization component, and the motor performance can not be deteriorated due to demagnetization. If the voltage component of d-axis is not considered, from uqThe motor is equivalent to a separately excited direct current motor, the d-axis winding is equivalent to open circuit, so that complete decoupling between the d-axis and the motor torque and between the d-axis and the stator winding is realized, a model of the motor is simplified, and control is facilitated.
However, this control strategy has certain disadvantages, such as that when the load increases, the leakage inductance voltage drop increases, and thus the power factor decreases; the lack of reluctance torque can also make the motor unable to reach the theoretical maximum torque, and at the same time unable to achieve flux weakening speed regulation, and the speed regulation range is very limited.
The control strategy is very suitable for a PMSM (permanent magnet synchronous motor) which is small in capacity and narrow in speed regulation range.
idThe PMSM state equation under the control strategy of 0 is
Figure BDA0002367053820000032
Here, the current i is needed for the rotation speed omega and the d axisdAnd q-axis current iqThese 3 amounts were controlled. PID controllers are commonly used in general industrial applications. The PID control technique has been proposed formally in 1936 by colander (a. Callender) and Stevenson (a. Stevenson) in the uk, and has a very important position in the automatic control technique. The overall control structure of the PMSM vector control system based on the PID controller is shown in fig. 1.
With the rapid development of scientific technology, the requirements on the controller are higher and higher, and the following defects of the traditional PID control technology are gradually revealed:
(1) the error extraction method may cause large initial error and easily cause overshoot;
(2) the differential part is easy to enlarge interference, and pollutes signal distortion to generate huge errors;
(3) the integral link in the classical PID control has obvious effect on inhibiting constant disturbance, but when no disturbance exists, the dynamic characteristic of the system is poor (the closed-loop system has slow reaction and is easy to generate side effects of oscillation and saturation of control quantity), and the inhibition capability of the integral link is not obvious for the disturbance changing along with time;
(4) classical PID control is a weighted combination (i.e., a simple linear combination) of the present (P), past (I), and future (D) of errors, which is not reasonable for many complex controlled objects.
In practical application, the current loop and the rotating speed loop of the PMSM control system based on the PID controller have the problems in many cases, and the PID controller is difficult to apply to occasions with higher requirements on the performance of the motor (such as hypersonic aircrafts, tank fire control systems and the like).
To overcome the defects of the PID Control technology, mr. hangul has formally proposed Auto-Disturbance Rejection Control (ADRC) in 2002. The active disturbance rejection control uses the 'error feedback to eliminate errors' in classical control, and uses the observation of system state quantity in modern control. By combining the essences of the two control methods, the active disturbance rejection controller can control various objects quickly, accurately, without overshoot and stably.
ADRC is used in many fields, such as hypersonic aircrafts, magnetic levitation, superconducting particle accelerators, large radio telescopes, tank fire control systems, nuclear power station cooling systems and the like, and industrial robots, servo motor drivers, unmanned plane flight control systems, sweeping robots, even website search engines and the like.
The active disturbance rejection controller mainly comprises a Tracking Differentiator (TD), a nonlinear state error feedback control law (NLSEF) and an Extended State Observer (ESO).
TD: the transition process of the closed-loop system is arranged through the TD, and the overshoot of the system is reduced. Because discontinuous and random noise input signals often appear in actual engineering, a complex input tracking signal can be transited into a stable continuous input signal through extracting a signal and arranging a reasonable transition process, so that the control quality of the system is improved.
NLSEF: the PID controller calculates the error signal by using a linear superposition method with the simplest form, and the control efficiency and effect are difficult to meet the requirements of people, so the active disturbance rejection controller performs nonlinear calculation on the error signal through NLSEF. The control performance of the linear superposition method for processing the error signal is not as good as that of a nonlinear algorithm of NLSEF through mathematical derivation and simulation.
ESO: the observer is difficult to directly observe the state of the system, so that the internal state and interference of the system are observed in real time through the ESO, and the controllable and observable performance of the system is greatly improved. The interference is external interference suffered by the system and interference factors in the system, and the internal state of the system and the state of an unmodeled part in the system are directly estimated from the output quantity of the system only by setting parameters in a proper extended state observer, so that the interference is estimated and compensated in real time.
In summary, the improved active disturbance rejection controllers in the three parts can improve the response speed, reduce overshoot, and control a nonlinear system which is difficult to control by other controllers with high precision, thereby improving the application capability of the nonlinear system in practical engineering.
When the controlled object is a first-order system, the state equation is assumed to be
Figure BDA0002367053820000041
Where x is a state variable, f (x) is a polynomial associated with the state variable x, u is a control quantity, and w is a disturbance quantity. Then the expressions for TD, NLSEF, and ESO of the active disturbance rejection controller are
Figure BDA0002367053820000051
Figure BDA0002367053820000052
Figure BDA0002367053820000053
Wherein
Figure BDA0002367053820000054
r is a velocity factor, determining the tracking velocity. The larger r, the faster the tracking speed, but at the same time, the overshoot is increased.
k1For feedback gain, increase k1The response speed of the system can be increased, but if the value is too large, the system can oscillate or even be unstable.
β1And beta2To be transportedAnd an error correction factor is generated, and the dynamic performance of the system is greatly influenced. Estimation of state variables is primarily governed by beta1The estimation of the system disturbance is mainly influenced by beta2The influence of (c). Beta is a1And beta2The larger the state variable and the faster the estimate of the system disturbance will converge; however, if the value is too large, the output of the ESO will generate oscillation divergence and generate a high frequency noise signal.
Alpha is a nonlinear factor, and the smaller alpha, the stronger the nonlinearity of the fal function. Alpha is alpha1、α2、α3Respectively TD nonlinear factor, NLSEF nonlinear factor and ESO nonlinear factor.
Delta is a filtering factor, and increasing delta can make the filtering effect better, but also increases the delay of tracking. Delta1、δ2、δ3TD filter factor, NLSEF filter factor, ESO filter factor.
Thus, when the controlled object is a first-order system, the active disturbance rejection controller contains r, α1,δ1,k1,α2,δ2,β1,β2,α3And delta3These 10 parameters to be adjusted (parameter b is determined by the controlled object and no adjustment is needed). The structure of the active disturbance rejection controller is shown in fig. 2.
And designing an active disturbance rejection controller for a rotating speed loop, a d-axis current loop and a q-axis current loop of the PMSM.
The rotating speed loop state equation of PMSM is
Figure BDA0002367053820000055
It is written in the form shown in formula (12)
Figure BDA0002367053820000056
Wherein
Figure BDA0002367053820000061
Thus, the expressions for the speed loop TD, NLSEF and ESO can be written as
Figure BDA0002367053820000062
Figure BDA0002367053820000063
Figure BDA0002367053820000064
Wherein r1, α 11, δ 11, k11, α 21, δ 21, β 11, β 21, α 31, and δ 31 are the speed factor, TD nonlinear factor, TD filter factor, feedback gain, NLSEF nonlinear factor, NLSEF filter factor, output error correction factor, ESO nonlinear factor, and ESO filter factor of the rotational speed active disturbance rejection controller, respectively.
The d-axis current loop state equation of PMSM is
Figure BDA0002367053820000065
It is written in the form shown in formula (12)
Figure BDA0002367053820000066
Wherein
Figure BDA0002367053820000067
Thus the expressions for d-axis current loop TD, NLSEF, and ESO can be written as
Figure BDA0002367053820000071
Figure BDA0002367053820000072
Figure BDA0002367053820000073
Wherein r2, α 12, δ 12, k12, α 22, δ 22, β 12, β 22, α 32, δ 32 are the speed factor, TD nonlinear factor, TD filter factor, feedback gain, NLSEF nonlinear factor, NLSEF filter factor, output error correction factor, ESO nonlinear factor, ESO filter factor of the d-axis current active disturbance rejection controller, respectively.
The q-axis current loop state equation of PMSM is
Figure BDA0002367053820000074
It is written in the form shown in formula (12)
Figure BDA0002367053820000075
Wherein
Figure BDA0002367053820000076
Thus the expressions for the q-axis current loop TD, NLSEF and ESO can be written as
Figure BDA0002367053820000077
Figure BDA0002367053820000078
Figure BDA0002367053820000081
Wherein r3, α 13, δ 13, k13, α 23, δ 23, β 13, β 23, α 33, and δ 33 are the speed factor, TD nonlinear factor, TD filter factor, feedback gain, NLSEF nonlinear factor, NLSEF filter factor, output error correction factor, ESO nonlinear factor, and ESO filter factor of the q-axis current active disturbance rejection controller, respectively.
The structure of the PMSM vector control system based on the active disturbance rejection control is shown in fig. 3.
Although the performance of the auto-disturbance-rejection controller has proved to be superior to that of the conventional PID controller in many experiments, its parameters to be tuned are excessive. There are 3 total active disturbance rejection controllers, each controller has 10 parameters to be adjusted, and then there are 30 parameters to be adjusted for 3 controllers, if the manual adjustment mode is adopted, the workload is very large, and the manual adjustment is difficult to achieve the optimal control effect.
Disclosure of Invention
Aiming at the technical defects in the prior art, the invention provides an auto-disturbance-rejection controller parameter automatic adjusting device for vector control of a permanent magnet synchronous motor. The device can realize more efficient parameter self-tuning of the controller based on the superior group intelligent algorithm such as the lion group algorithm, and the workload of operators is greatly reduced.
The invention relates to an auto-disturbance rejection controller parameter automatic regulating system for permanent magnet synchronous motor vector control, which comprises an auto-disturbance rejection controller parameter selection index calculating module (COM) and an auto-disturbance rejection controller parameter optimizing module (OPT);
the parameter selection index calculation module (COM) of the active disturbance rejection controller obtains a rotating speed error (e)1) Rotational speed auto-disturbance rejection controller (ADRC)1) Is output signal (u)1) D-axis current error (e)2) D-axis current auto-disturbance rejection controller (ADRC)2) Is output signal (u)2) Q-axis current error (e)3) Q-axis current auto-disturbance rejection controller (ADRC)3) Is output signal (u)3) (ii) a Output active disturbance rejection controlSelecting an index H from the device parameters; the working state of an auto-disturbance rejection controller parameter selection index calculation module (COM) is controlled by a driving signal;
the parameter optimizing module (OPT) of the active disturbance rejection controller obtains an active disturbance rejection controller parameter selection index H, and outputs 3 groups of optimized controller parameter groups, namely a parameter group (p) of the rotating speed active disturbance rejection controller1) D-axis current auto-disturbance rejection controller parameter set (p)2) Q-axis current auto-disturbance-rejection controller parameter set (p)3)。
As a preferred scheme of the present invention, the process of the auto-disturbance rejection controller parameter selection index calculation module (COM) for obtaining the auto-disturbance rejection controller parameter selection index H specifically includes:
setting the parameter selection index of the rotation speed active disturbance rejection controller as H1By selecting e in real time1、u1And a rise time t of the rotational speedu1To H1And (6) performing calculation. H1Is expressed as
Figure BDA0002367053820000091
Wherein, TsFor a sampling period, TaFor the total sampling time of the controller state in one iteration of the OPT, int is a floor function, w11、w21、w31Is the weight; e.g. of the type1(k) Is e1In discrete form u1(k) Is u1In discrete form.
Setting the parameter selection index of the d-axis current active disturbance rejection controller as H2By selecting e in real time2、u2And d-axis current rise time tu2To H2Is calculated, H2Is expressed as
Figure BDA0002367053820000092
Wherein, w12、w22、w32As a weight value, e2(k) Is e2In discrete form u2(k) Is u2In discrete form.
Setting the parameter selection index of the q-axis current active disturbance rejection controller as H3By selecting e in real time3、u3And q-axis current rise time tu3To H3Is calculated, H3Is expressed as
Figure BDA0002367053820000093
Wherein, w13、w23、w33As a weight value, e3(k) Is e3In discrete form u3(k) Is u3In discrete form.
COM then passes through pair H1、H2、H3Carrying out weighted summation to obtain a comprehensive evaluation index H capable of comprehensively measuring the performances of 3 active disturbance rejection controllers, wherein the expression of H is as follows
Figure BDA0002367053820000094
Wherein c is1、c2、c3Is a weight; the smaller H, the better the controller performance.
In a preferred embodiment of the present invention, the rotation speed rise time t isu1The time required for the first time the rotational speed ω reaches the steady state value from zero; d-axis current rise time tu2Is d-axis current idThe time required to reach a steady state value for the first time from zero; q-axis current rise time tu3Is q-axis current iqThe time required to reach the steady state value for the first time from time zero.
As a preferred embodiment of the present invention, the process of the auto disturbance rejection controller parameter optimization module (OPT) for obtaining 3 sets of optimized controller parameter sets specifically includes:
OPT takes H as an objective function to optimize 30 parameters of 3 controllers, and the controller parameter optimization problem is equivalent to a function maximum solving problem:
Figure BDA0002367053820000101
s.t.pimin≤pi≤pimax(i=1,2,3)
wherein p isimaxFor the upper limit, p, of the 3 ADRC parameter setsiminThe lower limit of the 3 active disturbance rejection controller parameter sets, i is 1,2 and 3; are respectively expressed as
pi=(ri1i1i,k1i2i2i1i2i3i3i)
pimax=(rimax1imax1imax,k1imax2imax2imax1imax2imax3imax3imax)
pimin=(rimin1imin1imin,k1imin2imin2imin1imin2imin3imin3imin)
r is a velocity factor, k1For feedback gain, beta1And beta2For outputting the error correction factor, alpha is a non-linear factor, alpha1、α2、α3TD nonlinear factor, NLSEF nonlinear factor and ESO nonlinear factor; delta is the filter factor, delta1、δ2、δ3TD filter factor, NLSEF filter factor, ESO filter factor.
As a preferred embodiment of the present invention, the OPT adjusts the controller parameters based on the lion group algorithm to combine 3 sets of parameters (p)1,p2,p3) As the position information of the lions, the fitness function for representing the degree of goodness of the position of each lion is taken as
Figure BDA0002367053820000102
Wherein epsilon0Is a constant greater than 0 such that the denominator of the fitness function is not 0. The larger f, the better the position, i.e. parameter, of each lion;
stopping optimization after the iteration end condition is reached, and outputting the final controller parameter group p1、p2、p3To the controller.
As a preferred scheme of the invention, the OPT is based on a lion group algorithm, so that the position x of each lioniIs composed of
xi=(p1,p2,p3)
The steps of adjusting the controller parameters are as follows:
step 9, initializing the position x of the lion in the lion groupiThe number N, the maximum iteration number T, the dimension space D and the scale factor beta of the adult lion in the lion group are calculated;
step 10, calculating the number of the lion king, the adult female lion and the young lion, setting the individual historical optimal position as the current position of each lion, and setting the initial group optimal position as the lion king position;
step 11, updating the position of the lion king and calculating a fitness value;
step 12, updating the position of the female lion;
step 13, updating the position of the young lion;
step 14, calculating a fitness value according to the position of the lion, updating the historical optimal position of the lion and the historical optimal position of the lion group, judging whether the algorithm meets an end condition, and turning to Step 8 if the algorithm meets the end condition, or turning to Step 7 if the algorithm does not meet the end condition;
step 15, reordering every certain iteration number, determining the positions of the lion king, the adult female lion and the young lion, and turning to Step 3;
and Step 16, outputting the position of the lion king, namely the optimal solution of the solved problem, and ending the algorithm.
Stopping optimization after the iteration end condition is reached, and outputting the final controller parameter group p1,p2,p3To the controller, the control system is then ready for use. The structure of the PMSM vector control system with automatic parameter adjustment is shown in fig. 4.
The beneficial technical effects of the invention are as follows:
(1) because the parameters of the active disturbance rejection controller are many, and the number of the active disturbance rejection controller is up to 30 in the present example, the invention can carry out automatic adjustment on the active disturbance rejection controller, thereby replacing the traditional manual adjustment, greatly reducing the workload of debugging personnel, avoiding accidental errors which are possibly generated by the manual adjustment, and being very suitable for being applied to a complex multi-parameter control system such as a PMSM vector control system based on the active disturbance rejection control.
(2) The auto-disturbance rejection controller parameter optimization module (OPT) in the invention is based on the lion group algorithm. Compared with the traditional swarm intelligence algorithms such as an ant colony algorithm, a particle swarm algorithm, an artificial fish swarm algorithm, a mixed frog-leaping algorithm, a firefly optimization algorithm, a wolf colony algorithm, an artificial bee colony algorithm and the like, the lion colony algorithm has higher convergence speed, higher precision and stronger global search capability. Therefore, the OPT based on the algorithm can theoretically greatly improve the parameter optimization efficiency compared with a parameter optimization module based on a traditional algorithm.
Drawings
FIG. 1 is a block diagram of a PMSM vector control system based on a PID controller;
FIG. 2 is a block diagram of an active disturbance rejection controller;
FIG. 3 is a block diagram of an active disturbance rejection controller based PMSM vector control system;
FIG. 4 is a block diagram of a PMSM vector control system with an automatic parameter adjustment system;
fig. 5 is a waveform diagram of a pulse type rotational speed setting signal.
Detailed Description
In order to describe the present invention more specifically, the following detailed description will be made of the technical solutions of the present invention and the related working principles.
As shown in fig. 1, an auto-disturbance rejection controller parameter automatic adjustment system for vector control of a permanent magnet synchronous motor includes an auto-disturbance rejection controller parameter selection index calculation module (COM).
COM obtains the rotation speed error (e)1) Rotational speed auto-disturbance rejection controller (ADRC)1) Is output signal (u)1) D-axis current error (e)2) D-axis current auto-disturbance rejection controller (ADRC)2) Is output signal (u)2) Q-axis current error (e)3) Q-axis current auto-disturbance rejection controller (ADRC)3) Is output signal (u)3) (ii) a Outputting parameters of the active disturbance rejection controller to select an index H; the working state of an auto-disturbance rejection controller parameter selection index calculation module (COM) is controlled by a driving signal;
the COM calculates the obtained signals to obtain a comprehensive evaluation index H which can comprehensively measure the performance of the 3 active disturbance rejection controllers, and the expression of the H is as follows
Figure BDA0002367053820000121
The smaller H, the better the control performance of the controller. The control performance of 3 controllers is represented by an index H, so that the complexity of an evaluation system is greatly reduced.
As shown in fig. 1, an auto-disturbance rejection controller parameter automatic adjustment system for vector control of a permanent magnet synchronous motor includes an auto-disturbance rejection controller parameter optimization module (OPT).
OPT obtains the parameter selection index H of the active disturbance rejection controller, and outputs 3 groups of optimized controller parameter groups, namely the parameter group (p) of the rotating speed active disturbance rejection controller1) D-axis current auto-disturbance rejection controller parameter set (p)2) Q-axis current auto-disturbance-rejection controller parameter set (p)3)。
OPT takes H as an objective function to optimize 30 parameters of 3 controllers, and the controller parameter optimization problem is equivalent to a function maximum solving problem
Figure BDA0002367053820000122
s.t.pimin≤pi≤pimax(i=1,2,3)
Wherein p isi(i-1, 2,3) are a rotation speed active-disturbance-rejection controller, a d-axis current active-disturbance-rejection controller and a q-axis current active-disturbance-rejection controller respectivelyController 3 parameter sets of the active disturbance rejection controller, pimax(i is 1,2,3) is the upper limit, p, of the 3 sets of auto-disturbance-rejection controller parametersimin(i ═ 1,2,3) is the lower limit of the 3 sets of auto-disturbance-rejection controller parameters; are respectively expressed as
pi=(ri1i1i,k1i2i2i1i2i3i3i)
pimax=(rimax1imax1imax,k1imax2imax2imax1imax2imax3imax3imax)
pimin=(rimin1imin1imin,k1imin2imin2imin1imin2imin3imin3imin)
The OPT carries out parameter optimization based on a lion group algorithm. Position x of each lion in the algorithmiIs composed of
xi=(p1,p2,p3)
The steps of adjusting the controller parameters are as follows:
step 17, initializing the position x of the lion in the lion groupiThe number N, the maximum iteration number T, the dimension space D and the scale factor beta of the adult lion in the lion group are calculated;
step 18, calculating the number of the lion king, the adult female lion and the young lion, setting the individual historical optimal position as the current position of each lion, and setting the initial group optimal position as the lion king position;
step 19, updating the position of the lion king and calculating a fitness value;
step 20, updating the position of the female lion;
step 21, updating the position of the young lion;
step 22, calculating a fitness value according to the position of the lion, updating the historical optimal position of the lion and the historical optimal position of the lion group, judging whether the algorithm meets an end condition, and turning to Step 8 if the algorithm meets the end condition, or turning to Step 7 if the algorithm does not meet the end condition;
step 23, reordering every certain iteration number, determining the positions of the lion king, the adult female lion and the young lion, and turning to Step 3;
and Step 24, outputting the position of the lion king, namely the optimal solution of the solved problem, and ending the algorithm.
The OPT carries out parameter optimization based on a lion group algorithm. Compared with the traditional group intelligent algorithm, the lion group algorithm has higher convergence speed, higher precision and stronger global search capability. Therefore, the OPT based on the algorithm can theoretically greatly improve the parameter optimization efficiency compared with a parameter optimization module based on a traditional algorithm.
When the parameters of the active disturbance rejection controller are optimized, the automatic parameter adjusting system is started through the driving signal, and then the parameters of the controller can be automatically optimized by adopting a pulse type rotating speed setting signal. And outputting the final parameters to the controller after the optimization process is finished, and closing the automatic parameter adjusting system through a driving signal.
The pulse type rotational speed setting signal is shown in fig. 5. Wherein t is1Optimizing the time, t, for the parameter2For system recovery time, T is the pulse period. t is t1In a time period, the rotating speed given signal is equivalent to a step signal, the H value can be obtained by calculation after the state of 3 controllers is sampled through COM in the time period, and then the COM outputs the H value to OPT for parameter optimization; t is t2The rotation speed setting signal is set to 0 for a period of time long enough to allow the system to return to approximately the state at time 0. The period of the pulse is T, the OPT is iterated once after each period, the optimization is stopped until the end condition is met, and the final parameters are output to the controller and put into use. The automatic parameter adjusting system realizes automatic adjustment of parameters, reduces the workload of operators, avoids accidental errors which are possibly generated by manual adjustment, and is very suitable for being applied to a complex multi-parameter control system such as a PMSM vector control system based on an active disturbance rejection controller.

Claims (4)

1.一种用于永磁同步电机矢量控制的自抗扰控制器参数自动调节系统,其特征在于包括自抗扰控制器参数选取指标计算模块COM与自抗扰控制器参数寻优模块OPT;1. an automatic disturbance rejection controller parameter automatic adjustment system for permanent magnet synchronous motor vector control, is characterized in that comprising the active disturbance rejection controller parameter selection index calculation module COM and the active disturbance rejection controller parameter optimization module OPT; 所述的自抗扰控制器参数选取指标计算模块COM获取转速误差e1、转速自抗扰控制器ADRC1的输出信号u1、d轴电流误差e2、d轴电流自抗扰控制器ADRC2的输出信号u2、q轴电流误差e3、q轴电流自抗扰控制器ADRC3的输出信号u3;输出自抗扰控制器参数选取指标H;自抗扰控制器参数选取指标计算模块COM的工作状态由驱动信号控制;The parameter selection index calculation module COM of the active disturbance rejection controller obtains the rotational speed error e 1 , the output signal u 1 of the rotational speed active disturbance rejection controller ADRC 1 , the d-axis current error e 2 , and the d-axis current active disturbance rejection controller ADRC Output signal u 2 of 2 , q-axis current error e 3 , output signal u 3 of q-axis current ADRC 3 ; output ADRC parameter selection index H; ADRC parameter selection index calculation The working state of the module COM is controlled by the drive signal; 自抗扰控制器参数选取指标计算模块COM求取自抗扰控制器参数选取指标H的过程具体为:The process of calculating the parameter selection index H of the ADRC controller by the calculation module COM is as follows: 设转速自抗扰控制器参数选取指标为H1,通过实时选取e1、u1以及转速上升时间tu1来对H1进行计算,H1的表达式为Set the parameter selection index of the speed active disturbance rejection controller as H 1 , and calculate H 1 by selecting e 1 , u 1 and the speed rise time t u1 in real time. The expression of H 1 is:
Figure FDA0003019327230000011
Figure FDA0003019327230000011
其中,Ts为采样周期,Ta为OPT一次迭代过程中对控制器状态的总采样时间,int为向下取整函数,w11、w21、w31为权值;e1(k)为e1的离散形式,u1(k)为u1的离散形式;Among them, T s is the sampling period, Ta is the total sampling time of the controller state in one iteration of OPT, int is the round-down function, w 11 , w 21 , and w 31 are weights; e 1 ( k ) is the discrete form of e 1 , and u 1 (k) is the discrete form of u 1 ; 设d轴电流自抗扰控制器参数选取指标为H2,通过实时选取e2、u2以及d轴电流上升时间tu2对H2进行计算,H2的表达式为Assume that the parameter selection index of the d-axis current active disturbance rejection controller is H 2 , and calculate H 2 by selecting e 2 , u 2 and the d-axis current rise time t u2 in real time. The expression of H 2 is:
Figure FDA0003019327230000012
Figure FDA0003019327230000012
其中,w12、w22、w32为权值,e2(k)为e2的离散形式,u2(k)为u2的离散形式;Among them, w 12 , w 22 , and w 32 are weights, e 2 (k) is the discrete form of e 2 , and u 2 (k) is the discrete form of u 2 ; 设q轴电流自抗扰控制器参数选取指标为H3,通过实时选取e3、u3以及q轴电流上升时间tu3来对H3进行计算,H3的表达式为Assume that the parameter selection index of the q-axis current active disturbance rejection controller is H 3 , and calculate H 3 by selecting e 3 , u 3 and the q-axis current rise time t u3 in real time. The expression of H 3 is
Figure FDA0003019327230000013
Figure FDA0003019327230000013
其中,w13、w23、w33为权值,e3(k)为e3的离散形式,u3(k)为u3的离散形式;Among them, w 13 , w 23 , and w 33 are weights, e 3 (k) is the discrete form of e 3 , and u 3 (k) is the discrete form of u 3 ; 然后COM通过对H1、H2、H3进行加权求和得到可以综合衡量3个自抗扰控制器性能的综合评价指标H,H的表达式如下Then COM obtains a comprehensive evaluation index H that can comprehensively measure the performance of the three ADRC controllers by weighted summation of H 1 , H 2 , and H 3 , and the expression of H is as follows
Figure FDA0003019327230000014
Figure FDA0003019327230000014
其中c1、c2、c3为权重;H越小,控制器性能越好;where c 1 , c 2 , and c 3 are weights; the smaller H is, the better the performance of the controller; 所述的自抗扰控制器参数寻优模块OPT获取自抗扰控制器参数选取指标H,输出3组优化后的控制器参数组,即转速自抗扰控制器的参数组p1、d轴电流自抗扰控制器的参数组p2、q轴电流自抗扰控制器的参数组p3The ADRC parameter optimization module OPT obtains the parameter selection index H of the ADRC controller, and outputs 3 groups of optimized controller parameter groups, namely the parameter groups p 1 and d axes of the rotational speed ADRC controller. the parameter group p 2 of the current active disturbance rejection controller, the parameter group p 3 of the q-axis current active disturbance rejection controller ; 自抗扰控制器参数寻优模块OPT求取3组优化后的控制器参数组的过程具体为:The process of obtaining 3 groups of optimized controller parameter groups by the OPT controller parameter optimization module is as follows: OPT以H为目标函数对3个控制器的30个参数进行寻优,控制器参数优化问题等效为一个函数最值求解问题:OPT uses H as the objective function to optimize the 30 parameters of the three controllers. The controller parameter optimization problem is equivalent to a function maximum value solution problem:
Figure FDA0003019327230000021
Figure FDA0003019327230000021
s.t.pimin≤pi≤pimax(i=1,2,3)stp imin ≤ p i ≤ p imax (i=1,2,3) 其中,pimax为3个自抗扰控制器参数组的上限、pimin为3个自抗扰控制器参数组的下限,i=1、2、3,表达式分别为Among them, p imax is the upper limit of the three ADRC parameter groups, p imin is the lower limit of the three ADRC parameter groups, i=1, 2, 3, the expressions are pi=(ri1i1i,k1i2i2i1i2i3i3i)p i =(ri ,α 1i1i ,k 1i2i2i , β 1i , β 2i3i3i ) pimax=(rimax1imax1imax,k1imax2imax2imax1imax2imax3imax3imax)p imax =(r imax1imax1imax ,k 1imax2imax2imax1imax2imax3imax3imax ) pimin=(rimin1imin1imin,k1imin2imin2imin1imin2imin3imin3imin)p imin = ( rimin1imin1imin ,k 1imin2imin2imin1imin2imin3imin3imin ) r为速度因子,k1为反馈增益,β1与β2为输出误差校正因子,α为非线性因子,α1、α2、α3分别为TD非线性因子、NLSEF非线性因子、ESO非线性因子;δ为滤波因子,δ1、δ2、δ3分别为TD滤波因子、NLSEF滤波因子、ESO滤波因子。r is the speed factor, k 1 is the feedback gain, β 1 and β 2 are the output error correction factors, α is the nonlinear factor, α 1 , α 2 , and α 3 are the TD nonlinear factor, NLSEF nonlinear factor, and ESO nonlinear factor, respectively. Linear factor; δ is the filter factor, δ 1 , δ 2 , and δ 3 are the TD filter factor, the NLSEF filter factor, and the ESO filter factor, respectively.
2.根据权利要求1所述的用于永磁同步电机矢量控制的自抗扰控制器参数自动调节系统,其特征在于所述的转速上升时间tu1为转速ω从零时刻第一次到达稳态值所需的时间;d轴电流上升时间tu2为d轴电流id从零时刻第一次到达稳态值所需的时间;q轴电流上升时间tu3为q轴电流iq从零时刻第一次到达稳态值所需的时间。2. The automatic disturbance rejection controller parameter automatic adjustment system for permanent magnet synchronous motor vector control according to claim 1, characterized in that the speed rise time t u1 is that the speed ω reaches the steady state for the first time from zero time. The time required for the state value; the d-axis current rise time t u2 is the time required for the d -axis current id to reach the steady state value for the first time from zero time; the q-axis current rise time t u3 is the q-axis current i q from zero The time it takes for a moment to reach its steady-state value for the first time. 3.根据权利要求2所述的用于永磁同步电机矢量控制的自抗扰控制器参数自动调节系统,其特征在于OPT基于狮群算法对控制器参数进行调节,以3组参数的组合(p1,p2,p3)作为狮子的位置信息,表征每个狮子位置好坏程度的适应度函数取为3. the automatic disturbance rejection controller parameter automatic adjustment system that is used for permanent magnet synchronous motor vector control according to claim 2, it is characterized in that OPT is based on lion group algorithm to adjust the controller parameter, with the combination of 3 groups of parameters ( p 1 , p 2 , p 3 ) are used as the location information of the lions, and the fitness function that characterizes the position of each lion is taken as
Figure FDA0003019327230000022
Figure FDA0003019327230000022
其中,ε0为一个大于0的常数,使得适应度函数的分母不为0,f越大,每个狮子为位置即参数越好;Among them, ε 0 is a constant greater than 0, so that the denominator of the fitness function is not 0, and the larger the f is, the better the parameter is for each lion; 达到迭代结束条件后,停止优化,输出最终的控制器参数组p1、p2、p3给控制器。After reaching the iteration end condition, stop the optimization, and output the final controller parameter set p 1 , p 2 , p 3 to the controller.
4.根据权利要求3所述的用于永磁同步电机矢量控制的自抗扰控制器参数自动调节系统,其特征在于所述的OPT是基于狮群算法的,则每个狮子的位置xi4. the automatic disturbance rejection controller parameter automatic adjustment system for permanent magnet synchronous motor vector control according to claim 3, is characterized in that described OPT is based on lion group algorithm, then the position x i of each lion for xi=(p1,p2,p3)x i =(p 1 ,p 2 ,p 3 ) 对控制器参数进行调节的步骤如下:The steps to adjust the controller parameters are as follows: Step 1:初始化狮群中狮子的位置xi以及数目N,最大迭代次数T,维度空间D,成年狮占狮群的比例因子β;Step 1: Initialize the position x i and the number N of lions in the lion group, the maximum number of iterations T, the dimension space D, and the proportional factor β of adult lions in the lion group; Step 2:计算狮王、成年母狮以及幼狮的个数,将个体历史最优位置设置为各狮子的当前位置,初始群体最优位置设置为狮王位置;Step 2: Calculate the number of lion kings, adult lionesses and cubs, set the individual historical optimal position as the current position of each lion, and the initial group optimal position as the lion king position; Step 3:更新狮王的位置,并计算适应度值;Step 3: Update the position of the lion king and calculate the fitness value; Step 4:更新母狮的位置;Step 4: Update the position of the lioness; Step 5:更新幼狮的位置;Step 5: Update the position of the lion cub; Step 6:根据狮子的位置计算适应度值,更新自身历史最优位置以及狮群历史最优位置,判断算法是否满足结束条件,满足转Step 8,否则转Step 7;Step 6: Calculate the fitness value according to the position of the lion, update its own historical optimal position and the historical optimal position of the lion group, and judge whether the algorithm satisfies the end condition, and then go to Step 8, otherwise, go to Step 7; Step 7:每隔一定迭代次数,重新排序,确定狮王、成年母狮以及幼狮的位置,转Step3;Step 7: Re-order every certain number of iterations to determine the positions of the lion king, adult lioness and cub, and go to Step 3; Step 8:输出狮王的位置,即所求问题的最优解,算法结束。Step 8: Output the position of the lion king, that is, the optimal solution to the problem, and the algorithm ends.
CN202010038985.2A 2020-01-14 2020-01-14 Automatic Disturbance Rejection Controller Parameter Automatic Adjustment System for Permanent Magnet Synchronous Motor Vector Control Active CN111211718B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010038985.2A CN111211718B (en) 2020-01-14 2020-01-14 Automatic Disturbance Rejection Controller Parameter Automatic Adjustment System for Permanent Magnet Synchronous Motor Vector Control

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010038985.2A CN111211718B (en) 2020-01-14 2020-01-14 Automatic Disturbance Rejection Controller Parameter Automatic Adjustment System for Permanent Magnet Synchronous Motor Vector Control

Publications (2)

Publication Number Publication Date
CN111211718A CN111211718A (en) 2020-05-29
CN111211718B true CN111211718B (en) 2021-06-08

Family

ID=70785790

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010038985.2A Active CN111211718B (en) 2020-01-14 2020-01-14 Automatic Disturbance Rejection Controller Parameter Automatic Adjustment System for Permanent Magnet Synchronous Motor Vector Control

Country Status (1)

Country Link
CN (1) CN111211718B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114879502B (en) * 2022-05-23 2023-06-30 中国科学院光电技术研究所 Parameter self-tuning method for position loop active disturbance rejection controller

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9160260B1 (en) * 2014-06-16 2015-10-13 Moog Inc. Adaptive actuator control system
CN106647283A (en) * 2017-01-23 2017-05-10 无锡信捷电气股份有限公司 Auto-disturbance rejection position servo system optimization design method based on improved CPSO
CN107800343A (en) * 2017-11-24 2018-03-13 西安科技大学 The design method of asynchronous machine automatic disturbance rejection controller
CN110095981A (en) * 2019-04-02 2019-08-06 南京交通职业技术学院 A kind of setting method, device and the electronic equipment of automatic disturbance rejection controller parameter
CN110165962A (en) * 2019-05-15 2019-08-23 浙江工业大学 A kind of direct-drive permanent magnetism synchronous wind generating system and its full Auto-disturbance-rejection Control
CN110635738A (en) * 2019-10-12 2019-12-31 东南大学 A real-time identification method for permanent magnet synchronous motor stator resistance and motor temperature

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9160260B1 (en) * 2014-06-16 2015-10-13 Moog Inc. Adaptive actuator control system
CN106647283A (en) * 2017-01-23 2017-05-10 无锡信捷电气股份有限公司 Auto-disturbance rejection position servo system optimization design method based on improved CPSO
CN107800343A (en) * 2017-11-24 2018-03-13 西安科技大学 The design method of asynchronous machine automatic disturbance rejection controller
CN110095981A (en) * 2019-04-02 2019-08-06 南京交通职业技术学院 A kind of setting method, device and the electronic equipment of automatic disturbance rejection controller parameter
CN110165962A (en) * 2019-05-15 2019-08-23 浙江工业大学 A kind of direct-drive permanent magnetism synchronous wind generating system and its full Auto-disturbance-rejection Control
CN110635738A (en) * 2019-10-12 2019-12-31 东南大学 A real-time identification method for permanent magnet synchronous motor stator resistance and motor temperature

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种群体智能算法—狮群算法;刘生建等;《模式识别与人工智能》;20180515;第31卷(第5期);第431-440页 *
基于人群搜索算法的自抗扰控制器参数优化;周志刚等;《热力发电》;20180910;第47卷(第9期);第1-9页 *

Also Published As

Publication number Publication date
CN111211718A (en) 2020-05-29

Similar Documents

Publication Publication Date Title
CN110429881B (en) A kind of active disturbance rejection control method of permanent magnet synchronous motor
CN103401501A (en) Permanent magnet synchronous motor (PMSM) servo system control method based on fuzzy and active disturbance rejection control
CN107800343A (en) The design method of asynchronous machine automatic disturbance rejection controller
Zhou et al. An improved MTPA control based on amplitude-adjustable square wave injection
CN110492804A (en) A kind of permanent magnet synchronous motor Second Order Sliding Mode Control method based on novel disturbance observer
CN110707981A (en) Permanent magnet synchronous motor speed controller based on novel extended state observer
CN110995102A (en) Direct torque control method and system for permanent magnet synchronous motor
CN113556067A (en) A Disturbance Suppression Method for Low-speed Direct Drive Motor Based on Sliding Mode and Disturbance Compensation
CN110011582A (en) A vector control method for permanent magnet synchronous motor
Li et al. Improved active disturbance rejection control of permanent-magnet synchronous motor based on BP neural network
Rong et al. A new PMSM speed modulation system with sliding mode based on active-disturbance-rejection control
Hadi et al. Speed control of an SPMSM using a tracking differentiator-PID controller scheme with a genetic algorithm
Zhang et al. An enhanced nonlinear ADRC speed control method for electric propulsion system: Modeling, analysis, and validation
WO2023237018A1 (en) Nonlinear predictive position control method suitable for dual-axis permanent magnet servo system
Guo et al. Research on a new adaptive integral sliding mode controller based on a small BLDC
CN110492814A (en) The method of particle swarm algorithm optimization synovial membrane structure changes permanent magnet synchronous motor control parameter
Li et al. Overview of sliding mode control technology for permanent magnet synchronous motor system
Du et al. Research on a permanent magnet synchronous motor sensorless anti-disturbance control strategy based on an improved sliding mode observer
CN111211718B (en) Automatic Disturbance Rejection Controller Parameter Automatic Adjustment System for Permanent Magnet Synchronous Motor Vector Control
CN118659694A (en) A speed compensation control method for permanent magnet synchronous servo system
Zhao et al. Disturbance rejection enhancement of vector controlled PMSM using second-order nonlinear ADRC
CN114448308B (en) Variable gain sliding mode control method of permanent magnet synchronous motor speed regulation system
Li et al. Improved Non-Singular Fast Terminal Sensor-Less Sliding Mode Control of IPMSM Considering External Disturbance and Parameter Perturbation.
Zhao et al. Design of MRAC and Modified MRAC for the Turntable
Sonkriwal et al. Analysis of Sliding Mode, FUZZY-PI and PI Control Strategies for Permanent Magnetic Synchronous Motor Drive

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