Online identification method for operating efficiency of swing spray pipe
Technical Field
The invention relates to an on-line identification method for the operating efficiency of a swing spray pipe, and belongs to the field of flight control.
Background
The swing spray pipe generates a control moment by deviating the thrust from the axis of the aircraft, the control moment is mainly influenced by the swing angle and the thrust of the spray pipe, and the overall parameters such as mass center, moment of inertia and the like can also influence the spray pipe operating efficiency. Under different engine working temperature states, the engine thrust and working time are different, and meanwhile, the total parameters such as mass, mass center and rotational inertia are strong in discretization along with time change, so that the swing spray pipe operation in the flight process needs to be identified, and the control quality of the whole flight stage is improved.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, and designs an on-line identification method for the operating efficiency of the swing spray pipe by analyzing the relevance between the overall parameters and the flight state of the aircraft based on the theoretical flight trajectory, so that the adaptability of the aircraft to the discrete difference of the operating efficiency of the swing spray pipe and the flight control quality are improved.
The invention adopts the technical scheme that:
An on-line identification method for the operating efficiency of a swing spray pipe comprises the following steps:
step one, calculating a theoretical flight track according to engine working thrust curves under different working temperature states obtained by engine ground test;
determining the relation between the overall parameters of the aircraft in the flight process and key flight state quantities according to the obtained theoretical flight track, wherein the overall parameters of the aircraft comprise the mass, mass center, moment of inertia and thrust of the aircraft, and the key flight state quantities comprise flight altitude, mach number, flight attack angle, dynamic pressure and acceleration;
Fitting a function of the overall parameters of the aircraft along with the change of the key flight state quantity according to the relation between the overall parameters of the aircraft and the key flight state quantity;
calculating the change rule of the elastic modal frequency of the aircraft under different working temperature states, and designing a filter according to the change rule to filter overload information measured in the flight process;
And fifthly, establishing a mathematical model of the swing spray pipe operating efficiency changing along with the key flight state quantity, so as to realize online identification of the swing spray pipe operating efficiency.
Further, the theoretical flight trajectory includes changes in flight altitude, mach number, angle of attack, dynamic pressure, mass and moment of inertia during flight.
Further, the determining the relationship between the overall parameters of the aircraft and the critical flight state quantity in the flight process according to the obtained theoretical flight trajectory specifically comprises:
and respectively establishing fitting functions of the flight state quantity and the overall parameters of the aircraft, wherein the fitting functions are linear functions, comparing the fitting functions, and selecting the flight state quantity with the highest linearity and fitting degree with the overall parameters of the aircraft as a key flight state quantity.
Further, according to the relationship between the overall parameters of the aircraft and the critical flight state quantity, a function of the overall parameters of the aircraft, which changes along with the critical flight state quantity, is fitted, specifically:
Xg=k1×Ma+m1
Jz=k2×Ma+m2
Mass=k3×Ma+m3
Ca=k4×Ma+m4
Wherein the aircraft centroid Xg, moment of inertia Jz, mass, and drag coefficient Ca are fitted by mach number Ma, where k i、mi is the result of the parameters of the fit, i=1, 2,3,4.
Further, according to the overall parameters and the resistance coefficient obtained by the function of the overall parameters of the aircraft along with the change of the key flight state quantity, the Thrust is obtained through calculation by combining the real axial acceleration Ax obtained by the measurement of the aircraft.
Further, calculating the change rule of the elastic modal frequency of the aircraft under different working temperature states, and designing a filter based on the change rule to filter the real axial acceleration Ax measured in the flight process, wherein the method specifically comprises the following steps:
design of notch filter G(s):
Wherein T 1、T2 is a filtering time constant, and ζ 1、ζ2 is a damping coefficient of the structural filter;
The filter time constant T 1、T2 is adjusted by the Mach number Ma to make the change of the filter center frequency consistent with the change of the elastic mode frequency.
Further, the filter time constant T 1、T2:
T1=1/(m5+k5×Ma)
T2=1/(m6+k6×Ma)
Wherein m 5、m6、k5、k6 is a parameter result obtained by fitting.
Further, a mathematical model of the swing nozzle manipulation efficiency changing along with the key flight state quantity is established, so that the on-line identification of the swing nozzle manipulation efficiency a 3bp is realized, specifically:
Thrust=Ca×Q×S+Mass×Ax
a3bp=Thrust×(L-Xg)/Jz
wherein S is the characteristic area of the aircraft, L is the length of the aircraft, and Q is the dynamic pressure in real time in flight.
In a second aspect, the present invention also proposes a non-volatile storage medium comprising a computer program product which, when executed, performs the swing nozzle manipulation efficiency online identification method.
In a third aspect, the present invention also proposes a computer program product which, when executed by a processor, implements the steps of the method for online identification of operating efficiency of a swinging nozzle.
Compared with the prior art, the invention has the beneficial effects that:
(1) According to the invention, the flight state quantity with strong correlation of the overall parameters such as the thrust, the mass center, the rotational inertia and the like of the engine is analyzed and confirmed according to the theoretical flight trajectory, a mathematical model of the swing spray pipe operating efficiency changing along with the key flight state quantity is established, and a series of problems caused by the thrust characteristics of the aircraft engine under different temperature work are solved by on-line identification of the influence of the swing spray pipe operating efficiency on the improvement of the thrust discrete characteristics.
(2) According to the invention, theoretical flight trajectories under different engine working states are analyzed, the selected key flight state quantity is used as a basis for identification, the physical concept is clear, and in order to restrain the influence of elastic vibration in the flight process, a filter is designed according to the rule that the elastic mode is regulated along with the key flight state, so that the adaptability to the elastic vibration characteristics of the engine under different working states is good.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a diagram of theoretical flight trajectories of an aircraft under different engine operating conditions calculated from a flight mission;
FIG. 3 is a schematic diagram of Mach number versus overall aircraft parameters.
Detailed Description
The invention discloses an on-line identification method for the operating efficiency of a swing jet pipe by using a solid rocket engine, which is characterized in that a theoretical flight track is calculated according to the ground test result of the engine, the relevance of overall parameters such as the thrust, the mass center, the rotational inertia and the like of an aircraft under the condition of high and low normal temperature and different flight states is analyzed, so that a mathematical model of the operating efficiency of the swing jet pipe changing along with the key flight state quantity is established, and the on-line identification of the operating efficiency is completed.
As shown in FIG. 1, the invention provides an on-line identification method for the operating efficiency of a swing spray pipe, which comprises the following steps:
step one, calculating a theoretical flight track according to engine working thrust curves under different working temperature states obtained by engine ground test;
The theoretical flight trajectory comprises the changes of flight altitude, mach number, attack angle, dynamic pressure, mass and rotational inertia in the flight process.
Determining the relation between the overall parameters of the aircraft in the flight process and key flight state quantities according to the obtained theoretical flight track, wherein the overall parameters of the aircraft comprise the mass, mass center, moment of inertia and thrust of the aircraft, and the key flight state quantities comprise flight altitude, mach number, flight attack angle, dynamic pressure and acceleration;
And respectively establishing fitting functions of the flight state quantity and the overall parameters of the aircraft, wherein the fitting functions are linear functions, comparing the fitting functions, and selecting the flight state quantity with the highest linearity and fitting degree with the overall parameters of the aircraft as a key quantity.
Fitting a function of the overall parameters of the aircraft along with the change of the key flight state quantity according to the relation between the overall parameters of the aircraft and the key flight state quantity, wherein the function comprises the following specific steps:
Xg=k1×Ma+m1
Jz=k2×Ma+m2
Mass=k3×Ma+m3
Ca=k4×Ma+m4
Wherein the aircraft centroid Xg, moment of inertia Jz, mass, and drag coefficient Ca are fitted by mach number Ma, where k i、mi is the result of the parameters of the fit, i=1, 2,3,4.
And according to the overall parameters and the resistance coefficient obtained by the function of the overall parameters of the aircraft along with the change of the key flight state quantity, combining the real axial acceleration Ax measured by the aircraft, and calculating to obtain the Thrust.
Calculating the change rule of the elastic modal frequency of the aircraft under different working temperature states, and designing a filter according to the change rule to filter overload information measured in the flight process;
The filter is designed according to the method to carry out filter processing on the real axial acceleration Ax measured in the flying process, and the method specifically comprises the following steps:
design of notch filter G(s):
Wherein T 1、T2 is a filtering time constant, and ζ 1、ζ2 is a damping coefficient of the structural filter;
The filter time constant T 1、T2 is adjusted by the Mach number Ma to make the change of the filter center frequency consistent with the change of the elastic mode frequency.
Filter time constant T 1、T2:
T1=1/(m5+k5×Ma)
T2=1/(m6+k6×Ma)
Wherein m 5、m6、k5、k6 is a parameter result obtained by fitting.
And fifthly, establishing a mathematical model of the swing spray pipe operating efficiency changing along with the key flight state quantity, so as to realize online identification of the swing spray pipe operating efficiency.
The on-line identification of the swing spray pipe operating efficiency a 3bp is realized, and specifically:
Thrust=Ca×Q×S+Mass×Ax
a3bp=Thrust×(L-Xg)/Jz
wherein S is the characteristic area of the aircraft, L is the length of the aircraft, and Q is the dynamic pressure in real time in flight.
Examples:
according to the flight tasks, theoretical flight trajectories under different engine working temperature states are calculated respectively, as shown in fig. 2, and each piece of flight trajectory data comprises the change condition of data such as flight altitude, mach number, attack angle, dynamic pressure, mass, rotational inertia and the like in the flight process.
By analyzing the theoretical flight trajectory, the Mach number is selected as a key flight state quantity in the embodiment, as shown in fig. 3, the correlation between the thrust and the Mach number is poor, and the correlation between the aircraft mass, the centroid, the moment of inertia, the resistance and the Mach number is strong and the linearity is good in the embodiment. According to the dynamic relationship of the projectile, the axial resultant force born by the aircraft can be similar to the resultant force of the thrust and the resistance of the engine, so that the axial acceleration information measured by the aircraft is also required to be introduced to participate in the identification of the operating efficiency of the swing jet pipe.
The aircraft centroid Xg, moment of inertia Jz, mass and drag coefficient Ca are fitted using mach number Ma as shown in equation 1. Where k i(i=1,2,3,4…),mi (i=1, 2,3, 4.) is the result of the parameters obtained by fitting.
Xg=k1×Ma+m1
Jz=k2×Ma+m2
Mass=k3×Ma+m3
Ca=k4×Ma+m4 (1)
According to the overall parameters and the resistance coefficient obtained in the formula 1, the Thrust can be obtained through calculation according to the real axial acceleration Ax measured by the aircraft. Since there may be elastic vibration in the axial acceleration Ax obtained by actual measurement, if it is not filtered, it may affect the accuracy of a 3bp identification, so a notch filter is designed, and the expression of the filter is as shown in the formula (2).
Since the elastic mode frequency is directly related to the Mass of the aircraft, and the relation between Mass and Ma is shown in formula 1, the filter time constant can be adjusted by Ma, so that the change of the filter center frequency is consistent with the change of the elastic mode frequency.
Wherein T 1、T2 is a filtering time constant, and ζ 1、ζ2 is a damping coefficient of the structural filter.
T1=1/(m5+k5×Ma)
T2=1/(m6+k6×Ma) (3)
And then the identification of the swing nozzle manipulation efficiency a 3bp is completed as shown in formula 4.
Thrust=Ca×Q×S+Mass×Ax
a3bp=Thrust×(L-Xg)/Jz (4)
Wherein S is the characteristic area of the aircraft, L is the length of the aircraft, and Q is the dynamic pressure in real time in flight.
Although the present invention has been described in terms of the preferred embodiments, it is not intended to be limited to the embodiments, and any person skilled in the art can make any possible variations and modifications to the technical solution of the present invention by using the methods and technical matters disclosed above without departing from the spirit and scope of the present invention, so any simple modifications, equivalent variations and modifications to the embodiments described above according to the technical matters of the present invention are within the scope of the technical matters of the present invention.
The invention is not described in detail in the field of technical personnel common knowledge.