CN111522241A - Active fault-tolerant control method and device based on fixed time observer - Google Patents
Active fault-tolerant control method and device based on fixed time observer Download PDFInfo
- Publication number
- CN111522241A CN111522241A CN202010381563.5A CN202010381563A CN111522241A CN 111522241 A CN111522241 A CN 111522241A CN 202010381563 A CN202010381563 A CN 202010381563A CN 111522241 A CN111522241 A CN 111522241A
- Authority
- CN
- China
- Prior art keywords
- representing
- control system
- observer
- fault
- rotorcraft
- 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
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Abstract
The invention discloses an active fault-tolerant control method and device based on a fixed time observer, which are applied to the fault-tolerant control of a rotor type aircraft, and the method comprises the following steps: establishing a mathematical model of a rotary wing type aircraft attitude control system; establishing a fixed-time sliding-mode observer according to a mathematical model of a rotary wing type aircraft attitude control system; establishing a fault-tolerant controller according to a mathematical model of a rotary wing type aircraft attitude control system and a fixed time sliding mode observer; controlling the rotor type aircraft by using a fault-tolerant controller; the invention has the advantages that: the system stability is strong, and the modeling error is small.
Description
Technical Field
The invention relates to the technical field of fault-tolerant control of rotary wing type aircrafts, in particular to an active fault-tolerant control method and device based on a fixed time observer.
Background
Compared with a fixed-wing aircraft, the rotary-wing aircraft has the advantages of being capable of vertically taking off and landing and hovering in the air while having a high flying speed. In recent years, as the related technology of the unmanned rotary wing type aircraft is mature day by day, the unmanned rotary wing type aircraft is widely applied in various industries. For example, military monitoring and investigation tasks, rescue tasks in fire and earthquake, etc.; meanwhile, the method is also one of important tools and means for Mars detection.
The rotary wing type aircraft can efficiently and accurately complete tasks and is highly dependent on an internal flight control system. However, due to the complex working environment, the external atmospheric conditions and the activities of other objects can generate certain disturbances on the aircraft, which affects the control performance of the aircraft. Meanwhile, the high-intensity operation inevitably causes abrasion of machine body parts and damage of an actuator and a sensor. Therefore, the research on fault-tolerant control methods for rotor-type aircrafts is particularly important to ensure the efficient completion of the control tasks thereof.
Chinese patent publication No. CN110502027A discloses a four-rotor unmanned aerial vehicle attitude fault-tolerant control method based on a self-adaptive terminal sliding mode, and in order to solve a time-varying fault, firstly, mathematical models of a position subsystem and an attitude subsystem of a four-rotor unmanned aerial vehicle are established, and an error function of the position subsystem is defined. Then, under the condition that the attitude subsystem pitch angle takes place the time-varying trouble, application unmanned aerial vehicle attitude angle error is calmed to zero, realizes the fault-tolerant of self-adaptation, but this patent application system control's stability does not solve.
Chinese patent publication No. CN105353615B discloses an active fault-tolerant control method of a four-rotor aircraft based on a sliding-mode observer, which designs the sliding-mode observer, performs linear transformation on a system, reconstructs an actuator fault based on the idea of equivalent error injection, adds compensation control in sliding-mode control by utilizing a reconstructed estimated value of the actuator fault, and finally forms a complete active fault-tolerant controller. According to the method, the sliding-mode observer is designed to reconstruct and estimate the fault, the online adjustment of the gain of the controller can be realized, the control law is optimal, the control precision and the response speed of the flight of the four-rotor aircraft are effectively improved, and the design basis of the fault-tolerant controller can be provided for the complex four-rotor aircraft with the actuator fault. The invention can realize the online adjustment of the gain of the controller, so that the provided control law is optimal. However, the actual mathematical model of the four-rotor aircraft is highly nonlinear and strongly coupled, and the fault-tolerant control scheme adopts a linear model, so that the modeling error is large; in addition, the problem of jitter of sliding mode control is not solved, so that the system stability is poor.
Disclosure of Invention
The invention aims to solve the technical problems of poor system stability and large modeling error of the control method of the rotor type aircraft in the prior art.
The invention solves the technical problems through the following technical means: an active fault-tolerant control method based on a fixed time observer is applied to fault-tolerant control of a rotor type aircraft, and comprises the following steps:
the method comprises the following steps: establishing a mathematical model of the attitude control system of the rotary wing aircraft according to the control input dynamic value of the control system, the nonlinear dynamic value of the control system and the sum of the comprehensive disturbance quantities;
step two: establishing a fixed-time sliding mode observer according to an input dynamic value, a nonlinear dynamic value, an error variable and a state quantity of the sliding mode observer in a mathematical model of the attitude control system of the rotor type aircraft;
step three: establishing a fault-tolerant controller according to a mathematical model of a rotary wing type aircraft attitude control system, a fixed time sliding mode observer and a tracking error;
step four: and controlling the rotor type aircraft by using a fault-tolerant controller.
The method comprises the steps of obtaining system fault information by using a fixed time sliding mode observer, designing a mathematical model of a rotor type aircraft attitude control system containing comprehensive disturbance, establishing a fault-tolerant controller according to the mathematical model of the rotor type aircraft attitude control system and the fixed time sliding mode observer, controlling the rotor type aircraft by using the fault-tolerant controller, considering factors such as tracking error and the comprehensive disturbance by using the fault-tolerant controller, having small modeling error, effectively processing faults of the system and having stability.
Preferably, the first step includes: using formulas
Establishing a mathematical model of a rotary wing type aircraft attitude control system;
wherein X represents a state variable of the control system; x1A first state quantity representing the control system,representing the first derivative of a first state quantity of the control system, and X1=[α,β]T,[]TRepresenting the transpose of the matrix, α representing the yaw angle of the rotorcraft, β representing the pitch angle of the rotorcraft;
X2a second state quantity representing the control system,represents a first derivative of a second state quantity of the control system, andrepresenting the first derivative of the yaw angle of the rotorcraft,representing the first derivative of the pitch angle of the rotorcraft;
f (X) represents a nonlinear dynamic value of the control system, andm represents the effective mass of the aircraft, g represents the gravitational acceleration constant, LaRepresenting the length of the center of the rotorcraft from the front motor;
g (X) represents a control input dynamic value of the control system, andKfrepresenting the coefficient of moment generated by the motor, JαRepresenting the moment of inertia of the yaw axis, JβIndicating the pitch axisRotational inertia of LhRepresenting the length of the rear motor from the center of the rotorcraft;
u represents a matrix form of an input amount of the control system, andVfrepresenting the voltage, V, applied to the front motor of a rotorcraftbRepresenting the voltage applied to the rear motor of the rotorcraft;
Δ represents the integrated disturbance amount, and Δ ═ Δ f (x)) - (g (x)) + Δ g (x)) (ρi-1)U(t)+D(t),ρiFailure coefficient and p representing occurrence of actuator failure of rotorcrafti∈ [0, 1), Δ f (x) represents a first uncertainty value of the control system, Δ g (x) represents a second uncertainty value of the control system, u (t) represents a functional form of an input quantity of the control system, d (t) represents a functional form of an unknown external disturbance.
Preferably, the second step includes: using formulas
Establishing a fixed time sliding-mode observer, wherein ξ1Representing a first quantity of state of the sliding-mode observer, ξ2A first state quantity of the sliding mode observer is represented,a first derivative of a first state quantity representing a sliding mode observer,denotes the first derivative of the first state quantity of the sliding-mode observer, e denotes the error variable, and e is ξ1-X1And | | represents a euclidean distance symbol, p is a power exponent and p > 1, sign () represents a sign function, λ1Is a preset observer firstParameter, λ2Is a predetermined observer second parameter, λ3Is a predetermined observer third parameter, and λ1,λ2,λ3Are all constants greater than zero.
Preferably, the third step includes: using formulas
A fault tolerant controller is established in which, among other things,is an estimate of the overall disturbance variable Delta, X2dTo control a second quantity of state X of the system2The reference signal of (a) is set,for virtual controlAn output through a differentiator;is a preset second controller parameter, E2Is a second tracking error, and E2=X2-X2d,Is a first auxiliary function.
Preferably, the third step further comprises: using formulas
Obtaining a functionWherein z represents the state variable of the differentiator, z1Expressed as a first state quantity of the differentiator, z2A second state quantity representing a differentiator, s (t) being an input to the differentiator and representing a perturbation parameterAnd > 0;
using formulas
using formulas
zout=z2
A differentiator is established in which, among other things,is the first derivative of the first state quantity of the differentiator,is the first derivative of the second state quantity of the differentiator, zoutIs the output of the differentiator.
Preferably, the third step further comprises: using formulas
Constructing a first auxiliary function, wherein E1Is a first tracking error, and E1=x1-x1d,X1dFor controlling a first quantity of state X of the system1The reference signal of (a);
is an auxiliary variable, andis the first derivative of the auxiliary variable,is a preset first controller parameter, wherein,for controlling a first quantity of state X of the system1The first derivative of the reference signal.
The invention also provides an active fault-tolerant control device based on the fixed time observer, which is applied to the fault-tolerant control of the rotor type aircraft, and the device comprises:
the model establishing module is used for establishing a mathematical model of the attitude control system of the rotary wing aircraft according to the control input dynamic value of the control system, the nonlinear dynamic value of the control system and the sum of the comprehensive disturbance quantities;
the sliding mode observer establishing module is used for establishing a fixed time sliding mode observer according to an input dynamic value, a nonlinear dynamic value, an error variable and a state quantity of the sliding mode observer in a mathematical model of the attitude control system of the rotor type aircraft;
the fault-tolerant controller establishing module is used for establishing a fault-tolerant controller according to a mathematical model of a rotary wing type aircraft attitude control system, a fixed time sliding mode observer and a tracking error;
and the control module is used for controlling the rotor type aircraft by utilizing the fault-tolerant controller.
Preferably, the model building module is further configured to: using formulas
Establishing a mathematical model of a rotary wing type aircraft attitude control system;
wherein X represents a state variable of the control system; x1A first state quantity representing the control system,representing the first derivative of a first state quantity of the control system, and X1=[α,β]T,[]TRepresenting the transpose of the matrix, α representing the yaw angle of the rotorcraft, β representing the pitch angle of the rotorcraft;
X2a second state quantity representing the control system,represents a first derivative of a second state quantity of the control system, andrepresenting the first derivative of the yaw angle of the rotorcraft,representing the first derivative of the pitch angle of the rotorcraft;
f (X) represents a nonlinear dynamic value of the control system, andm represents the effective mass of the aircraft, g represents the gravitational acceleration constant, LaRepresenting the length of the center of the rotorcraft from the front motor;
g (X) represents a control input dynamic value of the control system, andKfrepresenting the coefficient of moment generated by the motor, JαRepresenting the moment of inertia of the yaw axis, JβRepresenting the moment of inertia of the pitch axis, LhRepresenting the length of the rear motor from the center of the rotorcraft;
u represents a matrix form of an input amount of the control system, andVfrepresenting the voltage, V, applied to the front motor of a rotorcraftbIndicating application to rotor planeVoltage on the rear motor of the traveling device;
Δ represents the integrated disturbance amount, and Δ ═ Δ f (x)) - (g (x)) + Δ g (x)) (ρi-1)U(t)+D(t),ρiFailure coefficient and p representing occurrence of actuator failure of rotorcrafti∈ [0, 1), Δ f (x) represents a first uncertainty value of the control system, Δ g (x) represents a second uncertainty value of the control system, u (t) represents a functional form of an input quantity of the control system, d (t) represents a functional form of an unknown external disturbance.
Preferably, the sliding-mode observer establishing module is further configured to: using formulas
Establishing a fixed time sliding-mode observer, wherein ξ1Representing a first quantity of state of the sliding-mode observer, ξ2A first state quantity of the sliding mode observer is represented,a first derivative of a first state quantity representing a sliding mode observer,denotes the first derivative of the first state quantity of the sliding-mode observer, e denotes the error variable, and e is ξ1-X1And | | represents a euclidean distance symbol, p is a power exponent and p > 1, sign () represents a sign function, λ1Is a predetermined observer first parameter, λ2Is a predetermined observer second parameter, λ3Is a predetermined observer third parameter, and λ1,λ2,λ3Are all constants greater than zero.
Preferably, the fault-tolerant controller establishing module is further configured to: using formulas
A fault tolerant controller is established in which, among other things,is an estimate of the overall disturbance variable Delta, X2dTo control a second quantity of state X of the system2The reference signal of (a) is set,for virtual controlAn output through a differentiator;is a preset second controller parameter, E2Is a second tracking error, and E2=X2-X2d,Is a first auxiliary function.
Preferably, the fault-tolerant controller establishing module is further configured to: using formulas
Obtaining a functionWherein z represents the state variable of the differentiator, z1Expressed as a first state quantity of the differentiator, z2A second state quantity representing a differentiator, s (t) being the input of the differentiator, representing a perturbation parameter and > 0;
using formulas
using formulas
zout=z2
A differentiator is established in which, among other things,is the first derivative of the first state quantity of the differentiator,is the first derivative of the second state quantity of the differentiator, zoutIs the output of the differentiator.
Preferably, the fault-tolerant controller establishing module is further configured to: using formulas
Constructing a first auxiliary function, wherein E1Is a first tracking error, and E1=x1-x1d,X1dFor controlling a first quantity of state X of the system1The reference signal of (a);
is an auxiliary variable, andis the first derivative of the auxiliary variable,is a preset first controller parameter, wherein,for controlling a first quantity of state X of the system1The first derivative of the reference signal.
The invention has the advantages that: the method comprises the steps of obtaining system fault information by using a fixed time sliding mode observer, designing a mathematical model of a rotor type aircraft attitude control system containing comprehensive disturbance, establishing a fault-tolerant controller according to the mathematical model of the rotor type aircraft attitude control system and the fixed time sliding mode observer, controlling the rotor type aircraft by using the fault-tolerant controller, considering factors such as tracking error and the comprehensive disturbance by using the fault-tolerant controller, having small modeling error, effectively processing faults of the system and having stability.
Drawings
FIG. 1 is a flowchart of an active fault-tolerant control method based on a fixed time observer according to an embodiment of the present invention;
FIG. 2 is a block diagram of an active fault-tolerant control method based on a fixed time observer according to an embodiment of the present invention;
FIG. 3 is a simplified diagram of an experimental platform of an active fault-tolerant control method based on a fixed time observer according to an embodiment of the present invention;
FIG. 4 is a graph illustrating a change in attitude of a rotor-based aircraft at a yaw angle in an active fault-tolerant control method based on a fixed time observer according to an embodiment of the present invention;
fig. 5 is a graph of a change of attitude of a rotor type aircraft in a pitch angle in an active fault-tolerant control method based on a fixed time observer according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the embodiments of the present invention, and it is obvious that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1 and fig. 2, an active fault-tolerant control method based on a fixed time observer is applied to fault-tolerant control of a rotor aircraft, and the method includes:
step S1: establishing a mathematical model of the attitude control system of the rotary wing aircraft according to the control input dynamic value of the control system, the nonlinear dynamic value of the control system and the sum of the comprehensive disturbance quantities; as can be seen in fig. 3, the twin-rotor aircraft has two motors, so that only two attitude angles can be arbitrarily controlled. Therefore, only two attitude angles, namely a yaw angle and a pitch angle, are studied, and the roll angle is changed under the influence of the yaw angle, the pitch angle and the voltage, and the roll angle is not studied.
First, using a formula
Establishing a mathematical model of the attitude control system of the rotor type aircraft without considering the fault condition;
wherein X represents a state variable of the control system; x1A first state quantity representing the control system,representing the first derivative of a first state quantity of the control system, and X1=[α,β]T,[]TRepresenting the transpose of the matrix, α representing the yaw angle of the rotorcraft, β representing the pitch angle of the rotorcraft;
X2a second state quantity representing the control system,represents a first derivative of a second state quantity of the control system, andrepresenting rotorcraftThe first derivative of the yaw angle is,representing the first derivative of the pitch angle of the rotorcraft;
f (X) represents a nonlinear dynamic value of the control system, andm represents the effective mass of the aircraft, g represents the gravitational acceleration constant, LaRepresenting the length of the center of the rotorcraft from the front motor;
g (X) represents a control input dynamic value of the control system, andKfrepresenting the coefficient of moment generated by the motor, JαRepresenting the moment of inertia of the yaw axis, JβRepresenting the moment of inertia of the pitch axis, LhRepresenting the length of the rear motor from the center of the rotorcraft;
u represents a matrix form of an input amount of the control system, andVfrepresenting the voltage, V, applied to the front motor of a rotorcraftbRepresenting the voltage applied to the rear motor of the rotorcraft; the voltage applied to the front motor of the rotorcraft and the voltage applied to the rear motor of the rotorcraft are the objects of the present invention to be obtained, and the present invention finally obtains the voltage applied to the front motor of the rotorcraft and the voltage applied to the rear motor of the rotorcraft, that is, the input of the control system, by derivation in steps and steps.
When considering system uncertainties, external disturbances and actuator faults, the mathematical model of the attitude control system of a rotary wing aircraft can be rewritten in the form:
where Δ F (X) represents a first uncertainty value of the control system, Δ G (X) represents a second uncertainty value of the control system, D represents a matrix form of the unknown external disturbance, UfAn input quantity matrix representing the control system with the fault is represented as follows:
ρi(t) a functional form of a failure coefficient representing the occurrence of an actuator failure of a rotorcraft, tfIndicating the time at which a failure of an actuator of a rotorcraft has occurred, i being an index value of the input variable with the failed control system, uiThe ith failure information element of the control system with the failure.
Then, the integrated disturbance amount is designed, and is represented by a symbol Δ ═ Δ f (x) - (g (x)) + Δ g (x)) (ρ ═ f (x)), (g (x)) + Δ g (x)), (ρi-1)U(t)+D(t),ρiFailure coefficient and p representing occurrence of actuator failure of rotorcrafti∈ [0, 1), U (t) is a function of the input of the control system, D (t) is a function of the unknown external disturbance, and p is expanded by the comprehensive disturbance quantity formulaiU (t) andthe expressions are synonymous and are each an element of the input variable matrix of the control system with the fault.
Finally, the mathematical model of the attitude control system of the rotor aircraft considering the fault condition is as follows:
step S2: establishing a fixed-time sliding mode observer according to an input dynamic value, a nonlinear dynamic value, an error variable and a state quantity of the sliding mode observer in a mathematical model of the attitude control system of the rotor type aircraft; the specific process is as follows: using formulas
Establishing a fixed time sliding-mode observer, wherein ξ1Representing a first quantity of state of the sliding-mode observer, ξ2A first state quantity of the sliding mode observer is represented,a first derivative of a first state quantity representing a sliding mode observer,denotes the first derivative of the first state quantity of the sliding-mode observer, e denotes the error variable, and e is ξ1-X1And | | represents a euclidean distance symbol, p is a power exponent and p > 1, sign () represents a sign function, λ1Is a predetermined observer first parameter, λ2Is a predetermined observer second parameter, λ3Is a predetermined observer third parameter, and λ1,λ2,λ3Are all constants greater than zero.
λ1,λ2,λ3The following conditions are satisfied:
l represents an upper bound of the integrated disturbance quantity.
Step S3: establishing a fault-tolerant controller according to a mathematical model of a rotary wing type aircraft attitude control system, a fixed time sliding mode observer and a tracking error; the specific process is as follows: using formulas
B1=x1-x1d
E2=x2-x2d
Tracking error of design system, wherein E1For the first tracking error, X1dFor controlling a first quantity of state X of the system1The reference signal of (a); e2For the second tracking error, X2dTo control a second quantity of state X of the system2The reference signal of (a);
then, using the formula
Obtaining a functionWherein z represents the state variable of the differentiator, z1Expressed as a first state quantity of the differentiator, z2A second state quantity representing a differentiator, s (t) being the input of the differentiator, representing a perturbation parameter and > 0;
then, using the formula
then, using the formula
zout=z2
A differentiator is established in which, among other things,is a differentialThe first derivative of the first state quantity of the device,is the first derivative of the second state quantity of the differentiator, zoutIs the output of the differentiator.
To further reduce tracking error, equations are used
is an auxiliary variable, andis the first derivative of the auxiliary variable,is a preset first controller parameter, wherein,for controlling a first quantity of state X of the system1The first derivative of the reference signal.
for tracking error E2The derivation was performed with the following results:
through derivation of all the formulas, the formula is constructedAs a fault-tolerant controller, where U is the same as before, is in the form of a matrix of input quantities to the control system, andwill be provided withSubstitution intoThe left side of the front panel is just the right side,is an estimate of the integrated disturbance variable delta,for virtual controlAn output through a differentiator;is a preset second controller parameter,is a first auxiliary function.
Step S4: and controlling the rotor type aircraft by using a fault-tolerant controller.
The stability, i.e. the feasibility, of the scheme designed by the present invention is demonstrated below by the lyapunov function, which is designed as follows:
Wherein the content of the first and second substances,therefore, according to the lyapunov stability theorem, the stability of the system can be ensured when the system fails, and the tracking error can be finally converged into a small range.
FIG. 4 is a graph of attitude change of a rotary wing aircraft at yaw angle, C1The curve represents the reference signal for yaw angle; b is1The curve represents the actual change of the yaw angle under the condition that the controller contains fault compensation; a. the1The curve represents the change in yaw angle without fault compensation by the controller. Comparison B1And A1It can be seen from the curve that when faults and disturbances occur in the system with t > 20s, the fault compensation scheme provided by the invention can enable the yaw angle of the aircraft to be hardly influenced by the faults and disturbances, and when fault compensation is not performed, the yaw angle of the aircraft has large fluctuation.
FIG. 5 is a graph of the attitude change in pitch angle of the rotary wing aircraft of the present invention; b is2The curve represents a reference signal for the pitch angle; a. the2The curve represents the actual change of the depression elevation angle under the condition that the controller contains fault compensation; c2The curve represents the change in pitch angle without fault compensation for the controller. Comparison A2And C2The curve shows that when the system has faults and disturbance in t > 20s, the fault compensation scheme provided by the invention can ensure that the pitching angle of the aircraft is hardly influenced by the faults and disturbance, and when the system has no fault compensation, the pitching angle of the aircraft has large fluctuation.
According to the active fault-tolerant control method based on the fixed time observer, the fault-tolerant controller is established according to the mathematical model of the attitude control system of the rotor type aircraft and the fixed time sliding-mode observer, the rotor type aircraft is controlled by the fault-tolerant controller, faults occurring in the system are effectively processed, the modeling error is small, and the system has stability.
Example 2
Corresponding to embodiment 1 of the present invention, embodiment 2 of the present invention further provides an active fault-tolerant control device based on a fixed time observer, which is applied to fault-tolerant control of a rotor aircraft, and the device includes:
the model establishing module is used for establishing a mathematical model of the attitude control system of the rotary wing aircraft according to the control input dynamic value of the control system, the nonlinear dynamic value of the control system and the sum of the comprehensive disturbance quantities;
the sliding mode observer establishing module is used for establishing a fixed time sliding mode observer according to an input dynamic value, a nonlinear dynamic value, an error variable and a state quantity of the sliding mode observer in a mathematical model of the attitude control system of the rotor type aircraft;
the fault-tolerant controller establishing module is used for establishing a fault-tolerant controller according to a mathematical model of a rotary wing type aircraft attitude control system, a fixed time sliding mode observer and a tracking error;
and the control module is used for controlling the rotor type aircraft by utilizing the fault-tolerant controller.
Preferably, the model building module is further configured to: using formulas
Establishing a mathematical model of a rotary wing type aircraft attitude control system;
wherein X represents a state variable of the control system; x1A first state quantity representing the control system,representing the first derivative of a first state quantity of the control system, and X1=[α,β]T,[]TRepresenting the transpose of the matrix, α representing the yaw angle of the rotorcraft, β representing the pitch angle of the rotorcraft;
X2a second state quantity representing the control system,represents a first derivative of a second state quantity of the control system, andrepresenting the first derivative of the yaw angle of the rotorcraft,representing the first derivative of the pitch angle of the rotorcraft;
f (X) represents a nonlinear dynamic value of the control system, andm represents the effective mass of the aircraft, g represents the gravitational acceleration constant, LaRepresenting the length of the center of the rotorcraft from the front motor;
g (X) represents a control input dynamic value of the control system, andKfrepresenting the coefficient of moment generated by the motor, JαRepresenting the moment of inertia of the yaw axis, JβRepresenting the moment of inertia of the pitch axis, LhRepresenting the length of the rear motor from the center of the rotorcraft;
u represents a matrix form of an input amount of the control system, andVfrepresenting the voltage, V, applied to the front motor of a rotorcraftbIndicating application to the rear of a rotorcraftThe voltage across the motor;
Δ represents the integrated disturbance amount, and Δ ═ Δ f (x)) - (g (x)) + Δ g (x)) (p)i-1)U(t)+D(t),ρiFailure coefficient and p representing occurrence of actuator failure of rotorcrafti∈ [0, 1), Δ f (x) represents a first uncertainty value of the control system, Δ g (x) represents a second uncertainty value of the control system, u (t) represents a functional form of an input quantity of the control system, d (t) represents a functional form of an unknown external disturbance.
Specifically, the sliding-mode observer establishing module is further configured to: using formulas
Establishing a fixed time sliding-mode observer, wherein ξ1Representing a first quantity of state of the sliding-mode observer, ξ2A first state quantity of the sliding mode observer is represented,a first derivative of a first state quantity representing a sliding mode observer,denotes the first derivative of the first state quantity of the sliding-mode observer, e denotes the error variable, and e is ξ1-X1And | | represents a euclidean distance symbol, p is a power exponent and p > 1, sign () represents a sign function, λ1Is a predetermined observer first parameter, λ2Is a predetermined observer second parameter, λ3Is a predetermined observer third parameter, and λ1,λ2,λ3Are all constants greater than zero.
Specifically, the fault-tolerant controller establishing module is further configured to: using formulas
A fault tolerant controller is established in which, among other things,is an estimate of the overall disturbance variable Delta, X2dTo control a second quantity of state X of the system2The reference signal of (a) is set,for virtual controlAn output through a differentiator;is a preset second controller parameter, E2Is a second tracking error, and E2=X2-X2d,Is a first auxiliary function.
Specifically, the fault-tolerant controller establishing module is further configured to: using formulas
Obtaining a functionWherein z represents the state variable of the differentiator, z1Expressed as a first state quantity of the differentiator, z2A second state quantity representing a differentiator, s (t) being the input of the differentiator, representing a perturbation parameter and > 0;
using formulas
using formulas
zout=z2
A differentiator is established in which, among other things,is the first derivative of the first state quantity of the differentiator,is the first derivative of the second state quantity of the differentiator, zoutIs the output of the differentiator.
Specifically, the fault-tolerant controller establishing module is further configured to: using formulas
Constructing a first auxiliary function, wherein E1Is a first tracking error, and E1=x1-x1d,X1dFor controlling a first quantity of state X of the system1The reference signal of (a);
is an auxiliary variable, andis the first derivative of the auxiliary variable,is a preset first controller parameter, wherein,for controlling a first quantity of state X of the system1Reference signal ofThe first derivative of (a).
The above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. An active fault-tolerant control method based on a fixed time observer is applied to fault-tolerant control of a rotor type aircraft, and comprises the following steps:
the method comprises the following steps: establishing a mathematical model of the attitude control system of the rotary wing aircraft according to the control input dynamic value of the control system, the nonlinear dynamic value of the control system and the sum of the comprehensive disturbance quantities;
step two: establishing a fixed-time sliding mode observer according to an input dynamic value, a nonlinear dynamic value, an error variable and a state quantity of the sliding mode observer in a mathematical model of the attitude control system of the rotor type aircraft;
step three: establishing a fault-tolerant controller according to a mathematical model of a rotary wing type aircraft attitude control system, a fixed time sliding mode observer and a tracking error;
step four: and controlling the rotor type aircraft by using a fault-tolerant controller.
2. The active fault-tolerant control method based on the fixed time observer is characterized in that the first step comprises the following steps: using formulas
Establishing a mathematical model of a rotary wing type aircraft attitude control system;
wherein X represents a state variable of the control system; x1A first state quantity representing the control system,representing the first derivative of a first state quantity of the control system, and X1=[α,β]T,[]TRepresenting the transpose of the matrix, α representing the yaw angle of the rotorcraft, β representing the pitch angle of the rotorcraft;
X2a second state quantity representing the control system,represents a first derivative of a second state quantity of the control system, and representing the first derivative of the yaw angle of the rotorcraft,representing the first derivative of the pitch angle of the rotorcraft;
f (X) represents a nonlinear dynamic value of the control system, andm represents the effective mass of the aircraft, g represents the gravitational acceleration constant, LaRepresenting the length of the center of the rotorcraft from the front motor;
g (X) represents a control input dynamic value of the control system, andKfrepresenting the coefficient of moment generated by the motor, JαRepresenting the moment of inertia of the yaw axis, JβRepresenting the moment of inertia of the pitch axis, LhRepresenting the length of the rear motor from the center of the rotorcraft;
u represents a matrix form of an input amount of the control system, andVfrepresenting the voltage, V, applied to the front motor of a rotorcraftbRepresenting the voltage applied to the rear motor of the rotorcraft;
Δ represents the integrated disturbance amount, and Δ ═ Δ f (x)) - (g (x)) + Δ g (x)) (ρi-1)U(t)+D(t),ρiFailure coefficient and p representing occurrence of actuator failure of rotorcrafti∈ [0, 1), Δ f (x) represents a first uncertainty value of the control system, Δ g (x) represents a second uncertainty value of the control system, u (t) represents a functional form of an input quantity of the control system, d (t) represents a functional form of an unknown external disturbance.
3. The active fault-tolerant control method based on the fixed time observer according to claim 2, wherein the second step comprises: using formulas
Establishing a fixed time sliding-mode observer, wherein ξ1Representing a first quantity of state of the sliding-mode observer, ξ2A first state quantity of the sliding mode observer is represented,a first derivative of a first state quantity representing a sliding mode observer,denotes the first derivative of the first state quantity of the sliding-mode observer, e denotes the error variable, and e is ξ1-X1And | | represents a euclidean distance symbol, p is a power exponent and p > 1, sign () represents a sign function, λ1Is a predetermined observer first parameter, λ2Is a predetermined observer second parameter, λ3Is a predetermined observer third parameter, and λ1,λ2,λ3Are all constants greater than zero.
4. The active fault-tolerant control method based on the fixed time observer according to claim 3, wherein the third step comprises: using formulas
A fault tolerant controller is established in which, among other things,is an estimate of the overall disturbance variable Delta, X2dTo control a second quantity of state X of the system2The reference signal of (a) is set,for virtual controlAn output through a differentiator;is a preset second controller parameter, E2Is a second tracking error, and E2=X2-X2d,Is a first auxiliary function.
5. The active fault-tolerant control method based on the fixed time observer according to claim 4, wherein the third step further comprises: using formulas
Obtaining a functionWherein z represents the state variable of the differentiator, z1Expressed as a first state quantity of the differentiator, z2A second state quantity representing a differentiator, s (t) being the input of the differentiator, representing a perturbation parameter and > 0;
using formulas
using formulas
zout=z2
6. The active fault-tolerant control method based on the fixed time observer according to claim 5, wherein the third step further comprises: using formulas
Constructing a first auxiliary function, wherein E1Is a first tracking error, and E1=X1-X1d,X1dFor controlling a first quantity of state X of the system1The reference signal of (a);
7. An active fault-tolerant control device based on a fixed time observer, which is applied to fault-tolerant control of a rotor type aircraft, and comprises:
the model establishing module is used for establishing a mathematical model of the attitude control system of the rotary wing aircraft according to the control input dynamic value of the control system, the nonlinear dynamic value of the control system and the sum of the comprehensive disturbance quantities;
the sliding mode observer establishing module is used for establishing a fixed time sliding mode observer according to an input dynamic value, a nonlinear dynamic value, an error variable and a state quantity of the sliding mode observer in a mathematical model of the attitude control system of the rotor type aircraft;
the fault-tolerant controller establishing module is used for establishing a fault-tolerant controller according to a mathematical model of a rotary wing type aircraft attitude control system, a fixed time sliding mode observer and a tracking error;
and the control module is used for controlling the rotor type aircraft by utilizing the fault-tolerant controller.
8. The active fault-tolerant control device based on the fixed-time observer of claim 7, wherein the model building module is further configured to: using formulas
Establishing a mathematical model of a rotary wing type aircraft attitude control system;
wherein X represents a state variable of the control system; x1A first state quantity representing the control system,representing the first derivative of a first state quantity of the control system, and X1=[α,β]T,[]TRepresenting the transpose of the matrix, α representing the yaw angle of the rotorcraft, β representing the pitch angle of the rotorcraft;
X2a second state quantity representing the control system,represents a first derivative of a second state quantity of the control system, and representing the first derivative of the yaw angle of the rotorcraft,representing the first derivative of the pitch angle of the rotorcraft;
f (X) represents a nonlinear dynamic value of the control system, andm represents the effective mass of the aircraft, g represents the gravitational acceleration constant, LaRepresenting the length of the center of the rotorcraft from the front motor;
g (X) represents a control input dynamic value of the control system, andKfrepresenting the coefficient of moment generated by the motor, JαRepresenting the moment of inertia of the yaw axis, JβRepresenting the moment of inertia of the pitch axis, LhRepresenting the length of the rear motor from the center of the rotorcraft;
u represents a matrix form of an input amount of the control system, andVfrepresenting the voltage, V, applied to the front motor of a rotorcraftbRepresenting the voltage applied to the rear motor of the rotorcraft;
Δ represents the integrated disturbance amount, and Δ ═ Δ f (x)) - (g (x)) + Δ g (x)) (ρi-1)U(t)+D(t),ρiFailure coefficient and p representing occurrence of actuator failure of rotorcrafti∈ [0, 1), Δ F (X) represents a first uncertainty value of the control system, Δ G (X) represents a second uncertainty value of the control system, U (t) represents the control systemD (t) represents a functional form of the unknown external disturbance.
9. The active fault-tolerant control device based on the fixed-time observer according to claim 8, wherein the sliding-mode observer establishing module is further configured to: using formulas
Establishing a fixed time sliding-mode observer, wherein ξ1Representing a first quantity of state of the sliding-mode observer, ξ2A first state quantity of the sliding mode observer is represented,a first derivative of a first state quantity representing a sliding mode observer,denotes the first derivative of the first state quantity of the sliding-mode observer, e denotes the error variable, and e is ξ1-X1And | | represents a euclidean distance symbol, p is a power exponent and p > 1, sign () represents a sign function, λ1Is a predetermined observer first parameter, λ2Is a predetermined observer second parameter, A3Is a predetermined observer third parameter, and λ1,λ2,λ3Are all constants greater than zero.
10. The active fault-tolerant control device based on the fixed-time observer of claim 9, wherein the fault-tolerant controller establishing module is further configured to: using formulas
A fault tolerant controller is established in which, among other things,is an estimate of the overall disturbance variable Delta, X2dTo control a second quantity of state X of the system2The reference signal of (a) is set,for virtual controlAn output through a differentiator;is a preset second controller parameter, E2Is a second tracking error, and E2=X2-X2d,Is a first auxiliary function.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010381563.5A CN111522241B (en) | 2020-05-08 | 2020-05-08 | Active fault-tolerant control method and device based on fixed time observer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010381563.5A CN111522241B (en) | 2020-05-08 | 2020-05-08 | Active fault-tolerant control method and device based on fixed time observer |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111522241A true CN111522241A (en) | 2020-08-11 |
CN111522241B CN111522241B (en) | 2020-12-29 |
Family
ID=71905177
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010381563.5A Active CN111522241B (en) | 2020-05-08 | 2020-05-08 | Active fault-tolerant control method and device based on fixed time observer |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111522241B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112650303A (en) * | 2021-01-20 | 2021-04-13 | 南京航空航天大学 | Unmanned aerial vehicle formation fault-tolerant control method based on fixed time technology |
CN113064349A (en) * | 2021-03-22 | 2021-07-02 | 中国人民解放军国防科技大学 | Fixed time convergence nonlinear platform sliding mode control method, device and system |
CN113204193A (en) * | 2021-05-06 | 2021-08-03 | 北京航空航天大学 | Fault-tolerant control method and device for aircraft and electronic equipment |
CN113325717A (en) * | 2021-06-10 | 2021-08-31 | 哈尔滨工业大学 | Optimal fault-tolerant control method, system, processing equipment and storage medium based on interconnected large-scale system |
CN113342033A (en) * | 2021-05-26 | 2021-09-03 | 南京航空航天大学 | Fixed-time-technology-based attitude synchronization fault-tolerant control method for fixed-wing unmanned aerial vehicle |
CN113359447A (en) * | 2021-06-03 | 2021-09-07 | 南京航空航天大学 | Compound rapid nonsingular terminal sliding mode control method for atmospheric entrance section of Mars detector |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030229407A1 (en) * | 2002-06-06 | 2003-12-11 | Honda Giken Kogyo Kabushiki Kaisha | Control system for plant |
EP0915399B1 (en) * | 1997-09-16 | 2005-01-19 | Honda Giken Kogyo Kabushiki Kaisha | Plant control system |
CA2690318A1 (en) * | 2009-01-16 | 2010-07-16 | Messier-Bugatti | Process for managing the attitude control of aircraft landing gear |
CN103135553A (en) * | 2013-01-21 | 2013-06-05 | 南京航空航天大学 | Four-rotor aircraft fault-tolerant control method |
CN103838145A (en) * | 2014-01-17 | 2014-06-04 | 南京航空航天大学 | Vertical take-off and landing airplane robust fault-tolerant control system and method based on cascaded observers |
CN104102225A (en) * | 2014-06-27 | 2014-10-15 | 金陵科技学院 | Unmanned aircraft posture decentralized fault-tolerant control method based on actuator dynamic state |
CN104965414A (en) * | 2015-06-30 | 2015-10-07 | 天津大学 | Tolerant control method for partial failure of four-rotor unmanned aerial vehicle actuator |
CN105353615A (en) * | 2015-11-10 | 2016-02-24 | 南京航空航天大学 | Active fault tolerance control method of four-rotor aircraft based on sliding-mode observer |
CN105843240A (en) * | 2016-04-08 | 2016-08-10 | 北京航空航天大学 | Spacecraft attitude integral sliding mode fault tolerance control method taking consideration of performer fault |
CN106647584A (en) * | 2017-01-17 | 2017-05-10 | 南京航空航天大学 | Fault tolerant control method of four-rotor-wing aircraft based on optimal sliding mode |
CN108628167A (en) * | 2018-05-08 | 2018-10-09 | 南京邮电大学 | The fault tolerant control method of flight control system based on self-adjusted block technology |
CN109144084A (en) * | 2018-07-11 | 2019-01-04 | 哈尔滨工业大学 | A kind of VTOL Reusable Launch Vehicles Attitude tracking control method based on set time Convergence monitoring device |
CN109188910A (en) * | 2018-09-28 | 2019-01-11 | 浙江工业大学 | A kind of fault-tolerant tracking and controlling method of the adaptive neural network of rigid aircraft |
CN109343369A (en) * | 2018-11-19 | 2019-02-15 | 南京邮电大学 | A kind of quadrotor fault controller method based on nonlinear observer |
CN110333656A (en) * | 2019-07-04 | 2019-10-15 | 南京航空航天大学 | A kind of flexible spacecraft fault tolerant control method based on interacted system method |
-
2020
- 2020-05-08 CN CN202010381563.5A patent/CN111522241B/en active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0915399B1 (en) * | 1997-09-16 | 2005-01-19 | Honda Giken Kogyo Kabushiki Kaisha | Plant control system |
US20030229407A1 (en) * | 2002-06-06 | 2003-12-11 | Honda Giken Kogyo Kabushiki Kaisha | Control system for plant |
CA2690318A1 (en) * | 2009-01-16 | 2010-07-16 | Messier-Bugatti | Process for managing the attitude control of aircraft landing gear |
CN103135553A (en) * | 2013-01-21 | 2013-06-05 | 南京航空航天大学 | Four-rotor aircraft fault-tolerant control method |
CN103838145A (en) * | 2014-01-17 | 2014-06-04 | 南京航空航天大学 | Vertical take-off and landing airplane robust fault-tolerant control system and method based on cascaded observers |
CN104102225A (en) * | 2014-06-27 | 2014-10-15 | 金陵科技学院 | Unmanned aircraft posture decentralized fault-tolerant control method based on actuator dynamic state |
CN104965414A (en) * | 2015-06-30 | 2015-10-07 | 天津大学 | Tolerant control method for partial failure of four-rotor unmanned aerial vehicle actuator |
CN105353615A (en) * | 2015-11-10 | 2016-02-24 | 南京航空航天大学 | Active fault tolerance control method of four-rotor aircraft based on sliding-mode observer |
CN105843240A (en) * | 2016-04-08 | 2016-08-10 | 北京航空航天大学 | Spacecraft attitude integral sliding mode fault tolerance control method taking consideration of performer fault |
CN106647584A (en) * | 2017-01-17 | 2017-05-10 | 南京航空航天大学 | Fault tolerant control method of four-rotor-wing aircraft based on optimal sliding mode |
CN108628167A (en) * | 2018-05-08 | 2018-10-09 | 南京邮电大学 | The fault tolerant control method of flight control system based on self-adjusted block technology |
CN109144084A (en) * | 2018-07-11 | 2019-01-04 | 哈尔滨工业大学 | A kind of VTOL Reusable Launch Vehicles Attitude tracking control method based on set time Convergence monitoring device |
CN109188910A (en) * | 2018-09-28 | 2019-01-11 | 浙江工业大学 | A kind of fault-tolerant tracking and controlling method of the adaptive neural network of rigid aircraft |
CN109343369A (en) * | 2018-11-19 | 2019-02-15 | 南京邮电大学 | A kind of quadrotor fault controller method based on nonlinear observer |
CN110333656A (en) * | 2019-07-04 | 2019-10-15 | 南京航空航天大学 | A kind of flexible spacecraft fault tolerant control method based on interacted system method |
Non-Patent Citations (5)
Title |
---|
PENG CHENG: "Active fault tolerant control design for UAV using nonsingular fast terminal sliding mode approach", 《2018 CHINESE CONTROL AND DECISION CONFERENCE (CCDC)》 * |
XIAOJUN XING: "Fault-tolerant flight control of quad-rotor UAV based on sliding mode theory", 《 2018 CHINESE CONTROL AND DECISION CONFERENCE (CCDC)》 * |
翟丽相: "无人机姿态控制系统滑模动态面容错控制方法", 《系统仿真技术》 * |
郑佳静: "采用滑模观测器的四旋翼无人机执行器加性故障容错控制", 《华侨大学学报(自然科学版)》 * |
韩嘉俊: "故障自适应滑模再入姿态控制设计", 《现代防御技术》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112650303A (en) * | 2021-01-20 | 2021-04-13 | 南京航空航天大学 | Unmanned aerial vehicle formation fault-tolerant control method based on fixed time technology |
CN113064349A (en) * | 2021-03-22 | 2021-07-02 | 中国人民解放军国防科技大学 | Fixed time convergence nonlinear platform sliding mode control method, device and system |
CN113204193A (en) * | 2021-05-06 | 2021-08-03 | 北京航空航天大学 | Fault-tolerant control method and device for aircraft and electronic equipment |
CN113342033A (en) * | 2021-05-26 | 2021-09-03 | 南京航空航天大学 | Fixed-time-technology-based attitude synchronization fault-tolerant control method for fixed-wing unmanned aerial vehicle |
CN113342033B (en) * | 2021-05-26 | 2022-10-04 | 南京航空航天大学 | Fixed time technology-based attitude synchronization fault-tolerant control method for fixed-wing unmanned aerial vehicle |
CN113359447A (en) * | 2021-06-03 | 2021-09-07 | 南京航空航天大学 | Compound rapid nonsingular terminal sliding mode control method for atmospheric entrance section of Mars detector |
CN113325717A (en) * | 2021-06-10 | 2021-08-31 | 哈尔滨工业大学 | Optimal fault-tolerant control method, system, processing equipment and storage medium based on interconnected large-scale system |
CN113325717B (en) * | 2021-06-10 | 2022-01-28 | 哈尔滨工业大学 | Optimal fault-tolerant control method, system, processing equipment and storage medium based on interconnected large-scale system |
Also Published As
Publication number | Publication date |
---|---|
CN111522241B (en) | 2020-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111522241B (en) | Active fault-tolerant control method and device based on fixed time observer | |
Luque-Vega et al. | Robust block second order sliding mode control for a quadrotor | |
Falcón et al. | Comparative analysis of continuous sliding-modes control strategies for quad-rotor robust tracking | |
Chen et al. | Design of flight control system for a novel tilt-rotor UAV | |
Lin et al. | Adaptive flight control design for quadrotor UAV based on dynamic inversion and neural networks | |
Bousbaine et al. | Modelling and simulation of a quad-rotor helicopter | |
Ali et al. | Attitude and altitude control of trirotor UAV by using adaptive hybrid controller | |
CN113296525A (en) | Data-driven nonlinear control method for tilting three-rotor unmanned aerial vehicle | |
Öner et al. | LQR and SMC stabilization of a new unmanned aerial vehicle | |
Ghasemi et al. | Control of quadrotor using sliding mode disturbance observer and nonlinear H∞ | |
CN112327896A (en) | Rotor fault-tolerant control method and device, computer storage medium and computer equipment | |
Wu et al. | Robust adaptive finite‐time trajectory tracking control of a quadrotor aircraft | |
Su et al. | Modeling and control of a class of urban air mobility tiltrotor aircraft | |
Hegde et al. | Transition flight modeling and robust control of a VTOL unmanned quad tilt-rotor aerial vehicle | |
Song et al. | Research on attitude control of quadrotor uav based on active disturbance rejection control | |
Idrissi et al. | Mathematical & physical modelling of a quadrotor UAV | |
Mokhtari et al. | Robust control for attitude tracking problem for a quadrotor unmanned aerial vehicle | |
Tang et al. | Second-order sliding mode attitude controller design of a small-scale helicopter | |
Staub et al. | Modeling and design of a lift-augmented quadcopter | |
Xu et al. | Modelling and hovering control of a novel multi-tandem ducted fan vehicle | |
Liu et al. | Identification of attitude flight dynamics for an unconventional UAV | |
Loubar et al. | Altitude back-stepping control of quadcopter | |
Hu et al. | Modeling and hovering control of 5-DoF tilt-birotor robot | |
Bilici | Modeling and Control of a Fixed-Wing High-Speed UAV | |
Fan et al. | System identification and attitude control of a small scale unmanned helicopter |
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 |