CN112046761B - Airplane icing on-line detection method based on statistical test and filtering - Google Patents

Airplane icing on-line detection method based on statistical test and filtering Download PDF

Info

Publication number
CN112046761B
CN112046761B CN202010771628.7A CN202010771628A CN112046761B CN 112046761 B CN112046761 B CN 112046761B CN 202010771628 A CN202010771628 A CN 202010771628A CN 112046761 B CN112046761 B CN 112046761B
Authority
CN
China
Prior art keywords
icing
data
aircraft
flight
overload
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010771628.7A
Other languages
Chinese (zh)
Other versions
CN112046761A (en
Inventor
丁娣
何开锋
汪清
钱炜祺
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202010771628.7A priority Critical patent/CN112046761B/en
Publication of CN112046761A publication Critical patent/CN112046761A/en
Application granted granted Critical
Publication of CN112046761B publication Critical patent/CN112046761B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64DEQUIPMENT FOR FITTING IN OR TO AIRCRAFT; FLIGHT SUITS; PARACHUTES; ARRANGEMENTS OR MOUNTING OF POWER PLANTS OR PROPULSION TRANSMISSIONS IN AIRCRAFT
    • B64D15/00De-icing or preventing icing on exterior surfaces of aircraft
    • B64D15/20Means for detecting icing or initiating de-icing

Abstract

The invention discloses an on-line detection method for aircraft icing based on statistical test and filtering, which at least comprises the following steps: s1: acquiring flight state measurement data, engine thrust data and rudder deflection input data in the flight process; s2: detecting the icing starting moment by using a generalized likelihood ratio test method; s3: generating a disturbance signal superimposed on the rudder deflection input data while using HThe filtering method carries out combined state estimation and icing-influenced pneumatic derivative identification on the flight state measurement data; s4: detecting the icing finishing time by using a generalized likelihood ratio test method; s5: stopping generating and superposing rudder deviation disturbance signals and continuously utilizing HThe filtering carries out joint state estimation and icing-influenced pneumatic derivative identification on the measured data until the flight is finished. By combining the rapid icing detection method with the parameter estimation method, the aircraft icing detection algorithm with stronger function is obtained.

Description

Airplane icing on-line detection method based on statistical test and filtering
Technical Field
The invention belongs to the technical field of online detection of aircraft icing, and particularly relates to an online detection method of aircraft icing based on statistical inspection and filtering.
Background
The online detection of the airplane icing means that the icing phenomenon generated in the flying process of the airplane is detected through a detection device or a detection method, and an icing alarm is given to start an anti-icing and deicing device of the airplane, so that the harm of the icing in the flying process to the safety of the airplane is avoided.
The icing on-line detection generally comprises two categories, namely hardware detection and algorithm detection, wherein the hardware detection is to detect icing parts, icing type parameters and the like through hardware devices such as sensors arranged on the surface of the airplane body, and the algorithm detection is to utilize a mathematical method and a physical model to carry out on-line estimation on the aerodynamic characteristic parameters of the airplane affected by icing through flight measurement data so as to obtain the key parameter information of the airplane icing. Compared with algorithm detection, hardware detection needs to arrange a large number of sensing devices on a machine body, so that the development cost of an airplane is generally increased, the problem that the sensor devices are not suitable to be arranged or are not sufficiently arranged exists, and the detection effect of icing hardware is greatly influenced; the online detection of the algorithm can greatly save development cost, the algorithm detection gets rid of the limitation of a sensor, the icing abnormity is found by estimating the characteristic parameters of the airplane body, and the online detection method has better adaptability and flexibility.
Aircraft icing algorithm detection is generally focused on the study of detection algorithms, and is less studied for some specific problems applied to flight procedures. Such as how to ensure the accuracy of the on-line estimation of parameters during flight, how to detect icing, how to add input disturbance signals at appropriate times, etc. Therefore, the icing on-line detection research is needed to be developed based on the actual flight process, the specific problems in algorithm detection are solved, and the airplane icing algorithm detection capability is improved.
Disclosure of Invention
The invention aims to provide a feasible aircraft icing on-line detection method based on statistical test and filtering, aiming at overcoming the technical problem of aircraft icing algorithm detection in actual flight application in the prior art. By combining the rapid icing detection method with the parameter estimation method, a powerful airplane icing detection algorithm is obtained; by adding the input disturbance signal at a proper time, a more accurate aerodynamic characteristic parameter estimation result is obtained under the condition of ensuring that the influence of the input disturbance on an aircraft flight system is small.
The purpose of the invention is realized by the following technical scheme:
an aircraft icing on-line detection method based on statistical test and filtering is characterized by at least comprising the following steps:
s1: acquiring flight state measurement data, engine thrust data and rudder deflection input data in the flight process;
s2: detecting the icing starting moment by using a generalized likelihood ratio test method;
s3: generating a disturbance signal superimposed on the rudder deflection input data while using HThe filtering method carries out combined state estimation and icing-influenced pneumatic derivative identification on the flight state measurement data;
s4: detecting the icing finishing time by using a generalized likelihood ratio test method;
s5: stopping generating and superposing rudder deviation disturbance signals and continuously utilizing HThe filtering carries out joint state estimation and icing-influenced pneumatic derivative identification on the measured data until the flight is finished.
According to a preferred embodiment, in step S1, the flight state measurement data at least includes position, speed, overload, angular speed, attitude angle and dynamic pressure data of the aircraft.
According to a preferred embodiment, said step S2 includes: sequentially extracting data segments from flight state measurement data to calculate overload margins, and judging whether the statistical characteristics of each segment of overload margins exceed a given threshold value by using a generalized likelihood ratio test method to judge the margins rzWhether or not an offset occurs, when rzAnd judging that icing has occurred if the deviation except the random noise exists, and taking the end time of the section of data as the icing start time.
According to a preferred embodiment, the calculation of the overload margin in step S2 includes: extracting a data segment of a given detection window length from the data collected in step S1, calculating an aircraft normal overload margin r from the difference between the estimated overload and the measured overloadz
Figure BDA0002616840530000021
Wherein N iszmIn order to measure the resulting normal overload,
Figure BDA0002616840530000022
to estimate the normal overload, and,
Figure BDA0002616840530000023
m is the aircraft mass, S is the aircraft characteristic area, g0Is the gravity acceleration constant, q, corresponding to sea levelFor measuring dynamic pressure, alpha is for measuring angle of attack, PzIn order to measure the normal thrust force,
Figure BDA0002616840530000024
and
Figure BDA0002616840530000025
the coefficient of drag and lift for the aircraft under the ice-free condition is obtained by interpolation of a ground test database or by calculation of a pneumatic derivative in combination with flight state parameters.
According to a preferred embodiment, in step S3, the generated disturbing signal is a "3211" disturbing signal.
According to a preferred embodiment, said step S4 includes: by aerodynamically stable derivative C of pitch directionThe identification result of (2) carries out convergence judgment to give the end time of icing.
According to a preferred embodiment, the step S4 specifically includes: first from HSequentially extracting C with given window length from the identification resultIdentification result using two adjacent CCalculating detection margin r of identification result data segmentCmα
Figure BDA0002616840530000026
i is 1,2, …, n, wherein i represents C extracted from the ith segmentIdentifying the result, reverse utilizing generalized likelihood ratio test, judging r when the margin statistical characteristic is less than or equal to the given threshold valueCmαIf there is no offset, then CThe identification result of (1) is converged, the icing accumulation phase is ended, and the end time of the (i + 1) th data segment is taken as the icing end time.
The main scheme and the further selection schemes can be freely combined to form a plurality of schemes which are all adopted and claimed by the invention; in the invention, the selection (each non-conflict selection) and other selections can be freely combined. The skilled person in the art can understand that there are many combinations, which are all the technical solutions to be protected by the present invention, according to the prior art and the common general knowledge after understanding the scheme of the present invention, and the technical solutions are not exhaustive herein.
The invention has the beneficial effects that:
1. the invention discloses an on-line detection method for aircraft icing, which is a generalized likelihood ratio method andHthe filtering method is combined, so that the change of the pneumatic derivative influenced by icing accumulation can be obtained while icing early warning is obtained, and compared with a single method, the algorithm can obtain more icing related information and has stronger online airplane icing detection capability;
2. the invention solves the problem of determining the action time of the input disturbance signal in the flight process, and the input disturbance signal is very important for ensuring the identification precision of the pneumatic derivative of the airplane in the icing accumulation process, so the solution provided by the invention can superpose the small disturbance signal of '3211' particularly in the icing accumulation process, can avoid unnecessary or meaningless input disturbance, and can reduce the disturbance of the disturbance signal to the flight system as much as possible.
3. The invention provides a method for judging the end of icing, which can obtain the conclusion whether icing is ended or not by judging whether the identification result of a pneumatic derivative affected by icing is converged or not, wherein the convergence is judged by reversely adopting a generalized likelihood ratio test method, and the method provides a basis for accurately applying a small-disturbance input signal.
Drawings
FIG. 1 is a flow chart of an online detection method for aircraft icing according to the present invention;
FIG. 2 is a diagram of the results of the on-line detection of icing on a flat aircraft wing in an embodiment of the invention.
Detailed Description
The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention. It is to be noted that the features in the following embodiments and examples may be combined with each other without conflict.
It should be noted that, 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 are clearly and completely described below, and it is obvious that the described embodiments are some embodiments of the present invention, but not all embodiments.
Thus, the following detailed description of the embodiments of the present invention is not intended to limit the scope of the invention as claimed, but is merely representative of selected embodiments of the invention. 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.
In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc. indicate orientations and positional relationships that are conventionally used in the products of the present invention, and are used merely for convenience in describing the present invention and for simplicity in description, but do not indicate or imply that the devices or elements referred to must have a particular orientation, be constructed in a particular orientation, and be operated, and therefore, should not be construed as limiting the present invention. Furthermore, the terms "first," "second," "third," and the like are used solely to distinguish one from another and are not to be construed as indicating or implying relative importance.
Furthermore, the terms "horizontal", "vertical", "overhang" and the like do not imply that the components are required to be absolutely horizontal or overhang, but may be slightly inclined. For example, "horizontal" merely means that the direction is more horizontal than "vertical" and does not mean that the structure must be perfectly horizontal, but may be slightly inclined.
In the description of the present invention, it should also be noted that, unless otherwise explicitly specified or limited, the terms "disposed," "mounted," "connected," and "connected" are to be construed broadly and may, for example, be fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
In addition, it should be noted that, in the present invention, if the specific structures, connection relationships, position relationships, power source relationships, and the like are not written in particular, the structures, connection relationships, position relationships, power source relationships, and the like related to the present invention can be known by those skilled in the art without creative work on the basis of the prior art.
Example 1:
referring to fig. 1, the invention discloses an aircraft icing on-line detection method based on statistical test and filtering, which is characterized by at least comprising the following steps:
step S1: and acquiring flight state measurement data, engine thrust data and rudder deflection input data in the flight process.
Preferably, in step S1, the flight state measurement data at least includes position, speed, overload, angular speed, attitude angle and dynamic pressure data of the aircraft.
Step S2: and detecting the icing starting moment by using a generalized likelihood ratio test method.
Preferably, the step S2 includes: sequentially extracting data segments from flight state measurement data to calculate overload margins, and judging whether the statistical characteristics of each segment of overload margins exceed a given threshold value by using a generalized likelihood ratio test method to judge the margins rzWhether or not an offset occurs, when rzAnd judging that icing has occurred if the deviation except the random noise exists, and taking the end time of the section of data as the icing start time.
Further, the calculation process of the overload margin in step S2 includes: extracting a data segment of a given detection window length from the data collected in step S1, calculating an aircraft normal overload margin r from the difference between the estimated overload and the measured overloadz
Figure BDA0002616840530000051
Wherein N iszmIn order to measure the resulting normal overload,
Figure BDA0002616840530000052
to estimate normal overload, and,
Figure BDA0002616840530000053
m is the aircraft mass, S is the aircraft characteristic area, g0Is the gravity acceleration constant, q, corresponding to sea levelFor measuring dynamic pressure, alpha is for measuring angle of attack, PzIn order to measure the normal thrust force,
Figure BDA0002616840530000054
and
Figure BDA0002616840530000055
the coefficient of drag and lift for the aircraft under the ice-free condition is obtained by interpolation of a ground test database or by calculation of a pneumatic derivative in combination with flight state parameters.
Step S3: generating a disturbance signal superimposed on the rudder deflection input data while using HThe filtering method carries out combined state estimation and icing-influenced pneumatic derivative identification on the flight state measurement data.
Preferably, in the step S3, the generated disturbance signal is a "3211" disturbance signal.
The method solves the problem of determining the action time of the input disturbance signal in the flight process, and the input disturbance signal is very important for ensuring the identification precision of the pneumatic derivative of the airplane in the icing accumulation process.
Step S4: and detecting the icing ending moment by using a generalized likelihood ratio test method.
Preferably, the step S4 includes: by aerodynamically stable derivative C of pitch directionThe identification result of (2) carries out convergence judgment to give the end time of icing.
Further, the step S4 specifically includes: first from HSequentially extracting C with given window length from the identification resultIdentification result using two adjacent CIdentification result data segment calculationDetection margin rCmα
Figure BDA0002616840530000056
i is 1,2, …, n, wherein i represents C extracted from the ith segmentIdentifying the result, reverse utilizing generalized likelihood ratio test, judging r when the margin statistical characteristic is less than or equal to the given threshold valueCmαIf there is no offset, then CThe identification result of (1) is converged, the icing accumulation phase is ended, and the end time of the (i + 1) th data segment is taken as the icing end time.
Namely, the conclusion whether icing is finished or not can be obtained by judging whether the identification result of the pneumatic derivative affected by icing is converged or not, wherein the convergence is judged by reversely adopting a generalized likelihood ratio test method, and the method provides a basis for accurately applying a small-disturbance input signal.
Step S5: stopping generating and superposing rudder deviation disturbance signals and continuously utilizing HThe filtering carries out joint state estimation and icing-influenced pneumatic derivative identification on the measured data until the flight is finished.
Example 1
The method is used for carrying out icing on-line detection aiming at the longitudinal flight simulation data of a horizontal cruise section of a certain airplane in consideration of the icing accumulation process at the leading edge of the wing, and the specific flow is shown in figure 1.
The cruising altitude of the airplane is 5000m, the cruising speed is 0.3 Mach, the trim elevator deflection angle is 1.06 degrees, the attitude angle and the angular rate are kept to be zero in the cruising process, and if the leading edge of the wing begins to freeze 100s after the simulation begins, the icing accumulation time is Tcld200s, the cumulative velocity is given byice(Tcld) 1.0 and ηice(TcldAnd/2) is described as 0.7.
The concrete implementation steps of the calculation example are as follows:
the method comprises the following steps: flight data required to be acquired for online detection are acquired according to simulation of a longitudinal flight dynamics model of the airplane, change of aerodynamic derivatives under icing of wings is considered, and Gaussian white noise is added to a simulation result to serve as measurement noise.
Step two: will simulateThe sampling frequency of the data is 10ms, the length of a detection window is 5s, namely 500 data points are sequentially extracted from the simulation data, and the sampling frequency is utilized
Figure BDA0002616840530000061
Estimating normal overload N for each sampling pointzIs further compared with NzSubtracting the measured values to obtain the normal overload allowance r corresponding to each sampling pointz
Step three: calculating mean value of overload margin by using generalized likelihood ratio method
Figure BDA0002616840530000062
Sum mean square error
Figure BDA0002616840530000063
To obtain a statistical test quantity Tx
Figure BDA0002616840530000064
Where N is the number of data points, the false alarm probability PFAThe small amount is 1e-6, and the judgment threshold gamma is 23.93 which can be calculated by chi-square distribution.
The method of the invention combines a generalized likelihood ratio method with HThe filtering method is combined, the change of the pneumatic derivative influenced by icing accumulation can be obtained while icing early warning is obtained, and compared with a single method, the algorithm can obtain more icing related information and has stronger airplane icing online detection capability.
Step four: the statistical test quantity T of each extracted data segmentxComparing with a decision threshold gamma if TxIf the value is less than or equal to gamma, returning to the step two, and considering that no icing phenomenon occurs, the pneumatic derivative affected by icing is not identified, and the pneumatic derivative of the airplane is the same as the non-icing state and is kept unchanged; if TxAnd if the data segment is more than gamma, judging that the icing is started, and taking the end time of the data segment as the icing start time.
Step five: when the icing is judged to be started, a small disturbance signal of '3211' with the amplitude of 5 degrees and the period of 5s is generated and is superposed on the matchingWhen the flight simulation result is input to an elevator rudder deflection signal in a flat mode, due to the influence of rudder deflection disturbance, the flight simulation result can generate periodic oscillation so as to be convenient for identifying the change of a pneumatic derivative; by simultaneous use of HThe filtering method carries out state estimation and pneumatic derivative identification on the flight simulation output data, and obtains the state variable and C of the airplane by reading in the flight simulation output data point by point、CAnd CEstimation of three aerodynamically stable derivatives that are more severely affected by icing; after judging the onset of icing, it is also necessary to aim at CCarrying out convergence judgment detection on the identification result to determine the icing finishing time, setting the length of a detection window to be 5s, and setting 500CSequentially extracting identification data, subtracting the front and rear data sections to obtain detection allowance, and obtaining the statistical detection quantity T of each section of detection allowance by using the method in step threex
Step six: the statistical detection quantity T of the identification result marginxComparing with a decision threshold gamma if TxIf the deviation is larger than gamma, returning to the step five, and if the deviation is still in the icing accumulation stage, continuously superposing the small rudder disturbance signal and identifying the pneumatic derivative; if TxIf gamma is less than or equal to gamma, then C is judgedThe data segment is taken as the icing ending time when the icing accumulation process is ended and the ending time of the data segment is taken as the icing ending time, the generation of the small disturbance rudder deflection signal is stopped, and H is continuously utilizedAnd (4) filtering and estimating a system state variable and identifying three aerodynamic stability derivatives until the flight process is finished.
Namely, the conclusion whether icing is finished or not can be obtained by judging whether the identification result of the pneumatic derivative affected by icing is converged or not, wherein the convergence is judged by reversely adopting a generalized likelihood ratio test method, and the method provides a basis for accurately applying a small-disturbance input signal.
Aiming at the wing icing calculation example of the level flight cruise section, the steps are adopted, the icing starting time obtained through final calculation is 174.99s, the icing finishing time is 304.99s, and the C is divided by three aerodynamic stability derivatives obtained through identificationThe oscillations are more pronounced due to measurement noise, CAnd CRoot mean square of relative truth of identification resultThe errors are both small, 8.15% and 7.78% of the true value, respectively, and follow the ice accumulation process well, as shown in fig. 2, where: etaiceIs the pneumatic derivative normalization result and represents the icing severity coefficient, t is time, CFor stable derivatives of lift, CFor steady derivative of pitching moment, tice0For the onset of icing, tice1Rms represents the root mean square error for the end of icing time. The result shows that the method has the capability of detecting the beginning and the end of icing and the capability of tracking the change of the aerodynamic derivative of the icing accumulation process, and can be conveniently popularized to similar online detection of the icing in the flight process of the airplane.
The foregoing basic embodiments of the invention and their various further alternatives can be freely combined to form multiple embodiments, all of which are contemplated and claimed herein. In the scheme of the invention, each selection example can be combined with any other basic example and selection example at will. Numerous combinations will be known to those skilled in the art.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.

Claims (4)

1. An aircraft icing on-line detection method based on statistical test and filtering is characterized by at least comprising the following steps:
s1: acquiring flight state measurement data, engine thrust data and rudder deflection input data in the flight process;
in step S1, the flight status measurement data at least includes position, speed, overload, angular speed, attitude angle and dynamic pressure data of the aircraft
S2: detecting the icing starting moment by using a generalized likelihood ratio test method;
the step S2 includes:
sequentially extracting data segments from flight state measurement data to calculate overload margin, and judging by using generalized likelihood ratio test methodJudging whether the statistical characteristic of each overload margin exceeds a given threshold value to judge the margin rzWhether or not a shift has occurred is determined,
when r iszJudging that icing occurs if offset except random noise exists, and taking the ending time of the section of data as the icing starting time;
the calculation process of the overload margin in step S2 includes:
extracting a data segment of a given detection window length from the data collected in step S1, calculating an aircraft normal overload margin r from the difference between the estimated overload and the measured overloadz
Figure FDA0003240730150000011
Wherein N iszmIn order to measure the resulting normal overload,
Figure FDA0003240730150000012
in order to estimate the normal overload, the overload is estimated,
and the number of the first and second electrodes,
Figure FDA0003240730150000013
m is the aircraft mass, S is the aircraft characteristic area, g0Is the gravity acceleration constant, q, corresponding to sea levelFor measuring dynamic pressure, alpha is for measuring angle of attack, PzIn order to measure the normal thrust force,
Figure FDA0003240730150000014
and
Figure FDA0003240730150000015
the drag and lift coefficients of the aircraft under the ice-free condition are obtained by interpolation of a ground test database or calculation by combining pneumatic derivative with flight state parameters;
s3: generating a disturbance signal superimposed on the rudder deflection input data while using HThe filtering method carries out combined state estimation and icing-influenced pneumatic derivative identification on the flight state measurement data;
s4: detecting the icing finishing time by using a generalized likelihood ratio test method;
s5: stopping generating and superposing rudder deviation disturbance signals and continuously utilizing HThe filtering carries out joint state estimation and icing-influenced pneumatic derivative identification on the measured data until the flight is finished.
2. The method for detecting icing condition of an aircraft according to claim 1, wherein the disturbance signal generated in step S3 is a "3211" disturbance signal.
3. An aircraft icing on-line detection method based on statistical test and filtering according to claim 1, wherein the step S4 includes: by aerodynamically stable derivative C of pitch directionThe identification result of (2) carries out convergence judgment to give the end time of icing.
4. The method for detecting icing on-line of an aircraft based on statistical testing and filtering as claimed in claim 3, wherein said step S4 specifically comprises:
first from HSequentially extracting C with given window length from the identification resultIdentification result using two adjacent CCalculating detection margin r of identification result data segmentCmα
Figure FDA0003240730150000021
Wherein i represents C extracted from the ith segmentIdentifying the result, reverse utilizing generalized likelihood ratio test, judging r when the margin statistical characteristic is less than or equal to the given threshold valueCmαIf there is no offset, then CThe identification result of (1) is converged, the icing accumulation phase is ended, and the end time of the (i + 1) th data segment is taken as the icing end time.
CN202010771628.7A 2020-08-04 2020-08-04 Airplane icing on-line detection method based on statistical test and filtering Active CN112046761B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010771628.7A CN112046761B (en) 2020-08-04 2020-08-04 Airplane icing on-line detection method based on statistical test and filtering

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010771628.7A CN112046761B (en) 2020-08-04 2020-08-04 Airplane icing on-line detection method based on statistical test and filtering

Publications (2)

Publication Number Publication Date
CN112046761A CN112046761A (en) 2020-12-08
CN112046761B true CN112046761B (en) 2021-10-15

Family

ID=73602413

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010771628.7A Active CN112046761B (en) 2020-08-04 2020-08-04 Airplane icing on-line detection method based on statistical test and filtering

Country Status (1)

Country Link
CN (1) CN112046761B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113777931B (en) * 2021-11-09 2022-02-11 中国空气动力研究与发展中心计算空气动力研究所 Icing wing type pneumatic model construction method, device, equipment and medium
CN114625159B (en) * 2022-01-21 2023-07-28 中国空气动力研究与发展中心计算空气动力研究所 Icing aircraft control method based on controlled variable

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5028929A (en) * 1990-04-30 1991-07-02 University Corporation For Atmospheric Research Icing hazard detection for aircraft
US6490510B1 (en) * 1999-04-30 2002-12-03 Thales Avionics S.A. Fixed multifunction probe for aircraft
US6757668B1 (en) * 1999-11-05 2004-06-29 General Electric Company Information fusion of classifiers in systems with partial redundant information
EP2829473A1 (en) * 2013-07-24 2015-01-28 Air China Limited Method and apparatus for monitoring a wing anti-icing valve
CN109911216A (en) * 2017-12-12 2019-06-21 波音公司 The system and method for hitting the position of the limit for correcting the water in aerofoil profile

Family Cites Families (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5195046A (en) * 1989-01-10 1993-03-16 Gerardi Joseph J Method and apparatus for structural integrity monitoring
FR2958622B1 (en) * 2010-04-09 2013-05-03 Airbus Operations Sas METHOD AND SYSTEM FOR CONTROLLING ICE FORMATION ON A FLIGHT AIRCRAFT.
CN102313512B (en) * 2010-07-02 2013-08-07 中国商用飞机有限责任公司 Image icing detection method
GB201014244D0 (en) * 2010-08-26 2010-10-13 Airbus Uk Ltd Laminar flow monitor
CN102139766B (en) * 2010-10-28 2013-01-23 中国科学院力学研究所 Method for obtaining design attack angle for reentry flight of suborbital vehicle
US8831884B2 (en) * 2011-10-25 2014-09-09 Rosemount Aerospace Inc. Methods of determining water droplet size distributions of clouds
CN102890734A (en) * 2012-09-18 2013-01-23 天津空中代码工程应用软件开发有限公司 Method for building flight icing reduced-order model in flight simulator
US9612163B2 (en) * 2013-10-10 2017-04-04 The Boeing Company Methods and apparatus for detecting ice formation on aircraft
US9297755B2 (en) * 2013-10-24 2016-03-29 The Regents Of The University Of Michigan Ice and supercooled water detection system
CN104166348B (en) * 2014-08-07 2016-08-17 北京理工大学 A kind of dynamic stability decision method of single sliding block Moving mass control aircraft
CN105629725B (en) * 2014-10-31 2019-01-08 北京临近空间飞行器***工程研究所 A kind of elastic movement modeling method of rear rudder glide vehicle
EP3314120B1 (en) * 2015-06-24 2023-06-07 Vestas Wind Systems A/S Blade load sensing system for a wind turbine
CN105389438A (en) * 2015-11-25 2016-03-09 中国航空工业集团公司沈阳飞机设计研究所 Gravity-center normal overload calculation method for ski-jump takeoff of aircraft
GB2554062A (en) * 2016-08-22 2018-03-28 Norwegian Univ Of Science And Technology Icing control system
CN107127190B (en) * 2017-05-31 2023-04-07 中国空气动力研究与发展中心低速空气动力研究所 Active and passive anti-icing device and method with coupled super-hydrophobic material and vibration deicing
US20190039742A1 (en) * 2017-08-01 2019-02-07 Honeywell International Inc. Managing response to icing threat
CN107702879B (en) * 2017-09-20 2019-06-18 中国空气动力研究与发展中心计算空气动力研究所 A kind of aircraft dynamic ice ice type microstructure features prediction technique
EP3704022A1 (en) * 2017-10-30 2020-09-09 Battelle Memorial Institute Ice protection system and controller
CN107976296B (en) * 2017-11-13 2019-10-22 北京临近空间飞行器***工程研究所 A kind of aerodynamic characteristics of vehicle on-line identification method based on backtracking adaptive algorithm
FR3074145B1 (en) * 2017-11-28 2022-04-29 Zodiac Aerotechnics METHOD FOR DETECTING THE PRESENCE OF ICE CRYSTALS IN AN ENVIRONMENT
US10160550B1 (en) * 2017-12-21 2018-12-25 United Technologies Corporation Ice detection and real time ice estimation using innovative sensor data learning approaches
FR3077059B1 (en) * 2018-01-24 2020-01-31 Airbus Operations METHOD FOR DETECTING ICING CONDITIONS FOR AN AIRCRAFT BY SUPERVISED AUTOMATIC LEARNING
KR101982071B1 (en) * 2018-08-24 2019-05-24 주식회사 에그 (Egg) method for autonomous icing protection apparatus
CN109878740B (en) * 2019-03-25 2024-02-27 上海工程技术大学 Anti-icing and deicing device and method for wings and tail wings of turboprop aircraft

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5028929A (en) * 1990-04-30 1991-07-02 University Corporation For Atmospheric Research Icing hazard detection for aircraft
US6490510B1 (en) * 1999-04-30 2002-12-03 Thales Avionics S.A. Fixed multifunction probe for aircraft
US6757668B1 (en) * 1999-11-05 2004-06-29 General Electric Company Information fusion of classifiers in systems with partial redundant information
EP2829473A1 (en) * 2013-07-24 2015-01-28 Air China Limited Method and apparatus for monitoring a wing anti-icing valve
CN109911216A (en) * 2017-12-12 2019-06-21 波音公司 The system and method for hitting the position of the limit for correcting the water in aerofoil profile

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
大飞机机翼结冰对飞行动力学特性影响研究;邵元培等;《飞行力学》;20180228;第36卷(第1期);12-19 *

Also Published As

Publication number Publication date
CN112046761A (en) 2020-12-08

Similar Documents

Publication Publication Date Title
EP2434296B1 (en) Airspeed sensing system for an aircraft
US6273370B1 (en) Method and system for estimation and correction of angle-of-attack and sideslip angle from acceleration measurements
CN103323624B (en) Method and device for verifying the consistency of the measurements of the angle-of-attack probe of an aircraft
CN112046761B (en) Airplane icing on-line detection method based on statistical test and filtering
US20070130096A1 (en) Fault detection in artificial intelligence based air data systems
US7257470B2 (en) Fault isolation method and apparatus in artificial intelligence based air data systems
CN108170127B (en) A kind of fault detection method of UAV Flight Control System
US20190193866A1 (en) Method and assistance system for detecting a degradation of light performance
CN106840093A (en) A kind of detection method of unmanned plane during flying height, device and unmanned plane
CN109885854A (en) The real-time forecasting system of Flutter Boundaries and prediction technique based on arma modeling
US20060212181A1 (en) Method and apparatus for extending useful range of air data parameter calculation in flush air data systems
US20110264308A1 (en) Method And Device For Automatically Estimating An Air Speed Of An Aircraft
CN111290428A (en) Flight control system for an aircraft and method of controlling the same
CN116127842A (en) Post-fault flight envelope online prediction method based on radial basis-counter propagation neural network
JPH08211082A (en) Method and system for determining anemobaroclinometric parameter of flying aircraft
EP3819641A2 (en) Airflow angle, static and total pressure probe
EP0629166A1 (en) Aerodynamic pressure sensor systems
RU2364846C1 (en) Method for automated assessment of aircraft engines total thrust in process of flight
JP2022032051A5 (en)
Siu et al. Flight test results of an angle of attack and angle of sideslip calibration method using Output-Error optimization
CN112528492A (en) Fault detection method and device under wing damage condition
RU2214582C1 (en) Method establishing aerodynamic corrections to readings of pitot-static tubes
WO2020141316A1 (en) Improvements in or relating to angle of attack sensing
Korsun et al. Detection of dynamic errors in aircraft flight data
Majeed et al. Multi sensor data fusion based approach for the calibration of airdata systems

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