CN111181468A - Method for determining control parameter stability domain of finite time control BBMC speed regulation system - Google Patents
Method for determining control parameter stability domain of finite time control BBMC speed regulation system Download PDFInfo
- Publication number
- CN111181468A CN111181468A CN202010064506.4A CN202010064506A CN111181468A CN 111181468 A CN111181468 A CN 111181468A CN 202010064506 A CN202010064506 A CN 202010064506A CN 111181468 A CN111181468 A CN 111181468A
- Authority
- CN
- China
- Prior art keywords
- bbmc
- control
- time
- phase
- speed regulating
- 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.)
- Granted
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P23/00—Arrangements or methods for the control of AC motors characterised by a control method other than vector control
- H02P23/0004—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
- H02P25/02—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P2207/00—Indexing scheme relating to controlling arrangements characterised by the type of motor
- H02P2207/01—Asynchronous machines
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
The invention discloses a method for determining a control parameter stability domain of a finite time control BBMC speed regulation system, which comprises the following steps: establishing a state differential equation of the BBMC speed regulation system by taking inductive current, capacitor voltage and output current in the BBMC as state variables; obtaining a dynamic equation of the system according to the state differential equation; designing a control function of the BBMC speed regulating system according to the dynamic equation of the system; obtaining the duty ratio of a power switch tube in the BBMC according to a control function and a finite time control principle; obtaining a discrete iteration mapping model of the BBMC inverter stage according to the duty ratio; establishing a discrete iteration mapping model of the three-phase asynchronous motor; obtaining a discrete iteration mapping model of the BBMC speed regulating system based on finite time control according to the discrete iteration mapping model of the BBMC inverter stage and the discrete iteration mapping model of the three-phase asynchronous motor; and obtaining the value range of the control parameter when the BBMC speed regulation system stably operates through numerical simulation according to the obtained discrete iterative mapping model.
Description
Technical Field
The invention relates to the field of asynchronous motor speed regulation, in particular to a method for determining a control parameter stability region of a finite time control BBMC speed regulation system.
Background
The Buck-Boost matrix converter (BBMC) is a novel power converter with high voltage transmission ratio and capable of directly outputting high-quality sine waves, and is suitable for being applied to an asynchronous motor speed regulating system as a power converter.
However, to realize high-performance speed regulation control of the BBMC-based asynchronous motor speed regulation system, an effective control strategy must be adopted for the research. The finite time control is a control method which can make a closed-loop control system converge in finite time, has extremely strong dynamic stability and disturbance resistance, and is very suitable for being used in the control of an asynchronous motor speed regulating system based on BBMC. When the speed regulating system is controlled by adopting limited time control, the value taking problem of a plurality of control parameters is involved, and if the control parameters are not selected properly, the expected control effect is difficult to achieve.
Disclosure of Invention
In order to solve the technical problem, the invention provides a method for determining a control parameter stability region of a limited time control BBMC speed regulation system, and the value range of relevant control parameters can be determined when the system stably operates.
The technical scheme for solving the technical problems comprises the following steps:
step (1): establishing a state differential equation of the BBMC speed regulation system by taking inductive current, capacitor voltage and output current in the BBMC as state variables;
step (2): obtaining a dynamic equation of the BBMC speed regulating system according to the state differential equation obtained in the step (1);
and (3): designing a control function of the BBMC speed regulating system according to the dynamic equation of the BBMC speed regulating system obtained in the step (2);
and (4): obtaining the duty ratio of a power switching tube in the BBMC according to the control function of the BBMC speed regulation system obtained in the step (3) and a finite time control principle;
and (5): obtaining a discrete iteration mapping model of the BBMC inverter stage according to the duty ratio of the power switch tube in the BBMC obtained in the step (4);
and (6): establishing a discrete iteration mapping model of the three-phase asynchronous motor;
and (7): obtaining a discrete iteration mapping model of the BBMC speed regulating system based on finite time control according to the discrete iteration mapping model of the BBMC inverter stage obtained in the step (5) and the discrete iteration mapping model of the three-phase asynchronous motor obtained in the step (6);
and (8): and (4) obtaining the value range of the control parameter when the BBMC speed regulation system stably operates through numerical simulation according to the discrete iterative mapping model of the BBMC speed regulation system based on the finite time control obtained in the step (7).
Preferably, in the step (1), the inductive current, the capacitor voltage and the output current in the BBMC are used as state variables, and a state differential equation of the asynchronous motor speed regulating system based on the BBMC is established, specifically:
wherein: i.e. iL、uCAnd ioRespectively representing inductive current, capacitor voltage and output current in BBMC, E is input voltage of BBMC inverter stage, u is output voltage of BBMC inverter stageDFor the common terminal voltage of three-phase stator winding of asynchronous motor, L is BBMC main circuit inductance, C is BBMC main circuit capacitance, R is1And L1Respectively is the equivalent resistance and the equivalent inductance of the single-phase winding of the asynchronous motor, d is the duty ratio of a power switch in the BBMC, and d belongs to [0,1 ]]。
Preferably, step (2) is to obtain a dynamic equation of the system according to the state differential equation obtained in step (1), specifically:
let the reference output voltage of BBMC be urefThe error of the output voltage of the BBMC is obtained as follows:
x1=uC-uref(2)
according to the formula (1) and the formula (2), the dynamic equation of the system is obtained as follows:
wherein: x is the number of1、x2、x3、iL、uCAnd ioAre all variables of time t.
Preferably, the step (3) designs a control function of the BBMC speed regulation system according to the dynamic equation of the system obtained in the step (2), specifically:
step (3-1): defining a time scale coordinate transformation as:
wherein: and t is Ks, and K is a coordinate transformation coefficient.
Step (3-2): and (3) according to the time scale coordinate transformation formula shown in the formula (4), transforming the formula (3) to obtain:
step (3-3): the control function of the BBMC governor system obtained according to equation (5) is:
preferably, the step (4) obtains the duty ratio of the power switching tube in the BBMC according to the control function and the finite time control principle of the BBMC speed control system obtained in the step (3), and specifically comprises:
firstly, according to the finite time control principle, a system control function f(s) is determined as follows:
wherein: satαAs a function of saturation, k1、k2、k3、α1、α2and alpha3Are control parameters.
Then, it is obtained from equations (6) and (7):
finally, according to the formula (4) and the formula (8), the duty ratio of the power switch tube in the BBMC is obtained as follows:
preferably, the step (5) obtains a discrete iterative mapping model of the BBMC inverter stage according to the duty ratio of the power switching tube in the BBMC obtained in the step (4), specifically:
the method comprises the following steps of establishing a differential equation of state by taking inductive current and capacitor voltage in BBMC as state variables:
according to the formula (9), the on-time T of the power switch tube in the nth switching period T can be obtainedonAnd off time toffRespectively is as follows:
wherein: dnAnd represents the duty ratio of the power switch tube in the BBMC in the nth switching period T.
According to equations (10) and (11), the discrete iterative mapping model for the BBMC inversion stage is obtained as follows:
wherein: u. ofC(n+1)And iL(n+1)Respectively representing the capacitor voltage and the inductor current i in the (n +1) T moment BBMCo(n)Denotes the output current of BBMC at time nT, M, N and ω are intermediate variables, anduC(n)and iL(n)Respectively representing the capacitor voltage and the inductor current at time nT BBMC.
Preferably, the step (6) of establishing a discrete iterative mapping model of the three-phase asynchronous motor specifically comprises:
establishing a state differential equation of the three-phase asynchronous motor on a two-phase static coordinate system as follows:
wherein: u. ofsαAnd usβRepresenting the stator voltage, i, of the motor in two-phase stationary coordinates, respectivelysαAnd isβRepresenting the stator currents of the motor in two stationary phases, ΨrαAnd ΨrβRespectively representing the rotor flux, T, of the motor in two-phase stationary coordinatesLRepresenting motor load torque, ωrRepresenting the angular velocity, L, of the rotor of the machines、Lr、Lm、Rs、Rr、npAnd J respectively represent the self inductance of the stator, the self inductance of the rotor, the mutual inductance of the stator and the rotor, the resistance of the stator, the resistance of the rotor, the pole pair number and the moment of inertia of the motor,representing the leakage coefficient, T, of the motorr=Lr/RrRepresenting the rotor electromagnetic time constant.
Discretizing the state differential equation shown in the formula (13) by a Runge-Kutta method to obtain:
wherein:is the state vector at time (n +1) T,is the state vector at time nT, K1、K2、K3And K4Are all intermediate variables, and K1=f(xn,yn),K4=f(xn+TK0,yn+TK3),
usα(n)And usβ(n)Respectively representing the stator voltage i of the motor at the moment nT on the two-phase stationary coordinatesα(n)And isβ(n)Representing the stator currents of the machine at the moment nT in two-phase stationary coordinates, Ψrα(n)And Ψrβ(n)Respectively representing the rotor flux linkage, omega, at time nT on the two-phase stationary coordinater(n)Representing the angular velocity of the rotor of the machine at time nT.
Preferably, the step (7) obtains a discrete iteration mapping model of the BBMC speed regulation system based on finite time control according to the discrete iteration mapping model of the BBMC inverter stage obtained in the step (5) and the discrete iteration mapping model of the three-phase asynchronous motor obtained in the step (6), and specifically comprises:
wherein:andrespectively representing the three-phase capacitor voltages at time (n +1) T BBMC,andrespectively represents three-phase inductive current in the BBMC at the (n +1) T moment, andrespectively represents the conduction time of the three-phase power switch tube in the n-th switch period T,andrespectively showing the turn-off time of the three-phase power switch tube in the BBMC in the nth switching period T,andrespectively representing the duty ratio of three-phase power switching tubes in the BBMC in the nth switching period T, MA、MB、MC、NA、NBAnd NCAre all intermediate variables, and
preferably, the step (8) obtains the value range of the control parameter when the BBMC speed regulation system operates stably through numerical simulation according to the discrete iterative mapping model of the BBMC speed regulation system based on finite time control obtained in the step (7), specifically:
step (8-1): setting system parameters, including: BBMC main circuit inductance L, main circuit capacitance C, motor stator self-inductance LsSelf-inductance of rotor LrStator and rotor mutual inductance LmStator resistor RsRotor resistance RrN number of pole pairspMoment of inertia J, load torque TLSwitching period T and maximum iteration number N of power switching tubemaxThree-phase capacitor reference voltage in BBMCPeriod T of reference voltage of capacitor0(and satisfy T)0kT, k being a positive integer), a control parameter k1、k2、k3、α1、α2and alpha3And increasing the change parameter delta X, and taking 1 as the initial value of the counting variable q.
Step (8-2): first order the controlled variable Y1~Y6Respectively representing control parameters k1、k2、k3、α1、α2and alpha3;
Step (8-3): let control variable YqThe initial value of the variable parameter X is set to be 0, and other control variables are kept unchanged;
step (8-4): calculating the state variable i at the time (n +1) T by equation (15)α(n+1)And iβ(n+1);
Step (8-5): judgment of iα(n+1)And iβ(n+1)Whether or not to satisfy i simultaneouslyα(n+1)=iα(n)And iβ(n+1)=iβ(n)If yes, the system is in a stable state, and the step (8-8) is executed; otherwise, executing the step (8-6);
step (8-6): judging whether the iteration number N is more than NmaxIf yes, executing the step (8-7); otherwise, adding 1 to the iteration number n, and returning to the step (8-4);
step (8-7): adding delta X to the change parameter X, enabling the iteration number n to return to 1, and returning to the step (8-4);
step (8-8): let the corresponding change parameter value X at this time be the lower limit value of the parameter stability domain, that is: xmin=X;
Step (8-9): adding delta X to the variation parameter X, and enabling the iteration number n to return to 1;
step (8-10): calculating the state variable i at the time (n +1) T by equation (15)α(n+1)And iβ(n+1);
Step (8-11): judgment of iα(n+1)And iβ(n+1)Whether or not to satisfy i simultaneouslyα(n+1)=iα(n)And iβ(n+1)=iβ(n)If yes, the system is in a stable state, and the step (8-12) is executed; otherwise, executing the step (8-13);
step (8-12): let variable XmaxX, then returning to step (8-9);
step (8-13): judging whether the iteration number N is more than NmaxIf yes, executing the step (8-14); otherwise, adding 1 to the iteration number n, and returning to the step (8-10);
step (8-14): let Yqmin=Xmin,Yqmax=XmaxObtaining the control variable Y of the speed regulating system when the speed regulating system operates stablyqHas a value range of (Y)qmin,Yqmax);
Step (8-15): judging whether the counting variable q is smaller than 6, if so, executing the step (8-16), otherwise, executing the step (8-17);
step (8-16): adding 1 to the counting variable q, and returning to the step (8-3);
step (8-17): the resulting controlled variable YqValue range (Y)qmin,Yqmax) (q 1, 2.. 6) corresponds to the control parameter k, respectively1、k2、k3、α1、α2and alpha3The specific value of each control parameter is determined according to the value range of each control parameter, and the stable operation of the BBMC speed regulating system based on the limited time control can be realized.
Compared with the prior art, the invention has the technical effects that: the method takes inductive current, capacitor voltage and output current in a Buck-Boost matrix converter (BBMC) as state variables, and establishes a state differential equation of an asynchronous motor speed regulating system based on the BBMC; obtaining a dynamic equation of the system according to the state differential equation; designing a control function of the BBMC speed regulating system according to the dynamic equation of the system; obtaining the duty ratio of a power switch tube in the BBMC according to the control function and the finite time control principle of the BBMC speed regulating system; obtaining a discrete iteration mapping model of a BBMC inverter stage according to the duty ratio of the power switch tube in the BBMC; establishing a discrete iteration mapping model of the three-phase asynchronous motor; obtaining a discrete iteration mapping model of the BBMC speed regulating system based on finite time control according to the discrete iteration mapping model of the BBMC inverter stage and the discrete iteration mapping model of the three-phase asynchronous motor; and obtaining the value range of the relevant control parameters when the BBMC speed regulating system based on the finite time control realizes stable operation through numerical simulation according to the discrete iterative mapping model of the BBMC speed regulating system based on the finite time control. The invention adopts a finite time control strategy aiming at an asynchronous motor speed regulating system taking BBMC as a power converter, determines the stable domain range of control parameters of the system on the premise of ensuring the stable operation of the system, and can ensure the stable operation of the BBMC speed regulating system.
Drawings
FIG. 1 is a topological structure diagram of a BBMC-based asynchronous motor speed regulation system of the present invention;
FIG. 2 is a flow chart of the present invention;
FIG. 3 is a flow chart of determining the value range of the control parameter of the BBMC speed control system discrete iterative mapping model based on finite time control in the present invention.
Detailed Description
Fig. 1 is a topology structure diagram of an asynchronous motor speed regulation system based on BBMC according to an embodiment of the present invention. The BBMC adopts the structural form of an AC-DC-AC two-stage converter, the rectification stage of the BBMC is a 3/2 phase matrix converter, and the inverter stage adopts the structural form of a three-phase Buck-Boost inverter and consists of three Buck-Boost DC/DC converters with the same structure; and three-phase stator windings of the three-phase asynchronous motor are respectively connected to three output ends of the BBMC.
Referring to fig. 2, a flow chart of a method for determining a control parameter stability region of a finite time control BBMC speed control system provided by the invention is shown. The method comprises the following steps:
step (1): the method comprises the following steps of establishing a state differential equation of an asynchronous motor speed regulating system based on BBMC by taking inductive current, capacitor voltage and output current in BBMC as state variables, wherein the state differential equation specifically comprises the following steps:
wherein: i.e. iL、uCAnd ioRespectively representing inductive current, capacitor voltage and output current in BBMC, E is input voltage of BBMC inverter stage, u is output voltage of BBMC inverter stageDFor the common terminal voltage of three-phase stator winding of asynchronous motor, L is BBMC main circuit inductance, C is BBMC main circuit capacitance, R is1And L1Respectively is the equivalent resistance and the equivalent inductance of the single-phase winding of the asynchronous motor, d is the duty ratio of a power switch in the BBMC, and d belongs to [0,1 ]]。
Step (2): obtaining a dynamic equation of the system according to the state differential equation obtained in the step (1), specifically:
let the reference output voltage of BBMC be urefThe error of the output voltage of the BBMC is obtained as follows:
x1=uC-uref(2)
according to the formula (1) and the formula (2), the dynamic equation of the BBMC speed regulating system is obtained as follows:
wherein: x is the number of1、x2、x3、iL、uCAnd ioAre all variables of time t.
And (3): designing a control function of the BBMC speed regulating system according to the dynamic equation of the BBMC speed regulating system obtained in the step (2), which specifically comprises the following steps:
step (3-1): defining a time scale coordinate transformation as:
wherein: and t is Ks, and K is a coordinate transformation coefficient.
Step (3-2): and (3) according to the time scale coordinate transformation formula shown in the formula (4), transforming the formula (3) to obtain:
step (3-3): the control function of the BBMC governor system obtained according to equation (5) is:
and (4): and (4) obtaining the duty ratio of a power switching tube in the BBMC according to the control function of the BBMC speed regulating system obtained in the step (3) and the finite time control principle, specifically:
firstly, according to the finite time control principle, a system control function f(s) is determined as follows:
wherein: satαAs a function of saturation, k1、k2、k3、α1、α2and alpha3Are control parameters.
Then, it is obtained from equations (6) and (7):
finally, according to the formula (4) and the formula (8), the duty ratio of the power switch tube in the BBMC is obtained as follows:
and (5): obtaining a discrete iteration mapping model of the BBMC inverter stage according to the duty ratio of the power switching tube in the BBMC obtained in the step (4), specifically:
the method comprises the following steps of establishing a differential equation of state by taking inductive current and capacitor voltage in BBMC as state variables:
according to the formula (9), the on-time T of the power switch tube in the nth switching period T can be obtainedonAnd off time toffRespectively is as follows:
wherein: dnAnd represents the duty ratio of the power switch tube in the BBMC in the nth switching period T.
According to equations (10) and (11), the discrete iterative mapping model for the BBMC inversion stage is obtained as follows:
wherein: u. ofC(n+1)And iL(n+1)Respectively representing the capacitor voltage and the inductor current i in the (n +1) T moment BBMCo(n)Denotes the output current of BBMC at time nT, M, N and ω are intermediate variables, anduC(n)and iL(n)Respectively representing the capacitor voltage and the inductor current at time nT BBMC.
And (6): establishing a discrete iteration mapping model of the three-phase asynchronous motor, specifically:
firstly, establishing a state differential equation of the three-phase asynchronous motor on a two-phase static coordinate system as follows:
wherein: u. ofsαAnd usβRepresenting the stator voltage, i, of the motor in two-phase stationary coordinates, respectivelysαAnd isβRepresenting the stator currents of the motor in two stationary phases, ΨrαAnd ΨrβRespectively representing the rotor flux, T, of the motor in two-phase stationary coordinatesLRepresenting motor load torque, ωrRepresenting the angular velocity, L, of the rotor of the machines、Lr、Lm、Rs、Rr、npAnd J respectively represent the self inductance of the stator, the self inductance of the rotor, the mutual inductance of the stator and the rotor, the resistance of the stator, the resistance of the rotor, the pole pair number and the moment of inertia of the motor,representing the leakage coefficient, T, of the motorr=Lrthe/R represents the rotor electromagnetic time constant.
Then, the differential equation of state shown in equation (13) is discretized by the Runge-Kutta method to obtain:
wherein:is the state vector at time (n +1) T,is the state vector at time nT, K1、K2、K3And K4Are all intermediate variables, and K1=f(xn,yn),K4=f(xn+TK0,yn+TK3),
usα(n)And usβ(n)Respectively representing the stator voltage i of the motor at the moment nT on the two-phase stationary coordinatesα(n)And isβ(n)Representing the stator currents of the machine at the moment nT in two-phase stationary coordinates, Ψrα(n)And Ψrβ(n)Respectively representing the rotor flux linkage, omega, at time nT on the two-phase stationary coordinater(n)Representing the angular velocity of the rotor of the machine at time nT.
And (7): obtaining a discrete iteration mapping model of the BBMC speed regulating system based on finite time control according to the discrete iteration mapping model of the BBMC inverter stage obtained in the step (5) and the discrete iteration mapping model of the three-phase asynchronous motor obtained in the step (6), and specifically:
wherein:andrespectively representing the three-phase capacitor voltages at time (n +1) T BBMC,andrespectively represents three-phase inductive current in the BBMC at the (n +1) T moment, andrespectively represents the conduction time of the three-phase power switch tube in the n-th switch period T,andrespectively showing the turn-off time of the three-phase power switch tube in the BBMC in the nth switching period T,andrespectively representing the duty ratio of three-phase power switching tubes in the BBMC in the nth switching period T, MA、MB、MC、NA、NBAnd NCAre all intermediate variables, and
and (8): and (4) obtaining the value range of the control parameter of the BBMC speed regulation system in stable operation through numerical simulation according to the discrete iterative mapping model of the BBMC speed regulation system based on the finite time control obtained in the step (7). Referring to fig. 3, a flowchart for determining a control parameter value range according to a finite time control-based BBMC speed regulation system discrete iteration mapping model according to an embodiment of the present invention includes the following specific steps:
step (8-1): setting system parameters, including: BBMC main circuit inductance L, main circuit capacitance C, motor stator self-inductance LsSelf-inductance of rotor LrStator and rotor mutual inductance LmStator resistor RsRotor resistance RrN number of pole pairspMoment of inertia J, load torque TLSwitching period T and maximum iteration number N of power switching tubemaxThree-phase capacitor reference voltage in BBMCPeriod T of reference voltage of capacitor0(and satisfy T)0kT, k being a positive integer), a control parameter k1、k2、k3、α1、α2and alpha3And increasing the change parameter delta X, and taking 1 as the initial value of the counting variable q.
Step (8-2): first order the controlled variable Y1~Y6Respectively representing control parameters k1、k2、k3、α1、α2and alpha3;
Step (8-3): let control variable YqThe initial value of the variable parameter X is set to be 0, and other control variables are kept unchanged;
step (8-4): calculating the state variable i at the time (n +1) T by equation (15)α(n+1)And iβ(n+1);
Step (8-5): judgment of iα(n+1)And iβ(n+1)Whether or not to satisfy i simultaneouslyα(n+1)=iα(n)And iβ(n+1)=iβ(n)If yes, the system is in a stable state, and the step (8-8) is executed; otherwise, executing the step (8-6);
step (8-6): judging whether the iteration number N is more than NmaxIf yes, executing the step (8-7); otherwise, adding 1 to the iteration number n, and returning to the step (8-4);
step (8-7): adding delta X to the change parameter X, enabling the iteration number n to return to 1, and returning to the step (8-4);
step (8-8): let the corresponding change parameter value X at this time be the lower limit value of the parameter stability domain, that is: xmin=X;
Step (8-9): adding delta X to the variation parameter X, and enabling the iteration number n to return to 1;
step (8-10): calculating the state variable i at the time (n +1) T by equation (15)α(n+1)And iβ(n+1);
Step (8-11): judgment of iα(n+1)And iβ(n+1)Whether or not to satisfy i simultaneouslyα(n+1)=iα(n)And iβ(n+1)=iβ(n)If yes, the system is in a stable state, and the step (8-12) is executed; otherwise, executing the step (8-13);
step (8-12): let an intermediate variable XmaxX, then returning to step (8-9);
step (8-13): judging whether the iteration number N is more than NmaxIf yes, executing the step (8-14); otherwise, adding 1 to the iteration number n, and returning to the step (8-10);
step (8-14): let Yqmin=Xmin,Yqmax=XmaxObtaining the control variable Y of the speed regulating system when the speed regulating system operates stablyqHas a value range of (Y)qmin,Yqmax);
Step (8-15): judging whether the counting variable q is smaller than 6, if so, executing the step (8-16), otherwise, executing the step (8-17);
step (8-16): adding 1 to the counting variable q, and returning to the step (8-3);
step (8-17): the resulting controlled variable YqValue range (Y)qmin,Yqmax) (q 1, 2.. 6) corresponds to the control parameter k, respectively1、k2、k3、α1、α2and alpha3The specific value of each control parameter is determined according to the value range of each control parameter, and the stable operation of the BBMC speed regulating system based on the limited time control can be realized.
Claims (9)
1. A method for determining a control parameter stability region of a finite time control BBMC speed regulation system is characterized by comprising the following steps:
step (1): establishing a state differential equation of the BBMC speed regulation system by taking inductive current, capacitor voltage and output current in the BBMC as state variables;
step (2): obtaining a dynamic equation of the BBMC speed regulating system according to the state differential equation obtained in the step (1);
and (3): designing a control function of the BBMC speed regulating system according to the dynamic equation of the BBMC speed regulating system obtained in the step (2);
and (4): obtaining the duty ratio of a power switching tube in the BBMC according to the control function of the BBMC speed regulation system obtained in the step (3) and a finite time control principle;
and (5): obtaining a discrete iteration mapping model of the BBMC inverter stage according to the duty ratio of the power switch tube in the BBMC obtained in the step (4);
and (6): establishing a discrete iteration mapping model of the three-phase asynchronous motor;
and (7): obtaining a discrete iteration mapping model of the BBMC speed regulating system based on finite time control according to the discrete iteration mapping model of the BBMC inverter stage obtained in the step (5) and the discrete iteration mapping model of the three-phase asynchronous motor obtained in the step (6);
and (8): and (4) obtaining the value range of the control parameter when the BBMC speed regulation system stably operates through numerical simulation according to the discrete iterative mapping model of the BBMC speed regulation system based on the finite time control obtained in the step (7).
2. The finite time control BBMC speed regulating system control parameter stability domain determining method as claimed in claim 1, wherein: the specific steps of the step (1) are as follows:
wherein: i.e. iL、uCAnd ioRespectively representing inductive current, capacitor voltage and output current in BBMC, E is input voltage of BBMC inverter stage, u is output voltage of BBMC inverter stageDFor the common terminal voltage of three-phase stator winding of asynchronous motor, L is BBMC main circuit inductance, C is BBMC main circuit capacitance, R is1And L1Respectively is the equivalent resistance and the equivalent inductance of the single-phase winding of the asynchronous motor, d is the duty ratio of a power switch in the BBMC, and d belongs to [0,1 ]]。
3. The finite time control BBMC speed regulating system control parameter stability domain determining method as claimed in claim 1, wherein: the specific steps of the step (2) are as follows:
let the reference output voltage of BBMC be urefThe error of the output voltage of the BBMC is obtained as follows:
x1=uC-uref(2)
according to the formula (1) and the formula (2), the dynamic equation of the BBMC speed regulating system is obtained as follows:
wherein: x is the number of1、x2、x3、iL、uCAnd ioAre all variables of time t.
4. The finite time control BBMC speed regulating system control parameter stability domain determining method of claim 3, characterized in that: the specific steps of the step (3) are as follows:
step (3-1): defining a time scale coordinate transformation as:
wherein: k is a coordinate transformation coefficient;
step (3-2): and (3) according to the time scale coordinate transformation formula shown in the formula (4), transforming the formula (3) to obtain:
step (3-3): the control function of the BBMC governor system obtained according to equation (5) is:
5. the finite time control BBMC speed regulating system control parameter stability domain determining method of claim 4, wherein: the specific steps of the step (4) are as follows:
according to the finite time control principle, determining a system control function f(s) as follows:
wherein: satαAs a function of saturation, k1、k2、k3、α1、α2and alpha3Is a control parameter;
according to the formulae (6) and (7):
finally, according to equations (4) and (8), the duty ratio of the power switch tube in the BBMC can be obtained as follows:
6. the finite time control BBMC speed regulating system control parameter stability domain determining method of claim 5, characterized in that: the specific steps of the step (5) are as follows:
the method comprises the following steps of establishing a differential equation of state by taking inductive current and capacitor voltage in BBMC as state variables:
according to the formula (9), the on-time T of the power switch tube in the nth switching period T can be obtainedonAnd off time toffRespectively is as follows:
wherein: dnThe duty ratio of the power switch tube in the BBMC in the nth switching period T is represented;
according to the equations (10) and (11), the discrete iterative mapping model of the BBMC inverter stage is obtained as follows:
wherein: u. ofC(n+1)And iL(n+1)Respectively representing the capacitor voltage and the inductor current i in the (n +1) T moment BBMCo(n)Denotes the output current of BBMC at time nT, M, N and ω are intermediate variables, anduC(n)and iL(n)Respectively representing the capacitor voltage and the inductor current at time nT BBMC.
7. The finite time control BBMC speed regulating system control parameter stability domain determining method as claimed in claim 1, wherein: the specific steps of the step (6) are as follows:
establishing a state differential equation of the three-phase asynchronous motor on a two-phase static coordinate system as follows:
wherein: u. ofsαAnd usβRepresenting the stator voltage, i, of the motor in two-phase stationary coordinates, respectivelysαAnd isβRepresenting the stator currents of the motor in two stationary phases, ΨrαAnd ΨrβRespectively representing the rotor flux, T, of the motor in two-phase stationary coordinatesLRepresenting motor load torque, ωrRepresenting the angular velocity, L, of the rotor of the machines、Lr、Lm、Rs、Rr、npAnd J respectively represents self-inductance of the stator, self-inductance of the rotor and fixed rotation of the motorSub-mutual inductance, stator resistance, rotor resistance, pole pair number and moment of inertia,representing the leakage coefficient, T, of the motorr=Lrthe/R represents the rotor electromagnetic time constant;
discretization of formula (13) by Runge-Kutta method gives:
wherein:is the state vector at time (n +1) T,is the state vector at time nT, K1、K2、K3And K4Are all intermediate variables, and K1=f(xn,yn),K4=f(xn+TK0,yn+TK3),
usα(n)And usβ(n)Respectively representing the stator voltage i of the motor at the moment nT on the two-phase stationary coordinatesα(n)And isβ(n)Representing the stator currents of the machine at the moment nT in two-phase stationary coordinates, Ψrα(n)And Ψrβ(n)Respectively representing the rotor flux linkage, omega, at time nT on the two-phase stationary coordinater(n)Representing the angular velocity of the rotor of the machine at time nT.
8. The finite time control BBMC speed regulating system control parameter stability domain determining method of claim 7, wherein: the specific steps of the step (7) are as follows:
wherein:andrespectively representing the three-phase capacitor voltages at time (n +1) T BBMC,andrespectively represents three-phase inductive current in the BBMC at the (n +1) T moment, andrespectively represents the conduction time of the three-phase power switch tube in the n-th switch period T,andrespectively showing the turn-off time of the three-phase power switch tube in the BBMC in the nth switching period T,andrespectively representing the duty ratio of three-phase power switching tubes in the BBMC in the nth switching period T, MA、MB、MC、NA、NBAnd NCAre all intermediate variables, and
9. the finite time control BBMC speed regulating system control parameter stability domain determining method as claimed in claim 1, wherein: the specific steps of the step (8) are as follows:
step (8-1): setting system parameters, including: BBMC main circuit inductance L, main circuit capacitance C, motor stator self-inductance LsSelf-inductance of rotor LrStator and rotor mutual inductance LmStator resistor RsRotor resistance RrN number of pole pairspMoment of inertia J, load torque TLSwitching period T and maximum iteration number N of power switching tubemaxThree-phase capacitor reference voltage in BBMCPeriod T of reference voltage of capacitor0And satisfy T0K is positive integer, control parameter k1、k2、k3、α1、α2and alpha3Changing parameter increment delta X, and taking 1 as the initial value of a counting variable q;
step (8-2): first order the controlled variable Y1~Y6Respectively representing control parameters k1、k2、k3、α1、α2and alpha3;
Step (8-3): let control variable YqThe initial value of the variable parameter X is set to be 0, and other control variables are kept unchanged;
step (8-4): calculating the state variable i at the time (n +1) T by equation (15)α(n+1)And iβ(n+1);
Step (8-5): judgment of iα(n+1)And iβ(n+1)Whether or not to satisfy i simultaneouslyα(n+1)=iα(n)And iβ(n+1)=iβ(n)If yes, the system is in a stable state, and the step (8-8) is executed; otherwise, executing the step (8-6);
step (8-6): judging whether the iteration number N is more than NmaxIf yes, executing the step (8-7); otherwise, adding 1 to the iteration number n, and returning to the step (8-4);
step (8-7): adding delta X to the change parameter X, enabling the iteration number n to return to 1, and returning to the step (8-4);
step (8-8): let the corresponding change parameter value X at this time be the lower limit value of the parameter stability domain, that is: xmin=X;
Step (8-9): adding delta X to the variation parameter X, and enabling the iteration number n to return to 1;
step (8-10): calculating the state variable i at the time (n +1) T by equation (15)α(n+1)And iβ(n+1);
Step (8-11): judgment of iα(n+1)And iβ(n+1)Whether or not to satisfy i simultaneouslyα(n+1)=iα(n)And iβ(n+1)=iβ(n)If yes, the system is in a stable state, and the step (8-12) is executed; otherwise, executing the step (8-13);
step (8-12): let variable XmaxX, then returning to step (8-9);
step (8-13): judging whether the iteration number N is more than NmaxIf yes, executing the step (8-14); otherwise, adding 1 to the iteration number n, and returning to the step (8-10);
step (8-14): let Yqmin=Xmin,Yqmax=XmaxObtaining the control variable Y of the speed regulating system when the speed regulating system operates stablyqHas a value range of (Y)qmin,Yqmax);
Step (8-15): judging whether the counting variable q is smaller than 6, if so, executing the step (8-16), otherwise, executing the step (8-17);
step (8-16): adding 1 to the counting variable q, and returning to the step (8-3);
step (8-17): the resulting controlled variable YqValue range (Y)qmin,Yqmax) Q is 1,2, 6, i.e. corresponding to the control parameter k, respectively1、k2、k3、α1、α2and alpha3And determining the specific value of each control parameter according to the value range of each control parameter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010064506.4A CN111181468B (en) | 2020-01-20 | 2020-01-20 | Method for determining control parameter stability domain of finite time control BBMC speed regulation system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010064506.4A CN111181468B (en) | 2020-01-20 | 2020-01-20 | Method for determining control parameter stability domain of finite time control BBMC speed regulation system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111181468A true CN111181468A (en) | 2020-05-19 |
CN111181468B CN111181468B (en) | 2022-01-04 |
Family
ID=70656422
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010064506.4A Active CN111181468B (en) | 2020-01-20 | 2020-01-20 | Method for determining control parameter stability domain of finite time control BBMC speed regulation system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111181468B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116522580A (en) * | 2023-02-22 | 2023-08-01 | 广东轻工职业技术学院 | Buck-Boost intermediate frequency inversion main circuit parameter optimization method |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101780798A (en) * | 2009-01-19 | 2010-07-21 | 比亚迪股份有限公司 | Dual-clutch gear-shifting control method and device thereof |
WO2016049135A1 (en) * | 2014-09-25 | 2016-03-31 | Bae Systems Controls Inc. | Balanced bidirectional buck-boost converters and associated systems and methods |
CN106788086A (en) * | 2017-02-20 | 2017-05-31 | 青岛大学 | Consider the asynchronous machine command filtering finite time fuzzy control method of input saturation |
CN106788046A (en) * | 2017-02-20 | 2017-05-31 | 青岛大学 | Permagnetic synchronous motor command filtering finite time fuzzy control method |
CN108809176A (en) * | 2018-06-22 | 2018-11-13 | 湖南科技大学 | A kind of asynchronous motor speed-regulating system control method based on Buck-Boost matrix converters |
CN109842344A (en) * | 2019-03-07 | 2019-06-04 | 湖南科技大学 | BBMC asynchronous motor speed-regulating system control parameter self-adapting regulation method |
CN110690842A (en) * | 2019-10-31 | 2020-01-14 | 湖南科技大学 | Method for determining main circuit parameter stability region of three-phase asynchronous motor speed regulating system |
-
2020
- 2020-01-20 CN CN202010064506.4A patent/CN111181468B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101780798A (en) * | 2009-01-19 | 2010-07-21 | 比亚迪股份有限公司 | Dual-clutch gear-shifting control method and device thereof |
WO2016049135A1 (en) * | 2014-09-25 | 2016-03-31 | Bae Systems Controls Inc. | Balanced bidirectional buck-boost converters and associated systems and methods |
CN106788086A (en) * | 2017-02-20 | 2017-05-31 | 青岛大学 | Consider the asynchronous machine command filtering finite time fuzzy control method of input saturation |
CN106788046A (en) * | 2017-02-20 | 2017-05-31 | 青岛大学 | Permagnetic synchronous motor command filtering finite time fuzzy control method |
CN108809176A (en) * | 2018-06-22 | 2018-11-13 | 湖南科技大学 | A kind of asynchronous motor speed-regulating system control method based on Buck-Boost matrix converters |
CN109842344A (en) * | 2019-03-07 | 2019-06-04 | 湖南科技大学 | BBMC asynchronous motor speed-regulating system control parameter self-adapting regulation method |
CN110690842A (en) * | 2019-10-31 | 2020-01-14 | 湖南科技大学 | Method for determining main circuit parameter stability region of three-phase asynchronous motor speed regulating system |
Non-Patent Citations (1)
Title |
---|
刘继 等: ""基于FTC的BBMC调速控制策略及参数优化"", 《自动化学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116522580A (en) * | 2023-02-22 | 2023-08-01 | 广东轻工职业技术学院 | Buck-Boost intermediate frequency inversion main circuit parameter optimization method |
CN116522580B (en) * | 2023-02-22 | 2023-12-29 | 广东轻工职业技术学院 | Buck-Boost intermediate frequency inversion main circuit parameter optimization method |
Also Published As
Publication number | Publication date |
---|---|
CN111181468B (en) | 2022-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104518678B (en) | High dynamic control apparatus for current source converter | |
Xiao et al. | An improved model predictive control scheme for the PWM rectifier-inverter system based on power-balancing mechanism | |
CN105577058A (en) | Novel fuzzy active disturbance rejection controller based five-phase fault-tolerant permanent magnet motor speed control method | |
CN110165953B (en) | PMSM speed regulation control method based on approximation law | |
CN109167538B (en) | Control method for inhibiting torque ripple of permanent magnet brushless direct current motor based on double-layer structure | |
CN111800056A (en) | Permanent magnet synchronous motor three-vector model predicted torque control method based on novel switch table | |
CN111181468B (en) | Method for determining control parameter stability domain of finite time control BBMC speed regulation system | |
CN107659231A (en) | A kind of ultrahigh speed permagnetic synchronous motor method for controlling number of revolution based on the switching of single electric current field weakening mode | |
CN107040178A (en) | Two-phase hybrid stepping motor closed-loop control device and method | |
CN110690842B (en) | Method for determining main circuit parameter stability region of three-phase asynchronous motor speed regulating system | |
CN102427301B (en) | Control method for three-phase pulse width modulation (PWM) rectifier | |
Chiang et al. | Integral variable structure controller with grey prediction for synchronous reluctance motor drive | |
CN113131830A (en) | Efficiency optimization control method, system, medium and processing terminal for linear induction motor | |
CN112003318A (en) | Wind power grid-connected inverter direct-current bus voltage control method | |
Ye et al. | Full discrete sliding mode controller for three phase PWM rectifier based on load current estimation | |
He et al. | Modeling and stability analysis of three-phase PWM rectifier | |
Liaw et al. | A fuzzy adapted field-oriented mechanism for induction motor drive | |
CN113992108A (en) | Direct torque control system and method of switched reluctance motor | |
CN109768750B (en) | Zero-sequence current suppression method of open-winding permanent magnet motor based on extended virtual three phases | |
CN107370402A (en) | A kind of method for handover control based on discrete-time Lyapunov function | |
CN110323813B (en) | Electric vehicle charging pile pre-stage PWM rectifier control method based on performance protection control | |
Sun et al. | Novel modulation method for torque ripple suppression of brushless DC motors based on SIMO DC–DC converter | |
Talavat et al. | Direct predictive control of asynchronous machine torque using matrix converter | |
Alharbi et al. | Adaptive backstepping control of pwm dc-dc boost power converter considering parasitic elements effect | |
Kumar et al. | Energy efficient drive system for domestic and agriculture applications: a comparative study of SPIM and SRM drives |
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 |