CN107070315B - Multi-switch reluctance motor rotating speed synchronous control device and control method - Google Patents

Multi-switch reluctance motor rotating speed synchronous control device and control method Download PDF

Info

Publication number
CN107070315B
CN107070315B CN201710398938.7A CN201710398938A CN107070315B CN 107070315 B CN107070315 B CN 107070315B CN 201710398938 A CN201710398938 A CN 201710398938A CN 107070315 B CN107070315 B CN 107070315B
Authority
CN
China
Prior art keywords
motor
torque
virtual
module
drosophila
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
CN201710398938.7A
Other languages
Chinese (zh)
Other versions
CN107070315A (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.)
Hunan University of Science and Technology
Original Assignee
Hunan University of Science and Technology
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 Hunan University of Science and Technology filed Critical Hunan University of Science and Technology
Priority to CN201710398938.7A priority Critical patent/CN107070315B/en
Publication of CN107070315A publication Critical patent/CN107070315A/en
Application granted granted Critical
Publication of CN107070315B publication Critical patent/CN107070315B/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
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/08Arrangements for controlling the speed or torque of a single motor
    • 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/08Reluctance motors

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

The invention discloses a synchronous control method for the rotating speed of a multi-switch reluctance motor, which comprises the following steps: calculating the rotation speed of the virtual main shaft and taking the virtual main shaft as a given rotation speed of each motor; calculating the weight average rotating speed of each motor and taking the weight average rotating speed as the reference rotating speed of each motor torque compensation module; calculating the reference torque of each motor, and comparing the reference torque with the actual torque of the motor to obtain corresponding torque deviation; calculating corresponding compensation torque of each motor, and summing the compensation torque and the torque deviation to obtain torque adjustment quantity of each motor; comparing the actual flux linkage of the motor with a given flux linkage of the system to obtain flux linkage deviation; the motors are controlled according to the torque adjustment quantity and the flux linkage deviation, so that the actual rotating speeds of the motors can accurately track the given rotating speeds of the system, and the aim of synchronous operation of the rotating speeds of multiple motors is fulfilled. The invention also discloses a device for synchronously controlling the rotating speeds of the multi-switch reluctance motors.

Description

Multi-switch reluctance motor rotating speed synchronous control device and control method
Technical Field
The invention relates to the field of control of switched reluctance motors, in particular to a synchronous control device and a synchronous control method for rotating speeds of multiple switched reluctance motors.
Background
The switched reluctance motor has been rapidly popularized and applied in recent years because of its series advantages of simple and firm structure, small starting current, large starting torque, high efficiency, and the like. However, in some complex multi-motor drive systems, synchronous operation of the multiple motors is involved, and synchronous control of the rotational speeds of the multiple motors is required. At present, a plurality of control methods such as master control, master-slave control, coupling control and virtual main shaft control are provided, and although a certain control effect is obtained, the defects still exist: if the main command controls each motion shaft to work in parallel, the motion shafts are mutually incoherent, when one shaft is disturbed, the motion shafts can be adjusted only by the shaft, and other shafts can not respond, so that the motion shafts are only suitable for the places with less disturbance; the master-slave control is similar to master control, but the motion axis is divided into a main shaft and an auxiliary shaft, the reference signal of the auxiliary shaft is output from the main shaft, when the main shaft is disturbed, the auxiliary shaft can make corresponding adjustment, and when the auxiliary shaft is disturbed, the main shaft can not make corresponding response; the coupling control solves the problem that the motors are not coupled in the method, but the overall stability of the system is poor due to the introduction of parameter coupling between shafts; the virtual spindle mimics the physical properties of a mechanically stiff shaft connection and thus has similar inherent synchronization properties, but the virtual spindle may be out of synchronization during load upsets, start-up or shut-down, etc. Therefore, the method for researching more effective synchronous control of the rotating speed aiming at the multi-switch reluctance motor has important practical significance.
Disclosure of Invention
In order to solve the technical problems, the invention provides a synchronous control device for the rotating speed of a multi-switch reluctance motor, which has a simple structure and is accurate to control, and provides a synchronous control method for the rotating speed of the multi-switch reluctance motor.
The technical scheme for solving the problems is as follows: the multi-switch reluctance motor rotating speed synchronous control device comprises a rotating speed given module, a virtual main shaft controller, a magnetic linkage given module and a plurality of motor control modules, wherein the output end of the rotating speed given module is connected with the input end of the virtual main shaft controller; each motor control module comprises a motor rotating speed detection module, a motor rotating speed comparison module, a motor active disturbance rejection control module, a motor torque estimation module, a motor actual torque and reference torque comparison module, a motor torque compensation module, an adder, a motor flux linkage estimation module, a motor flux linkage comparison module and a motor controller; in each motor control module, the input end of the motor rotation speed detection module is connected with a corresponding motor, the output end of the motor rotation speed detection module is connected with the input end of the motor rotation speed comparison module, the motor torque compensation module and the virtual main shaft controller, the output end of the virtual main shaft controller is connected with the input end of the motor rotation speed comparison module and the input end of the motor torque compensation module, the output end of the motor active disturbance rejection control module is connected with the input end of the motor actual torque and the reference torque comparison module, the output end of the motor torque compensation module is connected with the input end of the adder, the input end of the motor torque estimation module is connected with the input end of the reference torque comparison module, the output end of the motor actual torque and the reference torque comparison module is connected with the input end of the adder, the output end of the adder is connected with the input end of the motor controller, the input end of the motor magnetic link estimation module is connected with the corresponding motor, and the output end of the motor magnetic link estimation module is connected with the motor magnetic link comparison module; the output end of the flux linkage given module is connected with the input end of the flux linkage comparison module in each motor control module.
A synchronous control method for the rotating speed of a multi-switch reluctance motor comprises the following steps:
step one: setting a given rotational speed ω of the system * Simultaneously detecting the actual rotational speed omega of each motor i (i=1, 2 … n), n represents the number of motors, and the virtual spindle rotational speed ω is obtained by processing the number of motors by the virtual spindle controller r And takes the same as a given rotating speed of each motor; the virtual main shaft controller is used for controlling the virtual main shaft according to the actual rotation speed omega of each motor i (i=1, 2 … n) and the corresponding moment of inertia j i (i=1, 2 … n), and obtaining the weight average rotation speed omega of each motor through calculation w And taking the torque compensation module as a reference rotating speed of each motor torque compensation module;
step two: the actual rotation speed omega of each motor i And virtual spindle speed omega r Comparing, and obtaining the reference torque T of each motor after the deviation is processed by the operation of the active disturbance rejection controller i * (i=1, 2 … n) while the actual torque T of each motor is obtained by the torque estimation module i (i=1, 2 … n), and the actual torque T is calculated i With the above mentioned ginsengTest torque T i * Comparing to obtain corresponding torque deviation delta T i (i=1,2…n);
Step three: the torque compensation modules of the motors are based on the actual rotational speed omega i And weight average rotation speed omega w The corresponding compensation torque delta T is obtained through the operation processing i ' (i=1, 2 … n) and compensating the torque Δt i ' deviation from the torque DeltaT i Summing to obtain the torque adjustment quantity delta T of the motor i ″(i=1,2…n);
Step four: estimating the actual flux linkage ψ of each motor i (i=1, 2 … n) and is given a flux linkage ψ with the system * Comparing to obtain the flux linkage deviation delta phi of each motor i (i=1,2…n);
Step five: adjusting the torque by an amount DeltaT i "sum flux linkage bias Δψ i After the operation processing of the motor controller, the output control signal controls the switching state of the power switch in the power converter corresponding to the motor, thereby realizing the actual rotation speed of each motor to the given rotation speed omega of the system * The aim of synchronous control of the rotating speed of the multi-switch reluctance motor is achieved.
In the above method for synchronously controlling the rotation speeds of the multiple switched reluctance motors, the virtual spindle controller obtains the virtual spindle rotation speed omega in the first step r The method of (1) comprises:
1-1-1) setting the system given rotational speed ω * Simultaneously detecting the initial rotation speed omega output by the virtual main shaft controller r ' and correlating it with a system given rotational speed ω * Comparing, the deviation of the torque T is obtained by Hooke's law r As shown in formula (1):
T r =b r*r ′)+k r ∫(ω *r ′)dt (1)
wherein: b r Attenuation coefficient, k, of virtual principal axis r Is the virtual main shaft elastic coefficient;
1-1-2) the actual rotational speed ω of each motor i With a given rotational speed omega of the system * Comparing, the deviation of the motor is obtained by Hooke's law to obtain the drive of the virtual transmission shaft of each motorDynamic torque, as shown in formula (2):
T i ′=b i*i )+k i ∫(ω *i )dt (2)
wherein: b i Attenuation coefficient k of motor i transmission shaft i The elastic coefficient of the transmission shaft of the motor i;
1-1-3) comparing the driving torque of the virtual main shaft with the sum of the driving torques of the virtual transmission shafts of the motors, and obtaining the adjusted virtual main shaft rotating speed omega by the rigid fixed shaft rotating law r The method comprises the following steps:
Figure BDA0001309219380000041
wherein: j (j) r The moment of inertia of the virtual main shaft;
1-1-4) adjusting the virtual main shaft rotating speed omega r As a given rotational speed of each motor.
According to the synchronous control method for the rotating speeds of the multi-switch reluctance motor, the drosophila optimization algorithm is adopted to optimize the virtual main shaft attenuation coefficient b r And virtual principal axis elastic coefficient k r The method is characterized by comprising the following steps:
with virtual principal axis attenuation coefficient b r And virtual principal axis elastic coefficient k r To optimize the object, the rotational speed omega is set by the system * And virtual spindle speed omega r Deviation delta omega of (2) r System given rotational speed omega * And weight average rotation speed omega w Deviation delta omega of (2) w For optimization purposes, the optimization steps are as follows:
(1) the method comprises the following steps Setting optimization objects b respectively r And k is equal to r The initial individual position of the fruit fly is X c (x c0 ,y c0 ) And X c ′(x′ c0 ,y′ c0 ) (c is the number of drosophila individuals, c=1, 2 … n'), while setting the maximum number of iterations to Maxgen;
(2) the method comprises the following steps The searching direction and distance of randomly generated drosophila are respectively shown as the formula (4) and the formula (5):
Figure BDA0001309219380000051
Figure BDA0001309219380000052
wherein: d (D) xc 、D yc 、D′ xc And D' yc Is a random value;
(3) the method comprises the following steps Taking the reciprocal of the distance from the original point of the drosophila individual as the judgment value S of the taste concentration c And S is c ' its expression is:
Figure BDA0001309219380000053
Figure BDA0001309219380000054
(4) the method comprises the following steps Determining a set of taste concentration values
Figure BDA0001309219380000055
Substituting the synchronous control model of the virtual main shaft of the multi-switch reluctance motor to obtain corresponding delta omega rc And Δω wc
(5) The method comprises the following steps Constructing an optimization objective function H (c) as shown in a formula (8):
Figure BDA0001309219380000056
wherein: d, d 1 And d 2 As the weight coefficient, d 1 >0,d 2 >0 and d 1 +d 2 =1, take d 1 =d 2 =0.5;
(6) The method comprises the following steps Maintaining the maximum value of formula (8), i.e. the current highest taste concentration, and preserving the position of the drosophila population at that time;
(7) the method comprises the following steps Performing iterative optimization, repeating the steps (2) - (6), judging whether the taste concentration of a new drosophila individual is better than the current highest taste concentration value, and if so, updating the current highest taste concentration and the initial position of the drosophila population; otherwise, returning to the step (2), and executing the step (8) until the current iteration times are equal to the maximum iteration times Maxgen or the target precision requirement is reached;
(8) the method comprises the following steps After the iterative optimization is finished, the optimal taste concentration value and the position of the drosophila at the moment are reserved, and the optimal virtual main axis attenuation coefficient b is obtained r And virtual principal axis elastic coefficient k r
In the above method for synchronously controlling the rotation speeds of the multiple switched reluctance motors, the virtual spindle controller in the step (1) calculates the weight average rotation speed ω w The method of (1) is as follows:
1-2-1) determining the weight coefficient g of each motor i
Based on moment of inertia j of each motor i (i=1, 2 … n), and a weight coefficient g of each motor is determined i (i=1, 2 … n), i.e.:
Figure BDA0001309219380000061
1-2-2) calculating the weight average rotation speed omega w
According to the actual rotation speed omega of each motor i And its corresponding weight coefficient g i Obtaining the weight average rotating speed omega of the system w As shown in formula (10):
Figure BDA0001309219380000062
in the above method for synchronously controlling the rotation speeds of the multiple switch reluctance motors, the reference torque T is in the second step i * The acquisition method of (1) is as follows:
2-1) at the actual rotational speed ω of the motor i And virtual spindle speed omega r Deviation delta omega of (2) i As input of the motor i active disturbance rejection control module, the total disturbance y of the motor i i Performing real-time estimation as shown in formula (11):
y i =-β i1 fal(Δω i1 ,η) (11)
wherein: y is i Representing the total disturbance suffered by the motor i in the running process; beta i1 (i=1, 2 … n) is a gain coefficient; function fal (Deltaomega) i1 η) is represented by the expression (12):
Figure BDA0001309219380000063
wherein: parameter alpha 1 Is a constant of 0 to 1, and is generally alpha 1 =0.25; the parameter eta is a constant affecting the filtering effect, and eta=0.5 is taken;
2-2) from Δω i The nonlinear error feedback control law is obtained as follows:
u i0 (t)=β i2 fal(Δω i2 ,η) (13)
wherein: beta i2 (i=1, 2 … n) is a gain coefficient; function fal (Deltaomega) i2 η) is represented by the expression (14):
Figure BDA0001309219380000071
wherein: parameter alpha 2 Constant of 0-1, alpha is taken 2 =0.75;
2-3) obtaining a reference torque T of the motor i according to the formulas (11) and (13) i * The method comprises the following steps:
T i * =u i0 (t)-y i /e i (15)
wherein: e, e i (i=1, 2 … n) is a disturbance compensation coefficient.
According to the synchronous control method for the rotating speeds of the multi-switch reluctance motor, the gain coefficient beta is optimized by adopting a Drosophila optimization algorithm i1 、β i2 And disturbance compensation coefficient e i The method is characterized by comprising the following steps:
gain coefficient beta of ith motor active disturbance rejection controller i1 、β i2 And disturbance compensation coefficient e i For optimizing the object, the rotation speed deviation delta omega is used i For optimization purposes, youThe method comprises the following steps:
2-1-1): setting optimization parameters beta respectively i1 、β i2 And e i The initial individual position of the fruit fly is X e (x e0 ,y e0 ),X e ′(x′ e0 ,y′ e0 ) And X' e (x″ e0 ,y″ e0 ) (e is the number of drosophila individuals, e=1, 2 … n "), while setting the maximum number of iterations to Maxgen';
2-1-2): randomly generating search direction and distance of drosophila:
Figure BDA0001309219380000072
Figure BDA0001309219380000073
Figure BDA0001309219380000074
wherein: d (D) xe 、D ye 、D′ xe 、D′ ye 、D″ xe And D' ye Is a random value
2-1-3): taking the reciprocal of the distance from the original point of the drosophila individual as the judgment value S of the taste concentration e 、S e ' and S e The expressions are as follows:
Figure BDA0001309219380000081
Figure BDA0001309219380000082
Figure BDA0001309219380000083
2-1-4): determining the taste concentration
Figure BDA0001309219380000084
Substituting the taste concentration determination values into the active disturbance rejection controller, simulating the active disturbance rejection controller, and determining the taste concentration W corresponding to each taste concentration determination value e ,W e The expression of (2) is shown in the formula (22):
Figure BDA0001309219380000085
wherein:
Figure BDA0001309219380000086
representing taking a set of taste concentration determination values +.>
Figure BDA0001309219380000087
When the noise is generated, the error between the input and the output of the active disturbance rejection controller is generated;
2-1-5): maintaining the maximum value of formula (22), i.e., the current highest taste concentration, and maintaining the position of the drosophila at that time;
2-1-6): performing iterative optimization, repeatedly executing the steps 2-1-2) -2-1-5, judging whether the taste concentration of the drosophila individuals is better than the current highest taste concentration value, and if so, updating the current highest taste concentration of the drosophila and the initial position of the drosophila population; otherwise, returning to the step 2-1-2), and executing the step 2-1-7 until the current iteration times are equal to the maximum iteration times Maxgen' or the target required precision is reached;
step 2-1-7): after the iterative optimization is finished, the optimal taste concentration value and the position of the drosophila individual at the moment are reserved, and the optimal active disturbance rejection control parameter beta of the motor i is obtained i1 、β i2 、e i
Adopting a drosophila optimization algorithm to calculate a control parameter beta of an ith motor i1 、β i2 E i And then, obtaining control parameters of other motors by using a time scale method, wherein the method comprises the following steps of:
i: obtaining a time scale p of the ith motor according to a state equation of phase current and rotating speed of the ith motor i As shown in formula (23):
Figure BDA0001309219380000091
wherein: i.e in Rated current of the ith motor; n is n i0 The rated rotation speed of the ith motor; j (j) i The rotational inertia of the ith motor; b (B) i The friction coefficient of the ith motor; l (L) imin The minimum value of the inductance, namely the phase inductance, of the ith motor stator salient pole at the position where the stator salient pole coincides with the center of the rotor groove;
Figure BDA0001309219380000092
the change rate of the phase inductance of the ith motor along with the position angle is set;
II: then, the time scale p of the kth motor (k=1 to n, and k+.i) is calculated according to the formula (23) k By time scale p of motor i and motor k i And p k Control parameter beta of motor i i1 、β i2 、e i Obtaining the corresponding control parameter beta of the motor k k1 、β k2 、e k As shown in formulas (24) to (26):
Figure BDA0001309219380000093
Figure BDA0001309219380000094
e k =e i (26)。
in the above method for synchronously controlling the rotation speeds of the multiple switched reluctance motors, in the third step, the ith motor compensates the torque delta T i The' acquisition mode is as follows:
1) At the actual rotational speed omega of the motor i i And weight average rotation speed omega w As state variables, the deviations and the integral thereof are represented by the following formulas (27) and (28):
x i1 =ω wi (i=1,2…n) (27)
Figure BDA0001309219380000095
2) Setting an integral sliding mode surface function as shown in formula (29):
s i =x i1 +Cx i2 (i=1,2…n) (29)
wherein: c is a positive constant;
3) The torque compensator is designed according to the established sliding mode surface function by adopting an exponential approach law, and the adopted exponential approach law expression is as follows:
Figure BDA0001309219380000101
wherein: epsilon and K are normal numbers, sgn(s) i ) Is a sign function;
4) Determining the compensation torque delta T according to the motion equation of the switched reluctance motor and (29) and (30) i ' is:
Figure BDA0001309219380000102
the invention has the beneficial effects that: the invention aims at a rotating speed synchronous control system formed by multiple switch reluctance motors, adopts a virtual main shaft control method and combines torque compensation to realize synchronous operation of the multiple motors. The basic principle is as follows: setting a given rotating speed of the system, detecting the actual rotating speed of each motor at the same time, obtaining a virtual spindle rotating speed through operation processing of a virtual spindle controller, and taking the virtual spindle rotating speed as the given rotating speed of each motor; meanwhile, the virtual main shaft controller calculates the weight average rotating speed of the system according to the actual rotating speed and the rotating inertia of each motor, and takes the weight average rotating speed as the reference rotating speed of each motor torque compensation module; then, the deviation of the actual rotating speed of each motor and the given rotating speed is processed by an active disturbance rejection control module to obtain reference torque of each motor, and the reference torque is compared with the actual torque of the motor to obtain corresponding torque deviation; meanwhile, each motor torque compensation module obtains corresponding compensation torque according to the actual rotating speed of the motor and the average rotating speed of the system weight, and sums the compensation torque and the torque deviation to obtain the torque adjustment quantity of the motor; in addition, the actual flux linkage of each motor is compared with the given flux linkage of the system to obtain corresponding flux linkage deviation, and finally the motor controller controls the motors according to the torque adjustment quantity and the flux linkage deviation, so that the accurate tracking of the actual rotating speed of each motor to the given rotating speed of the system can be realized, and the aim of synchronous operation of the rotating speeds of multiple motors is fulfilled.
Drawings
Fig. 1 is a block diagram of a control device according to the present invention.
FIG. 2 is a flow chart of the virtual spindle controller control according to the present invention.
FIG. 3 is a flow chart of a control method of the present invention.
FIG. 4 is a flow chart of the control parameter optimization of the active disturbance rejection controller according to the present invention.
Detailed Description
The invention is further described below with reference to the drawings and examples.
The multi-switch reluctance motor rotating speed synchronous control device comprises a rotating speed given module, a virtual main shaft controller, a magnetic linkage given module and a plurality of motor control modules, wherein the output end of the rotating speed given module is connected with the input end of the virtual main shaft controller; each motor control module comprises a motor rotating speed detection module, a motor rotating speed comparison module, a motor active disturbance rejection control module, a motor torque estimation module, a motor actual torque and reference torque comparison module, a motor torque compensation module, an adder, a motor flux linkage estimation module, a motor flux linkage comparison module and a motor controller; in each motor control module, the input end of the motor rotation speed detection module is connected with a corresponding motor, the output end of the motor rotation speed detection module is connected with the input end of the motor rotation speed comparison module, the motor torque compensation module and the virtual main shaft controller, the output end of the virtual main shaft controller is connected with the input end of the motor rotation speed comparison module and the input end of the motor torque compensation module, the output end of the motor active disturbance rejection control module is connected with the input end of the motor actual torque and the reference torque comparison module, the output end of the motor torque compensation module is connected with the input end of the adder, the input end of the motor torque estimation module is connected with the input end of the reference torque comparison module, the output end of the motor actual torque and the reference torque comparison module is connected with the input end of the adder, the output end of the adder is connected with the input end of the motor controller, the input end of the motor magnetic link estimation module is connected with the corresponding motor, and the output end of the motor magnetic link estimation module is connected with the motor magnetic link comparison module; the output end of the flux linkage given module is connected with the input end of the flux linkage comparison module in each motor control module.
As shown in fig. 2-4, a method for synchronously controlling the rotation speed of a multi-switch reluctance motor comprises the following steps:
step one: setting a given rotational speed ω of the system * Simultaneously detecting the actual rotational speed omega of each motor i (i=1, 2 … n), n represents the number of motors, and the virtual spindle rotational speed ω is obtained by processing the number of motors by the virtual spindle controller r And takes the same as a given rotating speed of each motor; the virtual main shaft controller is used for controlling the virtual main shaft according to the actual rotation speed omega of each motor i (i=1, 2 … n) and the corresponding moment of inertia j i (i=1, 2 … n), and obtaining the weight average rotation speed omega of each motor through calculation w And takes the reference rotation speed as the reference rotation speed of each motor torque compensation module.
The virtual spindle controller obtains the virtual spindle rotation speed omega r The method of (1) comprises:
1-1-1) setting the system given rotational speed ω * Simultaneously detecting the initial rotation speed omega output by the virtual main shaft controller r ' and correlating it with a system given rotational speed ω * Comparing, the deviation of the torque T is obtained by Hooke's law r As shown in formula (1):
T r =b r*r ′)+k r ∫(ω *r ′)dt (1)
wherein the method comprises the steps of:b r Attenuation coefficient, k, of virtual principal axis r Is the virtual main shaft elastic coefficient;
1-1-2) the actual rotational speed ω of each motor i With a given rotational speed omega of the system * Comparing, the deviation of the torque is obtained by Hooke's law to obtain the driving torque of the virtual transmission shaft of each motor, as shown in the formula (2):
T i ′=b i*i )+k i ∫(ω *i )dt (2)
wherein: b i Attenuation coefficient k of motor i transmission shaft i The elastic coefficient of the transmission shaft of the motor i;
1-1-3) comparing the driving torque of the virtual main shaft with the sum of the driving torques of the virtual transmission shafts of the motors, and obtaining the adjusted virtual main shaft rotating speed omega by the rigid fixed shaft rotating law r The method comprises the following steps:
Figure BDA0001309219380000121
wherein: j (j) r The moment of inertia of the virtual main shaft;
1-1-4) adjusting the virtual main shaft rotating speed omega r As a given rotational speed of each motor.
Optimizing the virtual principal axis attenuation coefficient b by adopting a Drosophila optimization algorithm r And virtual principal axis elastic coefficient k r The method is characterized by comprising the following steps:
with virtual principal axis attenuation coefficient b r And virtual principal axis elastic coefficient k r To optimize the object, the rotational speed omega is set by the system * And virtual spindle speed omega r Deviation delta omega of (2) r System given rotational speed omega * And weight average rotation speed omega w Deviation delta omega of (2) w For optimization purposes, the optimization steps are as follows:
(1) the method comprises the following steps Setting optimization objects b respectively r And k is equal to r The initial individual position of the fruit fly is X c (x c0 ,y c0 ) And X c ′(x′ c0 ,y′ c0 ) (c is the number of drosophila individuals, c=1, 2 … n') while setting the maximumThe iteration number is Maxgen;
(2) the method comprises the following steps The searching direction and distance of randomly generated drosophila are respectively shown as the formula (4) and the formula (5):
Figure BDA0001309219380000131
Figure BDA0001309219380000132
wherein: d (D) xc 、D yc 、D′ xc And D' yc Is a random value;
(3) the method comprises the following steps Taking the reciprocal of the distance from the original point of the drosophila individual as the judgment value S of the taste concentration c And S is c ' its expression is:
Figure BDA0001309219380000133
Figure BDA0001309219380000134
(4) the method comprises the following steps Determining a set of taste concentration values
Figure BDA0001309219380000135
Substituting the synchronous control model of the virtual main shaft of the multi-switch reluctance motor to obtain corresponding delta omega rc And Δω wc ;Δω rc And Δω wc Respectively represent S c ,S c After' substitution into the model, the system gives the rotation speed omega * And virtual spindle speed omega r Deviation of (2) and system given rotational speed omega * And weight average rotation speed omega w Deviation of (2);
(5) the method comprises the following steps Constructing an optimization objective function H (c) representing the optimization objective of the optimization step, i.e. to achieve Δω rc And Δω wc The minimum purpose is as shown in formula (8):
Figure BDA0001309219380000141
wherein: d, d 1 And d 2 As the weight coefficient, d 1 >0,d 2 >0 and d 1 +d 2 =1, take d 1 =d 2 =0.5;
(6) The method comprises the following steps Maintaining the maximum value of formula (8), i.e. the current highest taste concentration, and preserving the position of the drosophila population at that time;
(7) the method comprises the following steps Performing iterative optimization, repeating the steps (2) - (6), judging whether the taste concentration of a new drosophila individual is better than the current highest taste concentration value, and if so, updating the current highest taste concentration and the initial position of the drosophila population; otherwise, returning to the step (2), and executing the step (8) until the current iteration times are equal to the maximum iteration times Maxgen or the target precision requirement is reached;
(8) the method comprises the following steps After the iterative optimization is finished, the optimal taste concentration value and the position of the drosophila at the moment are reserved, and the optimal virtual main axis attenuation coefficient b is obtained r And virtual principal axis elastic coefficient k r
The virtual main shaft controller calculates the weight average rotating speed omega w The method of (1) is as follows:
1-2-1) determining the weight coefficient g of each motor i
Based on moment of inertia j of each motor i (i=1, 2 … n), and a weight coefficient g of each motor is determined i (i=1, 2 … n), i.e.:
Figure BDA0001309219380000142
1-2-2) calculating the weight average rotation speed omega w
According to the actual rotation speed omega of each motor i And its corresponding weight coefficient g i Obtaining the weight average rotating speed omega of the system w As shown in formula (10):
Figure BDA0001309219380000151
step two: the actual rotation speed omega of each motor i And virtual spindle speed omega r Comparing, and obtaining the reference torque T of each motor after the deviation is processed by the operation of the active disturbance rejection controller i * (i=1, 2 … n) while the actual torque T of each motor is obtained by the torque estimation module i (i=1, 2 … n), and the actual torque T is calculated i With the above reference torque T i * Comparing to obtain corresponding torque deviation delta T i (i=1,2…n)。
Reference torque T i * The acquisition method of (1) is as follows:
2-1) at the actual rotational speed ω of the motor i And virtual spindle speed omega r Deviation delta omega of (2) i As input of the motor i active disturbance rejection control module, the total disturbance y of the motor i i Performing real-time estimation as shown in formula (11):
y i =-β i1 fal(Δω i1 ,η) (11)
wherein: y is i Representing the total disturbance suffered by the motor i in the running process; beta i1 (i=1, 2 … n) is a gain coefficient; function fal (Deltaomega) i1 η) is represented by the expression (12):
Figure BDA0001309219380000152
wherein: parameter alpha 1 Is a constant of 0 to 1, and is generally alpha 1 =0.25; the parameter eta is a constant affecting the filtering effect, and eta=0.5 is taken;
2-2) from Δω i The nonlinear error feedback control law is obtained as follows:
u i0 (t)=β i2 fal(Δω i2 ,η) (13)
wherein: beta i2 (i=1, 2 … n) is a gain coefficient; function fal (Deltaomega) i2 η) is represented by the expression (14):
Figure BDA0001309219380000153
wherein: parameter alpha 2 Is a constant of 0 to 1, and is generally alpha 2 =0.75;
2-3) obtaining a reference torque T of the motor i according to the formulas (11) and (13) i * The method comprises the following steps:
T i * =u i0 (t)-y i /e i (15)
wherein: e, e i (i=1, 2 … n) is a disturbance compensation coefficient.
As shown in fig. 4, a drosophila optimization algorithm is used to optimize the gain factor beta i1 、β i2 And disturbance compensation coefficient e i The method is characterized by comprising the following steps:
gain coefficient beta of ith motor active disturbance rejection controller i1 、β i2 And disturbance compensation coefficient e i For optimizing the object, the rotation speed deviation delta omega is used i For optimization purposes, the optimization steps are as follows:
2-1-1): setting optimization parameters beta respectively i1 、β i2 And e i The initial individual position of the fruit fly is X e (x e0 ,y e0 ),X e ′(x′ e0 ,y′ e0 ) And X' e (x″ e0 ,y″ e0 ) (e is the number of drosophila individuals, e=1, 2 … n "), while setting the maximum number of iterations to Maxgen';
2-1-2): randomly generating search direction and distance of drosophila:
Figure BDA0001309219380000161
Figure BDA0001309219380000162
Figure BDA0001309219380000163
wherein: d (D) xe 、D ye 、D′ xe 、D′ ye 、D″ xe And D' ye Is a random value
2-1-3): taking the reciprocal of the distance from the original point of the drosophila individual as the judgment value S of the taste concentration e 、S e ' and S e The expressions are as follows:
Figure BDA0001309219380000164
Figure BDA0001309219380000165
Figure BDA0001309219380000166
2-1-4): determining a set of taste concentration values
Figure BDA0001309219380000167
Substituting the taste concentration W into the active disturbance rejection controller, simulating the active disturbance rejection controller, and determining the taste concentration W corresponding to each taste concentration determination value according to the simulation result e ,W e The expression of (2) is shown in the formula (22):
Figure BDA0001309219380000171
wherein:
Figure BDA0001309219380000172
representing taking a set of taste concentration determination values +.>
Figure BDA0001309219380000173
When the noise is generated, the error between the input and the output of the active disturbance rejection controller is generated;
2-1-5): maintaining the maximum value of formula (22), i.e., the current highest taste concentration, and maintaining the position of the drosophila at that time;
2-1-6): performing iterative optimization, repeatedly executing the steps 2-1-2) -2-1-5, judging whether the taste concentration of the drosophila individuals is better than the current highest taste concentration value, and if so, updating the current highest taste concentration of the drosophila and the initial position of the drosophila population; otherwise, returning to the step 2-1-2), and executing the step 2-1-7 until the current iteration times are equal to the maximum iteration times Maxgen' or the target precision requirement is reached;
step 2-1-7): after the iterative optimization is finished, the optimal taste concentration value and the position of the drosophila individual at the moment are reserved, and the optimal active disturbance rejection control parameter beta of the motor i is obtained i1 、β i2 、e i
Adopting a drosophila optimization algorithm to calculate a control parameter beta of an ith motor i1 、β i2 E i And then, obtaining control parameters of other motors by using a time scale method, wherein the method comprises the following steps of:
i: obtaining a time scale p of the ith motor according to a state equation of phase current and rotating speed of the ith motor i As shown in formula (23):
Figure BDA0001309219380000174
wherein: i.e in Rated current of the ith motor; n is n i0 The rated rotation speed of the ith motor; j (j) i The rotational inertia of the ith motor; b (B) i The friction coefficient of the ith motor; l (L) imin The minimum value of the inductance, namely the phase inductance, of the ith motor stator salient pole at the position where the stator salient pole coincides with the center of the rotor groove;
Figure BDA0001309219380000181
the change rate of the phase inductance of the ith motor along with the position angle is set;
II: then, the time scale p of the kth motor (k=1 to n, and k+.i) is calculated according to the formula (23) k By time scale p of motor i and motor k i And p k Control parameter beta of motor i i1 、β i2 、e i Obtaining the corresponding control parameter beta of the motor k k1 、β k2 、e k As shown in formulas (24) to (26):
Figure BDA0001309219380000182
Figure BDA0001309219380000183
e k =e i (26)
step three: the torque compensation modules of the motors are based on the actual rotational speed omega i And weight average rotation speed omega w The corresponding compensation torque delta T is obtained through the operation processing i ' (i=1, 2 … n) and compensating the torque Δt i ' deviation from the torque DeltaT i Summing to obtain the torque adjustment quantity delta T of the motor i ″(i=1,2…n)。
Compensation torque delta T of ith motor i The' acquisition mode is as follows:
1) At the actual rotational speed omega of the motor i i And weight average rotation speed omega w As state variables, the deviations and the integral thereof are represented by the following formulas (27) and (28):
x i1 =ω wi (i=1,2…n) (27)
Figure BDA0001309219380000184
2) Setting an integral sliding mode surface function as shown in formula (29):
s i =x i1 +Cx i2 (i=1,2…n) (29)
wherein: c is a positive constant;
3) Selecting an exponential approach law according to the established sliding mode surface function to design a torque compensator, wherein the selected exponential approach law expression is as follows:
Figure BDA0001309219380000191
wherein: epsilon and K are normal numbers, sgn(s) i ) Is a sign function;
4) The state variables set for equations (27) and (28) are derived:
Figure BDA0001309219380000192
Figure BDA0001309219380000193
wherein: the equation of motion of the combined switched reluctance motor from equation (29) can be given by equation (34):
Figure BDA0001309219380000194
wherein: t (T) Li Representing the load torque.
The combination of formulas (29), (33), (34) can be obtained:
Figure BDA0001309219380000195
from the formulas (30) and (35), it is possible to obtain:
Figure BDA0001309219380000196
let the compensation torque delta T i ' is:
Figure BDA0001309219380000197
owing to the>
Figure BDA0001309219380000198
Part is far smaller than->
Figure BDA0001309219380000199
Thus can +.>
Figure BDA00013092193800001910
Omitted, finally outputQuantity DeltaT of i ' is:
Figure BDA00013092193800001911
step four: estimating the actual flux linkage ψ of each motor i (i=1, 2 … n) and is given a flux linkage ψ with the system * Comparing to obtain the flux linkage deviation delta phi of each motor i (i=1,2…n)。
Step five: the motor controller adjusts the quantity delta T according to the torque i "sum flux linkage bias Δψ i And the interval of the magnetic linkage queries the switching table to obtain a corresponding voltage vector, and then the switching state of a corresponding power switch in the power converter is determined according to the voltage vector, so that the actual rotating speed of each motor can realize the given rotating speed omega of the system * The aim of synchronous control of the rotating speed of the multi-switch reluctance motor is achieved.

Claims (6)

1. A synchronous control method for the rotating speed of a multi-switch reluctance motor is characterized by comprising the following steps: the multi-switch reluctance motor rotating speed synchronous control device comprises a rotating speed given module, a virtual main shaft controller, a flux linkage given module and a plurality of motor control modules, wherein the output end of the rotating speed given module is connected with the input end of the virtual main shaft controller; each motor control module comprises a motor rotating speed detection module, a motor rotating speed comparison module, a motor active disturbance rejection control module, a motor torque estimation module, a motor actual torque and reference torque comparison module, a motor torque compensation module, an adder, a motor flux linkage estimation module, a motor flux linkage comparison module and a motor controller; in each motor control module, the input end of the motor rotation speed detection module is connected with a corresponding motor, the output end of the motor rotation speed detection module is connected with the input end of the motor rotation speed comparison module, the motor torque compensation module and the virtual main shaft controller, the output end of the virtual main shaft controller is connected with the input end of the motor rotation speed comparison module and the input end of the motor torque compensation module, the output end of the motor active disturbance rejection control module is connected with the input end of the motor actual torque and the reference torque comparison module, the output end of the motor torque compensation module is connected with the input end of the adder, the input end of the motor torque estimation module is connected with the input end of the reference torque comparison module, the output end of the motor actual torque and the reference torque comparison module is connected with the input end of the adder, the output end of the adder is connected with the input end of the motor controller, the input end of the motor magnetic link estimation module is connected with the corresponding motor, and the output end of the motor magnetic link estimation module is connected with the motor magnetic link comparison module; the output end of the flux linkage given module is connected with the input end of the flux linkage comparison module in each motor control module;
the synchronous control method of the rotating speed of the multi-switch reluctance motor comprises the following steps:
step one: setting a given rotational speed ω of the system * Simultaneously detecting the actual rotational speed omega of each motor i I=1, 2 … n, n represents the number of motors, and the virtual spindle rotational speed ω is obtained after the virtual spindle controller operation processing r And takes the same as a given rotating speed of each motor; the virtual main shaft controller is used for controlling the virtual main shaft according to the actual rotation speed omega of each motor i I=1, 2 … n and their corresponding moment of inertia j i I=1, 2 … n, and obtaining the weight average rotation speed omega of each motor through calculation w And taking the torque compensation module as a reference rotating speed of each motor torque compensation module;
the virtual spindle controller obtains the virtual spindle rotation speed omega r The method of (1) comprises:
1-1-1) setting the system given rotational speed ω * Simultaneously detecting the initial rotation speed omega 'output by the virtual main shaft controller' r And compares it with a system given rotational speed omega * Comparing, the deviation of the torque T is obtained by Hooke's law r As shown in formula (1):
Figure FDA0004174581510000021
wherein: b r Attenuation coefficient, k, of virtual principal axis r Is the virtual main shaft elastic coefficient;
1-1-2) the actual rotational speed ω of each motor i With a given rotational speed omega of the system * Comparing, the deviation of the torque is obtained by Hooke's law to obtain the driving torque of the virtual transmission shaft of each motor, as shown in the formula (2):
Figure FDA0004174581510000022
wherein: b i Attenuation coefficient k of motor i transmission shaft i The elastic coefficient of the transmission shaft of the motor i;
1-1-3) comparing the driving torque of the virtual main shaft with the sum of the driving torques of the virtual transmission shafts of the motors, and obtaining the adjusted virtual main shaft rotating speed omega by the rigid fixed shaft rotating law r The method comprises the following steps:
Figure FDA0004174581510000023
wherein: j (j) r The moment of inertia of the virtual main shaft;
1-1-4) adjusting the virtual main shaft rotating speed omega r As a given rotational speed of each motor;
step two: the actual rotation speed omega of each motor i And virtual spindle speed omega r Comparing, and obtaining the reference torque T of each motor after the deviation is processed by the operation of the active disturbance rejection controller i * I=1, 2 … n, and the actual torque T of each motor is obtained by the torque estimation module i I=1, 2 … n, and the actual torque T is calculated i With the above reference torque T i * Comparing to obtain corresponding torque deviation delta T i ,i=1,2…n;
Step three: the torque compensation modules of the motors are based on the actual rotational speed omega i And weight average rotation speed omega w The corresponding compensation torque delta T is obtained through the operation processing i ' i=1, 2 … n, and the compensation torque Δt is set i ' deviation from the torque DeltaT i Summing to obtain the torque adjustment quantity delta T of the motor i ″,i=1,2…n;
Step four: estimating the actual flux linkage ψ of each motor i I=1, 2 … n and is given a flux linkage ψ with the system * Comparing to obtain the flux linkage deviation delta phi of each motor i ,i=1,2…n;
Step five: adjusting the torque by an amount DeltaT i "sum flux linkage bias Δψ i After the operation processing of the motor controller, the output control signal controls the switching state of the power switch in the power converter corresponding to the motor, thereby realizing the actual rotation speed of each motor to the given rotation speed omega of the system * The aim of synchronous control of the rotating speed of the multi-switch reluctance motor is achieved.
2. The method for synchronously controlling the rotation speeds of the multiple-switch reluctance motor according to claim 1, wherein: optimizing the virtual principal axis attenuation coefficient b by adopting a Drosophila optimization algorithm r And virtual principal axis elastic coefficient k r The method is characterized by comprising the following steps:
with virtual principal axis attenuation coefficient b r And virtual principal axis elastic coefficient k r To optimize the object, the rotational speed omega is set by the system * And virtual spindle speed omega r Deviation delta omega of (2) r System given rotational speed omega * And weight average rotation speed omega w Deviation delta omega of (2) w For optimization purposes, the optimization steps are as follows:
(1) the method comprises the following steps Setting optimization objects b respectively r And k is equal to r The initial individual position of the fruit fly is X c (x c0 ,y c0 ) And X' c (x′ c0 ,y′ c0 ) C is the number of drosophila individuals, c=1, 2 … n', and the maximum iteration number is set as Maxgen;
(2) the method comprises the following steps The searching direction and distance of randomly generated drosophila are respectively shown as the formula (4) and the formula (5):
Figure FDA0004174581510000031
Figure FDA0004174581510000041
wherein: d (D) xc 、D yc 、D′ xc And D' yc Is a random value;
(3) the method comprises the following steps Taking the reciprocal of the distance from the original point of the drosophila individual as the judgment value S of the taste concentration c And S' c The expressions are respectively:
Figure FDA0004174581510000042
Figure FDA0004174581510000043
(4) the method comprises the following steps Determining a set of taste concentration values
Figure FDA0004174581510000045
Substituting the synchronous control model of the virtual main shaft of the multi-switch reluctance motor to obtain corresponding delta omega rc And Δω wc
(5) The method comprises the following steps Constructing an optimized objective function H (c) as shown in formula (8):
Figure FDA0004174581510000044
wherein: d, d 1 And d 2 As the weight coefficient, d 1 >0,d 2 > 0 and d 1 +d 2 =1, take d 1 =d 2 =0.5;
(6) The method comprises the following steps Maintaining the maximum value of formula (8), i.e. the current highest taste concentration, and preserving the position of the drosophila population at that time;
(7) the method comprises the following steps Performing iterative optimization, repeating the steps (2) - (6), judging whether the taste concentration of a new drosophila individual is better than the current highest taste concentration value, and if so, updating the current highest taste concentration and the initial position of the drosophila population; otherwise, returning to the step (2), and executing the step (8) until the current iteration times are equal to the maximum iteration times Maxgen or the target precision requirement is reached;
(8) the method comprises the following steps After the iterative optimization is finished, the optimal taste concentration value and the position of the drosophila at the moment are reserved, and the optimal virtual main axis attenuation coefficient b is obtained r And virtual principal axis elastic coefficient k r
3. The method for synchronously controlling the rotation speeds of the multiple-switch reluctance motor according to claim 1, wherein: in the first step, the virtual spindle controller calculates a weight average rotating speed omega w The method of (1) is as follows:
1-2-1) determining the weight coefficient g of each motor i
Based on moment of inertia j of each motor i I=1, 2 … n, and determining the weight coefficient g of each motor i I=1, 2 … n, i.e.:
Figure FDA0004174581510000051
1-2-2) calculating the weight average rotation speed omega w
According to the actual rotation speed omega of each motor i And its corresponding weight coefficient g i Obtaining the weight average rotating speed omega of the system w As shown in formula (10):
Figure FDA0004174581510000052
4. the method for synchronously controlling the rotating speeds of the multi-switch reluctance motors according to claim 3, wherein: the reference torque T in the second step i * The acquisition method of (1) is as follows:
2-1) At the actual rotation speed omega of the motor i And virtual spindle speed omega r Deviation delta omega of (2) i As input of the motor i active disturbance rejection control module, the total disturbance y of the motor i i Performing real-time estimation as shown in formula (11):
y i =-β i1 fal(Δω i1 ,η) (11)
wherein: y is i Representing the total disturbance suffered by the motor i in the running process; beta i1 I=1, 2 … n as gain coefficient; function fal (Deltaomega) i1 η) is represented by the expression (12):
Figure FDA0004174581510000053
wherein: parameter alpha 1 Constant of 0-1, alpha is taken 1 =0.25; the parameter eta is a constant affecting the filtering effect, and eta=0.5 is taken;
2-2) from Δω i The nonlinear error feedback control law is obtained as follows:
u i0 (t)=β i2 fal(Δω i2 ,η) (13)
wherein: beta i2 I=1, 2 … n as gain coefficient; function fal (Deltaomega) i2 η) is represented by the expression (14):
Figure FDA0004174581510000061
wherein: parameter alpha 2 Constant of 0-1, alpha is taken 2 =0.75;
2-3) obtaining a reference torque T of the motor i according to the formulas (11) and (13) i * The method comprises the following steps:
T i * =u i0 (t)-y i /e i (15)
wherein: e, e i For the disturbance compensation coefficient, i=1, 2 … n.
5. The method for synchronously controlling the rotating speeds of the multiple-switch reluctance motors according to claim 4, wherein: optimization of the gain factor beta using a Drosophila optimization algorithm i1 、β i2 And disturbance compensation coefficient e i The method is characterized by comprising the following steps:
gain coefficient beta of ith motor active disturbance rejection controller i1 、β i2 And disturbance compensation coefficient e i For optimizing the object, the rotation speed deviation delta omega is used i For optimization purposes, the optimization steps are as follows:
2-1-1): setting optimization parameters beta respectively i1 、β i2 And e i The initial individual position of the fruit fly is X e (x e0 ,y e0 ),X′ e (x′ e0 ,y′ e0 ) And X' e (x″ e0 ,y″ e0 ) E is the number of drosophila individuals, e=1, 2 … n ", while the maximum number of iterations is set to Maxgen';
2-1-2): randomly generating search direction and distance of drosophila:
Figure FDA0004174581510000062
Figure FDA0004174581510000063
Figure FDA0004174581510000064
wherein: d (D) xe 、D ye 、D′ xe 、D′ ye 、D″ xe And D' ye Is a random value;
2-1-3): taking the reciprocal of the distance from the original point of the drosophila individual as the judgment value S of the taste concentration e 、S e ' and S e The expressions are as follows:
Figure FDA0004174581510000071
Figure FDA0004174581510000072
Figure FDA0004174581510000073
2-1-4): determining a set of taste concentration values
Figure FDA0004174581510000077
Substituting the taste concentration determination values into the active disturbance rejection controller, simulating the active disturbance rejection controller, and determining the taste concentration W corresponding to each taste concentration determination value e ,W e The expression of (2) is shown in the formula (22):
Figure FDA0004174581510000074
wherein:
Figure FDA0004174581510000075
representing taking a set of taste concentration determination values +.>
Figure FDA0004174581510000076
When the noise is generated, the error between the input and the output of the active disturbance rejection controller is generated;
2-1-5): maintaining the maximum value of formula (22), i.e., the current highest taste concentration, and maintaining the position of the drosophila at that time;
2-1-6): performing iterative optimization, repeatedly executing the steps 2-1-2) -2-1-5, judging whether the taste concentration of the drosophila individuals is better than the current highest taste concentration value, and if so, updating the current highest taste concentration of the drosophila and the initial position of the drosophila population; otherwise, returning to the step 2-1-2), and executing the step 2-1-7 until the current iteration times are equal to the maximum iteration times Maxgen' or the target required precision is reached;
step 2-1-7): after the iterative optimization is finished, the optimal taste concentration value and the position of the drosophila individual at the moment are reserved, and the optimal active disturbance rejection control parameter beta of the motor i is obtained i1 、β i2 、e i
Adopting a drosophila optimization algorithm to calculate a control parameter beta of an ith motor i1 、β i2 E i And then, obtaining control parameters of other motors by using a time scale method, wherein the method comprises the following steps of:
i: obtaining a time scale p of the ith motor according to a state equation of phase current and rotating speed of the ith motor i As shown in formula (23):
Figure FDA0004174581510000081
wherein: i.e in Rated current of the ith motor; n is n i0 The rated rotation speed of the ith motor; j (j) i The rotational inertia of the ith motor; b (B) i The friction coefficient of the ith motor; l (L) imin The minimum value of the inductance, namely the phase inductance, of the ith motor stator salient pole at the position where the stator salient pole coincides with the center of the rotor groove;
Figure FDA0004174581510000082
the change rate of the phase inductance of the ith motor along with the position angle is set;
II: then, according to the formula (23), calculating the time scale p of the kth motor k K=1 to n, and k+.i, by the time scales p of motor i and motor k i And p k Control parameter beta of motor i i1 、β i2 、e i Obtaining the corresponding control parameter beta of the motor k k1 、β k2 、e k As shown in formulas (24) to (26):
Figure FDA0004174581510000083
Figure FDA0004174581510000084
e k =e i (26)。
6. the method for synchronously controlling the rotating speeds of the multiple-switch reluctance motors according to claim 5, wherein: in the third step, the ith motor compensates the torque delta T i The' acquisition mode is as follows:
1) At the actual rotational speed omega of the motor i i And weight average rotation speed omega w As state variables, the deviations and the integral thereof are represented by the following formulas (27) and (28):
x i1 =ω wi (i=1,2…n) (27)
Figure FDA0004174581510000085
2) Setting an integral sliding mode surface function as shown in formula (29):
s i =x i1 +Cx i2 (i=1,2…n) (29)
wherein: c is a positive constant;
3) The torque compensator is designed according to the established sliding mode surface function by adopting an exponential approach law, and the adopted exponential approach law expression is as follows:
Figure FDA0004174581510000091
wherein: epsilon and K are normal numbers, sgn(s) i ) Is a sign function;
4) Determining the compensation torque delta T according to the motion equation of the switched reluctance motor and (29) and (30) i ' is:
Figure FDA0004174581510000092
/>
CN201710398938.7A 2017-05-31 2017-05-31 Multi-switch reluctance motor rotating speed synchronous control device and control method Active CN107070315B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710398938.7A CN107070315B (en) 2017-05-31 2017-05-31 Multi-switch reluctance motor rotating speed synchronous control device and control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710398938.7A CN107070315B (en) 2017-05-31 2017-05-31 Multi-switch reluctance motor rotating speed synchronous control device and control method

Publications (2)

Publication Number Publication Date
CN107070315A CN107070315A (en) 2017-08-18
CN107070315B true CN107070315B (en) 2023-05-12

Family

ID=59616083

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710398938.7A Active CN107070315B (en) 2017-05-31 2017-05-31 Multi-switch reluctance motor rotating speed synchronous control device and control method

Country Status (1)

Country Link
CN (1) CN107070315B (en)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108733887B (en) * 2018-04-12 2020-01-10 湖南科技大学 Switched reluctance motor thermal analysis method based on variable density symmetrical grid division
CN109510521B (en) * 2018-11-06 2022-01-14 湖南工业大学 Multi-motor total quantity cooperative and consistent fault-tolerant control method based on virtual main shaft
CN109474207B (en) * 2018-11-22 2021-11-02 天津大学 Control method for improving synchronization performance of multi-axis motion system
CN109861618B (en) * 2019-01-11 2020-11-20 江苏大学 Construction method of anti-interference composite controller of BSG alternating current motor for hybrid electric vehicle
CN109687772A (en) * 2019-01-17 2019-04-26 武汉菲仕运动控制系统有限公司 A kind of multi-axis synchronized control method and system
CN109968051B (en) * 2019-03-28 2021-03-12 沈阳机床成套设备有限责任公司 Machine tool design method based on numerical control double-spindle synchronous machining
CN110456657B (en) * 2019-07-18 2022-10-04 中山瑞信智能控制系统有限公司 Simulation method of reluctance motor, embedded system and operation method
CN111555665B (en) * 2020-05-29 2021-11-23 中国铁建重工集团股份有限公司 Motor control method
CN111665824A (en) * 2020-06-23 2020-09-15 中国北方车辆研究所 Comparison test method for PI controller and linear active disturbance rejection controller
CN111884558B (en) * 2020-07-14 2021-09-14 湖南科技大学 Active heave compensation control method for marine winch driven by switched reluctance motor
CN113452287B (en) * 2021-06-07 2023-03-28 江苏科技大学 Control method and control system for multiple permanent magnet synchronous motors of underwater vehicle
CN113844281B (en) * 2021-10-09 2023-05-19 四川鼎鸿智电装备科技有限公司 Rotation speed synchronous control method and device, electronic equipment and storage medium
CN113965108B (en) * 2021-11-19 2023-07-25 江苏科技大学 Multi-motor collaborative propulsion system of underwater robot and control method
CN114679103B (en) * 2022-05-30 2022-12-06 广东工业大学 Sliding mode current compensation control system and method for switched reluctance motor

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6316893B1 (en) * 1998-05-20 2001-11-13 Abb Industry Oy Method and arrangement for adaptive load change compensation
CN1889358A (en) * 2006-07-11 2007-01-03 合肥工业大学 Permasyn electric machine control system based on adaptive sliding moding structure
CN102035456A (en) * 2010-12-14 2011-04-27 长春工业大学 Direct torque control system of permanent magnet synchronous motor based on terminal sliding mode
CN102611381A (en) * 2012-03-12 2012-07-25 浙江工业大学 Direct torque control system of permanent-magnet synchronous motor
CN105515455A (en) * 2016-01-28 2016-04-20 湖南科技大学 Synchronous control method and device for biswitch reluctance motors
CN106208865A (en) * 2016-08-10 2016-12-07 天津工业大学 Many permagnetic synchronous motors Virtual-shaft control method based on Load Torque Observer

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6316893B1 (en) * 1998-05-20 2001-11-13 Abb Industry Oy Method and arrangement for adaptive load change compensation
CN1889358A (en) * 2006-07-11 2007-01-03 合肥工业大学 Permasyn electric machine control system based on adaptive sliding moding structure
CN102035456A (en) * 2010-12-14 2011-04-27 长春工业大学 Direct torque control system of permanent magnet synchronous motor based on terminal sliding mode
CN102611381A (en) * 2012-03-12 2012-07-25 浙江工业大学 Direct torque control system of permanent-magnet synchronous motor
CN105515455A (en) * 2016-01-28 2016-04-20 湖南科技大学 Synchronous control method and device for biswitch reluctance motors
CN106208865A (en) * 2016-08-10 2016-12-07 天津工业大学 Many permagnetic synchronous motors Virtual-shaft control method based on Load Torque Observer

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
一种用于冶金轧钢传动系统的转矩补偿方法;甘韦韦等;《大功率变流技术》;20121231(第6期);第41页第1段至第44页最后1段,图3 *
永磁同步电机直接转矩控制系统的研究;张露等;《变频器世界》;20161031(第10期);第65页第1段至第68页最后1段,图1-3 *
齐宏等.果蝇.《辐射传输逆问题的智能优化理论与应用》.2016, *

Also Published As

Publication number Publication date
CN107070315A (en) 2017-08-18

Similar Documents

Publication Publication Date Title
CN107070315B (en) Multi-switch reluctance motor rotating speed synchronous control device and control method
CN107425769B (en) Active disturbance rejection control method and system for permanent magnet synchronous motor speed regulating system
CN108092567B (en) Permanent magnet synchronous motor rotating speed control system and method
CN105680746B (en) A method of designing permasyn morot electric current and speed control PI parameters using online particle swarm optimization algorithm
CN107070341A (en) Permagnetic synchronous motor method for suppressing torque ripple based on Robust Iterative Learning Control
CN104980069B (en) A kind of brshless DC motor double loop speed-regulating system Multipurpose Optimal Method
CN112054727B (en) Elevator permanent magnet synchronous motor drive control system
Shamseldin et al. Different techniques of self-tuning FOPID control for Brushless DC Motor
Zair et al. MRAS speed sensorless vector control of induction motor drives using predictive adaptation mechanism
CN112783099B (en) Fractional order composite control method and permanent magnet synchronous motor speed servo system
CN112054728B (en) Permanent magnet synchronous motor drive control system of numerical control machine tool
CN110266227A (en) A kind of control system for permanent-magnet synchronous motor based on fuzzy synovial membrane structure changes
CN108512476A (en) A kind of induction machine turn count method based on novel Luenberger observer
CN106685290B (en) Method for self-tuning of rotating speed controller parameters of permanent magnet synchronous motor based on chaotic molecular dynamics theory algorithm
Sayouti et al. Sensor less low speed control with ANN MRAS for direct torque controlled induction motor drive
CN113037161A (en) Model reference self-adaptive permanent magnet synchronous motor position-sensorless vector control method based on super-distortion sliding mode algorithm
CN115276491A (en) Speed optimization control method for permanent magnet synchronous motor
Diab et al. Particle swarm optimization based vector control of permanent magnet synchronous motor drive
CN114900086A (en) Sliding mode observation method for load torque of permanent magnet synchronous motor
Okoro et al. Performance assessment of a model-based DC motor scheme
Aydogdu et al. Kalman State Estimation and LQR Assisted Adaptive Control Of a Variable Loaded Servo System.
CN110048653B (en) Sliding mode variable structure MRAS rotating speed identification method of bearingless asynchronous motor
Ganjewar Novel approaches in sensorless induction motor drive for industrial applications
Vo et al. Pulse-width modulation direct torque control induction motor drive with Kalman filter
Zhang et al. Force Ripple Suppression of Permanent Magnet Linear Synchronous Motor Based on Fuzzy Adaptive Kalman Filter

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