CN110985216A - Intelligent multivariable control method for aero-engine with online correction - Google Patents
Intelligent multivariable control method for aero-engine with online correction Download PDFInfo
- Publication number
- CN110985216A CN110985216A CN201911132301.9A CN201911132301A CN110985216A CN 110985216 A CN110985216 A CN 110985216A CN 201911132301 A CN201911132301 A CN 201911132301A CN 110985216 A CN110985216 A CN 110985216A
- Authority
- CN
- China
- Prior art keywords
- engine
- narma
- intelligent
- model
- multivariable
- 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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
- F02C9/00—Controlling gas-turbine plants; Controlling fuel supply in air- breathing jet-propulsion plants
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
- F05D2270/00—Control
- F05D2270/70—Type of control algorithm
- F05D2270/709—Type of control algorithm with neural networks
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses an intelligent multivariable control method for an aeroengine with online correction, which comprises the following steps: establishing an NARMA-L2 model of an engine rotating speed and pressure ratio control system by adopting a neural network group; and (3) designing a controller online correction module based on an NARMA-L2 linear compensation strategy and a gradient descent algorithm by combining an NARMA-L2 model to obtain the intelligent multivariable controller of the aero-engine with self-adaptive capacity within a certain working condition range in an envelope. The invention solves the problems that the traditional neural network prediction model based establishment and the online solution of the control law are difficult, and the control quality is reduced when the controller is expanded and applied in the envelope, is suitable for the multivariable control of the engine speed and the pressure ratio at different working points in a certain flight envelope, and has positive promotion effects on eliminating the steady-state error of an engine control system, and improving the application range and the control quality in the controller envelope.
Description
Technical Field
The invention belongs to the technical field of multi-variable control of aero-engines, and particularly relates to an intelligent multi-variable control method for an aero-engine with online correction.
Background
With the rapid development of modern aviation technology, aircraft require aircraft engines with stronger performance and higher reliability. An aircraft engine is a complex thermodynamic system with strong uncertainties and time-varying properties that make the design of the controller difficult. Classical control theory and modern control theory rely on accurate mathematical models since the subject of the study is mainly a linear time invariant system, which leads to great difficulties in the application of such strongly non-linear, time variant and complex control systems in aircraft engines. Researchers have therefore attempted to replace complex systems with simple approximation models and judge their value in terms of their ease of analysis, ease of handling, and likelihood of application in practical problems. Although one recognizes from the outset that practical dynamic systems are nonlinear, linearization within a small neighborhood of equilibrium states is mathematically tractable. In particular the superposition principle, makes the calculation in a linear manner sufficiently simple. Most importantly, the linearized model of most nonlinear systems can obtain satisfactory precision in a normal use range, and therefore, the nonlinear system is widely applied.
The intelligent control method is a control theory based on an artificial intelligence algorithm close to a human thinking mode, and the neural network method is an important branch of the intelligent control theory. In recent years, rapid advances in technology have led to systems that are often required to operate in regions of state space that are not sufficiently linearized. To solve this problem, neural networks are widely used due to their good approximation and generalization capabilities. A Nonlinear autoregressive Moving Average with Feedback Linearization (NARMA-L2) controller is an effective artificial neural network controller architecture. Under certain conditions, the input-output relationship of the nonlinear system can be identified by a NARMA-L2 model, and the control law can be obtained by simple mathematical transformation. The model was first introduced by Narendra and Mukhopadhyay and is now widely used in many non-linear systems.
The NARMA-L2 model is a neural network approximation model that has modeling errors and training errors. In order to improve the applicability of the controller in the flight envelope, the invention provides an intelligent multivariable control method for an aircraft engine with online correction. The method utilizes a gradient descent algorithm to carry out online correction on neural network parameters in the NARMA-L2 model, so that the neural network parameters have self-adaptive characteristics, and simultaneously provides a NARMA-L2 linear compensation strategy, thereby further improving the control quality of the controller in an envelope range.
Disclosure of Invention
Aiming at the technical problems, the invention provides an intelligent multivariable control method for an aero-engine with online correction, which solves the problems that the traditional neural network prediction model-based establishment and online solution of a control law are difficult, and the control quality is reduced due to the fact that a controller is expanded and applied in an envelope.
The technical scheme is as follows: in order to achieve the purpose, the invention adopts the technical scheme that:
an intelligent multivariable control method for an aircraft engine with online correction is characterized by comprising the following steps:
step A), constructing a NARMA-L2 model of the intelligent multivariable control system of the aircraft engine by adopting a neural network group according to the relationship between the main fuel quantity and the rotating speed of a high-pressure rotor of an engine model and the area of a throat of a tail nozzle and the pressure ratio of the whole aircraft engine;
and step B) combining with an NARMA-L2 model, designing a controller online correction module based on an NARMA-L2 linear compensation strategy and a gradient descent algorithm, and obtaining the intelligent multivariable controller of the aero-engine with self-adaptive capacity within a certain working condition range in an envelope.
Further, in the step A), according to the relation between the main fuel quantity and the high-pressure rotor rotating speed of the engine model and the relation between the area of the throat of the tail nozzle and the overall pressure ratio, the NARMA-L2 model of the intelligent multivariable control system of the aero-engine is constructed by adopting a neural network group, and the specific steps are as follows:
step A1), selecting the main fuel quantity of the engine as an input quantity u1[k]High pressure rotor speed as output y1[k](ii) a Selecting the area of the throat of the engine tail nozzle as an input quantity u2[k]The total pressure ratio is taken as the output y2[k]And constructing a double-loop NARMA-L2 model according to the selected input and output quantities as follows:
y1[k]=f1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])+g1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])u1[k]
y2[k]=f2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])+g2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])u2[k]
wherein f is1(·)、g1(·)、f2(·)、g2Can be connected by four feedforward neural networks f1、g1、f2And g2Obtaining an approximation;
step A2), f is performed1、g1、f2And g2And (4) randomly initializing four neural network structures, and performing off-line training on the neural networks by using a gradient descent method according to the training samples to obtain a NARMA-L2 model of the intelligent multivariable control system of the trained aircraft engine.
Further, the specific steps of designing a controller online correction module based on an NARMA-L2 linear compensation strategy and a gradient descent algorithm by combining the NARMA-L2 model in the step B) to obtain the intelligent multivariable controller of the aero-engine with self-adaptive capability within a certain working condition range in an envelope are as follows:
step B1), constructing a structure diagram of an aircraft engine intelligent multivariable control system of a neural network, wherein the structure diagram comprises an NARMA-L2 model, a rotating speed controller, a pressure ratio controller, an engine model and an online correction module;
and B2), designing an online correction module based on a gradient descent algorithm and a NARMA-L2 linear compensation strategy to obtain the intelligent multivariable controller of the aero-engine.
Further, the step B2) is to design an online correction module based on a gradient descent algorithm and a NARMA-L2 linear compensation strategy, and the specific steps of obtaining the intelligent multivariable controller of the aircraft engine are as follows:
step B2.1), respectively obtaining the initial control laws of the rotating speed and the pressure ratio according to the NARMA-L2 model obtained by the training of the step A2):
wherein, y1 *[k]And y2 *[k]Respectively outputting expected rotation speed of a high-pressure rotor of the engine and the overall pressure ratio; f. of1(k)、g1(k)、f2(k) And g2(k) The expression of (a) is:
f1(k)=f1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])
g1(k)=g1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])
f2(k)=f2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])
g2(k)=g2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])
step B2.2), based on the error between the desired output and the actual output, the u is compensated for by using the NARMA-L2 linear compensation strategy1[k]And u2[k]Compensating and modifying the neural network f according to a gradient descent algorithm1、g1、f2And g2The topology parameter of (2).
Further, in the step B2.2), the pair u is linearly compensated by using a NARMA-L2 linear compensation strategy according to the error between the expected output and the actual output1[k]And u2[k]Compensating and modifying the neural network f according to a gradient descent algorithm1、g1、f2And g2The topology parameters are as follows:
step B2.2.1), judging whether the current time is the initial time, if so, not compensating; if the current time is not the initial time, calculating a compensation value delta u [ k ] of the initial control quantity at the current time according to the actual input u [ k-1] and the actual output y [ k-1] of the engine at the previous time and the neural network outputs f (k-1) and g (k-1) in the NARMA-L2 model, wherein the specific formula is as follows:
α is compensation gain, and is determined by a multi-trial method;
step B2.2.2), obtaining the control quantity u according to the initial control law at the current timeoriginal[k]Compensation value delta u k with control quantity]And calculating to obtain the actual control quantity u [ k ]]The concrete formula is as follows:
u1[k]=u1_original[k]+Δu1[k]
u2[k]=u2_original[k]+Δu2[k]
step B2.2.3) applying the control quantity to the aircraft engine to obtain the actual output y [ k ] at the current moment]Using the error e between the desired output and the actual outputC[k]For the neural network group f1、g1、f2And g2The topological parameters are corrected on line, and the specific formula is as follows:
W(k+1)=W(k)+βΔW(k)
V(k+1)=V(k)+βΔV(k)
wherein, W and V are weight and threshold of the neural network group respectively, β is learning rate, the adjusting speed of the size influence parameter is determined by a multi-trial method, Δ W (k) and Δ V (k) are Jacobian matrixes which respectively represent increment of the weight and the threshold, and the updating of each element in the topological vector at the moment k is as follows:
wi(k+1)=wi(k)-βΔwi(k)
vi(k+1)=vi(k)-βΔvi(k)
wherein, T ═ y (k), y (k-1), y (k-2),. ·, y (k-n), u (k-1), u (k-2), …, u (k-n) }; w is ai、viThe ith elements representing W and V, respectively; g0Denotes g1Or g2;Δwi、ΔviRespectively, the i-th elements of Δ w (k) and Δ v (k).
And obtaining the rotating speed and pressure ratio controller at the k +1 moment according to the updated NARMA-L2 model. To this end, the intelligent multivariable controller design ends.
Has the advantages that: the invention solves the problems that the traditional neural network prediction model based establishment and the online solution of the control law are difficult, and the control quality is reduced when the controller is expanded and applied in the envelope, is suitable for the multivariable control of the engine speed and the pressure ratio at different working points in a certain flight envelope, and has positive promotion effects on eliminating the steady-state error of an engine control system, and improving the application range and the control quality in the controller envelope.
Drawings
FIG. 1 is a diagram of the intelligent multivariable control structure of the aeroengine based on the neural network.
FIG. 2 is a fuel quantity input data sample set distribution plot.
FIG. 3 is a sample set distribution plot of high pressure rotor speed output data.
FIG. 4 is a sample set distribution plot of jet nozzle throat area input data.
Fig. 5 is a plot of the whole press ratio output data sample set.
FIG. 6 is a schematic representation of H, Ma as a function of time.
Fig. 7 is a schematic diagram of the change in fuel amount.
FIG. 8 is a schematic representation of the change in the area of the jet nozzle throat.
FIG. 9 is a schematic of the high pressure rotor speed versus command signal tracking.
FIG. 10 is a schematic diagram of a complete press ratio versus a trace of a command signal.
Detailed Description
The following further describes embodiments of the present invention with reference to the drawings.
The invention discloses an intelligent multivariable control method for an aeroengine with online correction, which is characterized by comprising the following steps of:
step A), constructing a NARMA-L2 model of the intelligent multivariable control system of the aircraft engine by adopting a neural network group according to the relationship between the main fuel quantity and the rotating speed of a high-pressure rotor of an engine model and the area of a throat of a tail nozzle and the pressure ratio of the whole aircraft engine;
and step B) combining with an NARMA-L2 model, designing a controller online correction module based on an NARMA-L2 linear compensation strategy and a gradient descent algorithm, and obtaining the intelligent multivariable controller of the aero-engine with self-adaptive capacity within a certain working condition range in an envelope.
Step A1), selecting the main fuel quantity of the engine as an input quantity u1[k]High pressure rotor speed as output y1[k](ii) a Selecting the area of the throat of the engine tail nozzle as an input quantity u2[k]The total pressure ratio is taken as the output y2[k]And constructing a double-loop NARMA-L2 model according to the selected input and output quantities as follows:
y1[k]=f1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])+g1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])u1[k]
y2[k]=f2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])+g2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])u2[k]
wherein f is1(·)、g1(·)、f2(·)、g2Can be connected by four feedforward neural networks f1、g1、f2And g2Obtaining an approximation;
step A2), f is performed1、g1、f2And g2And (4) randomly initializing four neural network structures, and performing off-line training on the neural networks by using a gradient descent method according to the training samples to obtain a NARMA-L2 model of the intelligent multivariable control system of the trained aircraft engine.
Step B1), constructing a structure diagram of the neural network aeroengine intelligent multivariable control system, wherein the structure diagram comprises an NARMA-L2 model, a rotating speed control formula, a pressure ratio control formula, an engine model and an online correction module;
and B2), designing an online correction module based on a gradient descent algorithm and a NARMA-L2 linear compensation strategy to obtain the intelligent multivariable controller of the aero-engine.
Further, the step B2) is to design an online correction module based on a gradient descent algorithm and a NARMA-L2 linear compensation strategy, and the specific steps of obtaining the intelligent multivariable controller of the aircraft engine are as follows:
step B2.1), respectively obtaining the initial control laws of the rotating speed and the pressure ratio according to the NARMA-L2 model obtained by the training of the step A2):
wherein, y1 *[k]And y2 *[k]Respectively outputting expected rotation speed of a high-pressure rotor of the engine and the overall pressure ratio; f. of1(k)、g1(k)、f2(k) And g2(k) The expression of (a) is:
f1(k)=f1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])
g1(k)=g1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])
f2(k)=f2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])
g2(k)=g2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])
step B2.2), based on the error between the desired output and the actual output, the u is compensated for by using the NARMA-L2 linear compensation strategy1[k]And u2[k]Compensating and modifying the neural network f according to a gradient descent algorithm1、g1、f2And g2The topology parameter of (2).
Step B2.2.1), judging whether the current time is the initial time, if so, not compensating; if the current time is not the initial time, calculating a compensation value delta u [ k ] of the initial control quantity at the current time according to the actual input u [ k-1] and the actual output y [ k-1] of the engine at the previous time and the neural network outputs f (k-1) and g (k-1) in the NARMA-L2 model, wherein the specific formula is as follows:
α is compensation gain, and is determined by a multi-trial method;
step B2.2.2), obtaining the control quantity u according to the initial control law at the current timeoriginal[k]Compensation value delta u k with control quantity]And calculating to obtain the actual control quantity u [ k ]]The concrete formula is as follows:
u1[k]=u1_original[k]+Δu1[k]
u2[k]=u2_original[k]+Δu2[k]
step B2.2.3) applying the control quantity to the aircraft engine to obtain the actual output y [ k ] at the current moment]Using the error e between the desired output and the actual outputC[k]For the neural network group f1、g1、f2And g2The topological parameters are corrected on line, and the specific formula is as follows:
W(k+1)=W(k)+βΔW(k)
V(k+1)=V(k)+βΔV(k)
wherein, W and V are weight and threshold of the neural network group respectively, β is learning rate, the adjusting speed of the size influence parameter is determined by a multi-trial method, and Δ W (k) and Δ V (k) are Jacobian matrixes, the updating of each element in the topological vector at the k moment is as follows:
wi(k+1)=wi(k)-βΔwi(k)
vi(k+1)=vi(k)-βΔvi(k)
wherein, T ═ y (k), y (k-1), y (k-2),.., y (k-n), u (k-1), u (k-2), …, u (k-n) }.
And obtaining the rotating speed and pressure ratio controller at the k +1 moment according to the updated NARMA-L2 model. To this end, the intelligent multivariable controller design ends.
In order to verify the effectiveness of the neural network-based intelligent multivariable control method for the aircraft engine, digital simulation of the takeoff, acceleration and deceleration processes of the engine in a certain envelope is carried out in an MATLAB environment.
The invention adopts a nonlinear part model of a certain type double-rotor turbofan engine with small bypass ratio as a controlled object. The model is constructed by an object-oriented programming method, comprises important parts of an aircraft engine such as an air inlet channel, a fan, an air compressor, a combustion chamber, a turbine, a tail nozzle and the like, and is easy to call in an MATLAB environment.
Before the simulation verification of the performance of the controller is carried out, the two sets of NARMA-L2 models are respectively identified by utilizing a neural network. Each model comprises two neural networks which are single hidden layers, and the number of the neurons of the input layer, the hidden layers and the output layer is 6, 10 and 1 respectively. Before training, an aircraft engine model is used for generating a main fuel quantity-high-pressure rotor rotating speed data set containing 20000 groups of training samples under the flight conditions of H & lt0 & gt and Ma & lt0 & gt, and 20000 groups of tail pipe fuel quantity-total machine pressure ratio data sets are generated by the same method. The ranges of the individual parameters in the data set are shown in fig. 2-5, and the NARMA-L2 model of the neural network set is trained based on the generated training data.
In order to verify the applicability of the method provided by the invention in the flight envelope, the engine running track is selected in the full envelope range, as shown in fig. 6, when the engine is in a state, and the change of the control quantity and the rotating speed-pressure ratio tracking simulation result are shown in fig. 7-10.
The simulation result shows that the engine rotating speed and the pressure ratio command signal change along with the change of the flight condition, and the intelligent multivariable control method provided by the invention has good applicability in the envelope, so that the rotating speed of the high-pressure rotor of the engine and the pressure ratio of the whole engine can track the command signal, and the dynamic response is good. Simulation results show that the intelligent multivariable control method for the aero-engine with online correction has good control quality in a flight envelope.
The invention designs an intelligent multivariable control method for an aircraft engine with online correction, which solves the problems that the traditional neural network prediction model based establishment and online solution of a control law are difficult, the controller is expanded and applied in a covered line, and the control quality is reduced, is suitable for multivariable control of the engine speed and the pressure ratio of different working points in a certain flight covered line, and has positive promotion effects on eliminating steady-state errors of an engine control system, and improving the application range and the control quality in the covered line of the controller.
It should be noted that the above mentioned embodiments are only specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any changes and substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the appended claims.
Claims (5)
1. An intelligent multivariable control method for an aircraft engine with online correction is characterized by comprising the following steps:
step A), constructing a NARMA-L2 model of the intelligent multivariable control system of the aircraft engine by adopting a neural network group according to the relationship between the main fuel quantity and the rotating speed of a high-pressure rotor of an engine model and the area of a throat of a tail nozzle and the pressure ratio of the whole aircraft engine;
and step B) combining with an NARMA-L2 model, designing a controller online correction module based on an NARMA-L2 linear compensation strategy and a gradient descent algorithm, and obtaining the intelligent multivariable controller of the aero-engine with self-adaptive capacity within a certain working condition range in an envelope.
2. The intelligent multivariable control method for the aero-engine with online correction as claimed in claim 1, wherein: in the step A), according to the relation between the main fuel quantity and the rotating speed of a high-pressure rotor of an engine model and the relation between the area of a throat of a tail nozzle and the pressure ratio of the whole aircraft, a neural network group is adopted to construct a NARMA-L2 model of the intelligent multivariable control system of the aircraft engine, and the specific steps are as follows:
step A1), selecting the main fuel quantity of the engine as an input quantity u1[k]High pressure rotor speed as output y1[k]K represents time; selecting the area of the throat of the engine tail nozzle as an input quantity u2[k]The total pressure ratio is taken as the output y2[k]And constructing a dual-loop NARMA-L2 model according to the selected input and output quantities as follows:
y1[k]=f1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])+g1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])u1[k]
y2[k]=f2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])+g2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])u2[k]
wherein f is1(·)、g1(·)、f2(·)、g2Through four feedforward neural networks f1、g1、f2And g2Obtaining an approximation;
step A2), f is performed1、g1、f2And g2Four feedforward neural network structures are initialized randomly, and a gradient descent method is used for a neural network group f according to training samples1、g1、f2And g2And performing off-line training to obtain a NARMA-L2 model of the trained aircraft engine intelligent multivariable control system.
3. The intelligent multivariable control method for the aero-engine with online correction as claimed in claim 2, wherein: the step B) is combined with an NARMA-L2 model, a controller online correction module based on an NARMA-L2 linear compensation strategy and a gradient descent algorithm is designed, and the specific steps of obtaining the intelligent multivariable controller of the aircraft engine with self-adaptive capacity in a certain working condition range in an envelope are as follows:
step B1), constructing a structure diagram of an aircraft engine intelligent multivariable control system of a neural network, wherein the structure diagram comprises an NARMA-L2 model, a rotating speed controller, a pressure ratio controller, an engine model and an online correction module;
and B2), designing an online correction module based on a gradient descent algorithm and a NARMA-L2 linear compensation strategy to obtain the intelligent multivariable controller of the aero-engine.
4. An intelligent multivariable control method for an aircraft engine with online correction as claimed in claim 3, wherein: the step B2) is designed into an online correction module based on a gradient descent algorithm and a NARMA-L2 linear compensation strategy, and the specific steps of obtaining the intelligent multivariable controller of the aero-engine are as follows:
step B2.1), respectively obtaining the initial control laws of the rotating speed and the pressure ratio according to the NARMA-L2 model obtained by the training of the step A2):
wherein, y1 *[k]And y2 *[k]Respectively outputting expected rotation speed of a high-pressure rotor of the engine and the overall pressure ratio; f. of1(k)、g1(k)、f2(k) And g2(k) The expression of (a) is:
f1(k)=f1(y1[k-1],y1[k-2],…,y1[k-n+1],u1[k-1],...,u1[k-n+1])
g1(k)=g1(y1[k-1],y1[k-2],...,y1[k-n+1],u1[k-1],...,u1[k-n+1])
f2(k)=f2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])
g2(k)=g2(y2[k-1],y2[k-2],...,y2[k-n+1],u2[k-1],...,u2[k-n+1])
step B2.2), based on the error between the desired output and the actual output, the u is compensated for by using the NARMA-L2 linear compensation strategy1[k]And u2[k]Compensating and modifying the neural network f according to a gradient descent algorithm1、g1、f2And g2The topology parameter of (2).
5. The intelligent multivariable control method for the aero-engine with online correction as claimed in claim 4, wherein: said step B2.2) of using a NARMA-L2 linear compensation strategy for u depending on the error between the desired output and the actual output1[k]And u2[k]Compensating and modifying the neural network f according to a gradient descent algorithm1、g1、f2And g2The topology parameters are as follows:
step B2.2.1), judging whether the current time is the initial time, if so, not compensating; if the current time is not the initial time, calculating a compensation value delta u [ k ] of the initial control quantity at the current time according to the actual input u [ k-1] and the actual output y [ k-1] of the engine at the previous time and the neural network outputs f (k-1) and g (k-1) in the NARMA-L2 model, wherein the specific formula is as follows:
α is compensation gain, and is determined by a multi-trial method;
step B2.2.2), obtaining the control quantity u according to the initial control law at the current timeoriginal[k]Compensation value delta u k with control quantity]And calculating to obtain the actual control quantity u [ k ]]The concrete formula is as follows:
u1[k]=u1_original[k]+Δu1[k]
u2[k]=u2_original[k]+Δu2[k]
step B2.2.3) applying the actual control quantity to the aircraft engine to obtain the actual output y [ k ] at the current moment]Using the error e between the desired output and the actual outputC[k]For the neural network group f1、g1、f2And g2Topology ofThe parameters are corrected on line, and the specific formula is as follows:
W(k+1)=W(k)+βΔW(k)
V(k+1)=V(k)+βΔV(k)
wherein, W and V are weight and threshold of the neural network group respectively, β is learning rate, the adjusting speed of the size influence parameter is determined by a multi-trial method, Δ W (k) and Δ V (k) are Jacobian matrixes which respectively represent the increment of the weight and the threshold, and the updating of each element in the topological vector at the moment k is as follows:
wi(k+1)=wi(k)-βΔwi(k)
vi(k+1)=vi(k)-βΔvi(k)
wherein, T ═ y (k), y (k-1), y (k-2),. ·, y (k-n), u (k-1), u (k-2), …, u (k-n) }; w is ai、viThe ith elements representing W and V, respectively; g0Denotes g1Or g2;Δwi、ΔviThe i-th elements respectively representing Δ w (k) and Δ v (k);
obtaining a rotating speed and pressure ratio controller at the k +1 moment according to the updated NARMA-L2 model; to this end, the intelligent multivariable controller design ends.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911132301.9A CN110985216B (en) | 2019-11-19 | 2019-11-19 | Intelligent multivariable control method for aero-engine with online correction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911132301.9A CN110985216B (en) | 2019-11-19 | 2019-11-19 | Intelligent multivariable control method for aero-engine with online correction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110985216A true CN110985216A (en) | 2020-04-10 |
CN110985216B CN110985216B (en) | 2022-06-28 |
Family
ID=70084827
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911132301.9A Active CN110985216B (en) | 2019-11-19 | 2019-11-19 | Intelligent multivariable control method for aero-engine with online correction |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110985216B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112286047A (en) * | 2020-09-28 | 2021-01-29 | 南京航空航天大学 | NARMA-L2 multivariable control method based on neural network |
CN113485117A (en) * | 2021-07-28 | 2021-10-08 | 沈阳航空航天大学 | Multivariable reinforcement learning control method for aircraft engine based on input and output information |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050038532A1 (en) * | 2003-08-12 | 2005-02-17 | Cheng George Shu-Xing | Apparatus and method of controlling single-input-multi-output systems |
CA2762484A1 (en) * | 2009-05-29 | 2010-12-02 | Qingsheng Q. Zheng | Apparatus and method for model quality estimation and model adaptation in multivariable process control |
US20160170384A1 (en) * | 2014-12-11 | 2016-06-16 | University Of New Brunswick | Model predictive controller and method with correction parameter to compensate for time lag |
CN106647253A (en) * | 2016-09-28 | 2017-05-10 | 南京航空航天大学 | Aero-engine distributed control system multi-performance robust tracking control method |
CN108762089A (en) * | 2018-06-15 | 2018-11-06 | 大连理工大学 | A kind of aero-engine on-line optimization and multivariable Control design method based on model prediction |
CN109162813A (en) * | 2018-10-09 | 2019-01-08 | 南京航空航天大学 | One kind being based on the modified Aeroengine Smart method for controlling number of revolution of iterative learning |
-
2019
- 2019-11-19 CN CN201911132301.9A patent/CN110985216B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050038532A1 (en) * | 2003-08-12 | 2005-02-17 | Cheng George Shu-Xing | Apparatus and method of controlling single-input-multi-output systems |
CA2762484A1 (en) * | 2009-05-29 | 2010-12-02 | Qingsheng Q. Zheng | Apparatus and method for model quality estimation and model adaptation in multivariable process control |
US20160170384A1 (en) * | 2014-12-11 | 2016-06-16 | University Of New Brunswick | Model predictive controller and method with correction parameter to compensate for time lag |
CN106647253A (en) * | 2016-09-28 | 2017-05-10 | 南京航空航天大学 | Aero-engine distributed control system multi-performance robust tracking control method |
CN108762089A (en) * | 2018-06-15 | 2018-11-06 | 大连理工大学 | A kind of aero-engine on-line optimization and multivariable Control design method based on model prediction |
CN109162813A (en) * | 2018-10-09 | 2019-01-08 | 南京航空航天大学 | One kind being based on the modified Aeroengine Smart method for controlling number of revolution of iterative learning |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112286047A (en) * | 2020-09-28 | 2021-01-29 | 南京航空航天大学 | NARMA-L2 multivariable control method based on neural network |
WO2022037157A1 (en) * | 2020-09-28 | 2022-02-24 | 南京航空航天大学 | Narma-l2 multi-variable control method based on neural network |
CN112286047B (en) * | 2020-09-28 | 2022-04-08 | 南京航空航天大学 | NARMA-L2 multivariable control method based on neural network |
CN113485117A (en) * | 2021-07-28 | 2021-10-08 | 沈阳航空航天大学 | Multivariable reinforcement learning control method for aircraft engine based on input and output information |
CN113485117B (en) * | 2021-07-28 | 2024-03-15 | 沈阳航空航天大学 | Multi-variable reinforcement learning control method for aeroengine based on input and output information |
Also Published As
Publication number | Publication date |
---|---|
CN110985216B (en) | 2022-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108762089B (en) | Model prediction-based aero-engine online optimization and multivariable control design method | |
CN109162813B (en) | One kind being based on the modified Aeroengine Smart method for controlling number of revolution of iterative learning | |
US11454177B2 (en) | Method of aero-engine on-line optimization and multivariable control based on model prediction | |
WO2019144337A1 (en) | Deep-learning algorithm-based self-adaptive correction method for full-envelope model of aero-engine | |
CN112286047B (en) | NARMA-L2 multivariable control method based on neural network | |
CN110502840B (en) | Online prediction method for gas circuit parameters of aero-engine | |
CN108828947B (en) | Modeling method for time-lag-containing uncertain fuzzy dynamic model of aircraft engine | |
CN106483850A (en) | The Fuzzy Self-adaptive PID method for designing that a kind of aero-engine is feedovered based on RBF neural | |
CN110219736B (en) | Aero-engine direct thrust control method based on nonlinear model predictive control | |
CN111042928B (en) | Variable cycle engine intelligent control method based on dynamic neural network | |
CN109800449B (en) | Neural network-based aeroengine compression component characteristic correction method | |
CN111006843B (en) | Continuous variable speed pressure method of temporary impulse type supersonic wind tunnel | |
CN110985216B (en) | Intelligent multivariable control method for aero-engine with online correction | |
CN110221537A (en) | Control method, device and aero-engine limitation protection control method, device | |
Zheng et al. | A study of aero-engine control method based on deep reinforcement learning | |
CN111680357A (en) | Component-level non-iterative construction method of variable-cycle engine airborne real-time model | |
CN107357176B (en) | Modeling method for test run data of aero-engine | |
CN113267314A (en) | Supersonic flow field total pressure control system of temporary-impulse wind tunnel | |
CN111255574A (en) | Autonomous control method for thrust recession under inlet distortion of aircraft engine | |
Yangjing et al. | Neural network-based model predictive control with fuzzy-SQP optimization for direct thrust control of turbofan engine | |
Sun et al. | Event-triggered intelligent critic control with input constraints applied to a nonlinear aeroelastic system | |
CN106597847B (en) | Maneuvering load controller based on recurrent neural network and control method thereof | |
CN114527654A (en) | Turbofan engine direct thrust intelligent control method based on reinforcement learning | |
CN111856929B (en) | Two-degree-of-freedom H-infinity controller for fault-tolerant gain scheduling of aero-engine | |
CN114815616A (en) | Intelligent regulation and control method and system for mode conversion of turbine stamping combined type engine |
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 |