CN111257593A - Atmospheric data estimation and state monitoring method fusing navigation data - Google Patents
Atmospheric data estimation and state monitoring method fusing navigation data Download PDFInfo
- Publication number
- CN111257593A CN111257593A CN202010090161.XA CN202010090161A CN111257593A CN 111257593 A CN111257593 A CN 111257593A CN 202010090161 A CN202010090161 A CN 202010090161A CN 111257593 A CN111257593 A CN 111257593A
- Authority
- CN
- China
- Prior art keywords
- data
- atmospheric data
- aircraft
- fusion
- atmospheric
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P13/00—Indicating or recording presence, absence, or direction, of movement
- G01P13/02—Indicating direction only, e.g. by weather vane
- G01P13/025—Indicating direction only, e.g. by weather vane indicating air data, i.e. flight variables of an aircraft, e.g. angle of attack, side slip, shear, yaw
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
The invention discloses an atmospheric data estimation and state monitoring method fusing navigation data, which comprises the following steps: collecting atmospheric data and inertial navigation data, and analyzing the relation between the atmospheric data and the inertial navigation data; establishing a virtual atmospheric data system model, and calculating a fusion attack angle and a fusion sideslip angle; when the atmospheric data system works normally, the mean value of the difference between the fusion attack angle and the measurement attack angle is calculatedAnd the mean square error σ; and (3) performing different threshold judgment by adopting a variance v of a difference value between the fused attack angle and the measured attack angle within a period of time, and judging that the atmospheric data system is not abnormal when | v | < n σ. The method and the device can effectively monitor the data abnormity in time, and are beneficial to enhancing the safety and reliability of the aviation aircraft.
Description
Technical Field
The invention belongs to the field of aviation aircraft data monitoring, and particularly relates to an atmospheric data estimation and state monitoring method.
Background
The angle of attack, the angle of sideslip, and the vacuum speed are very important atmospheric parameters, and therefore other atmospheric information can be obtained. The atmospheric parameters represent the stress and heating conditions of the aircraft, are instruction information of a flight control system of the aerospace aircraft, and the measurement precision and reliability of the atmospheric parameters are directly related to the normal work and performance exertion of systems such as aircraft control and the like.
With the rapid development of the aviation industry, the security and reliability of data systems applied to aviation become more and more important. In order to ensure the flight safety of an aviation aircraft and the accurate measurement of an attack angle and a sideslip angle, an anomaly monitoring system needs to be configured for the attack angle and the sideslip angle in the aircraft.
Aviation aircraft are sometimes required to make extensive high-speed maneuvers or fly under severe conditions, and thus there is a possibility of failure of the air data system. For safety and reliability reasons, an aircraft must have atmospheric parameter anomaly monitoring functionality. In an aircraft, data abnormality is generally monitored by means of hardware redundancy, sensor redundancy, voting modules and the like, so that the accuracy of data information is ensured. However, considering the influence of factors such as equipment price, the anomaly monitoring of the atmospheric data needs to be realized in a set of INS/GNSS/ADS equipment. Therefore, three common anomaly monitoring modes cannot be used, so that a new method needs to be found to realize the function of monitoring the data anomaly of the general aviation aircraft.
Disclosure of Invention
In order to solve the technical problems mentioned in the background art, the invention provides an atmospheric data estimation and state monitoring method fusing navigation data.
In order to achieve the technical purpose, the technical scheme of the invention is as follows:
an atmospheric data estimation and state monitoring method fusing navigation data comprises the following steps:
(1) collecting atmospheric data and inertial navigation data, and analyzing the relation between the atmospheric data and the inertial navigation data;
(2) establishing a virtual atmospheric data system model, and calculating a fusion attack angle and a fusion sideslip angle;
(3) when the atmospheric data system works normally, the fusion attack angle and measurement are calculatedMean of difference between angles of attackAnd the mean square error σ; and (3) performing different threshold judgment by adopting a variance v of a difference value between the fused attack angle and the measured attack angle within a period of time delta T, judging that the atmospheric data system is not abnormal when | v | < n σ, and otherwise, judging that the atmospheric data system is abnormal, wherein n is a set coefficient.
Further, in step (1), the relationship between the atmospheric data and the inertial navigation data is as follows:
in the above formula, M is the flying Mach number of the aircraft; t issIs the Kelvin temperature; vTIs the speed of the aircraft relative to air, i.e. the vacuum speed;the speed of the air under the machine body system relative to the geographical system, namely the wind speed under the machine body system;a posture transfer matrix for the navigation system to the body system;the speed of the aircraft under the engine system relative to the geographic system, namely the ground speed of the engine system;the machine body is under vacuum speed;the vacuum velocity is the component of the machine system along the transverse axis of the machine body to the right,the component of the vacuum speed of the machine system along the longitudinal axis of the machine body is forward,the component of the vacuum speed along the vertical axis of the machine body under the machine system, α is a fusion attack angle, and β is a fusion sideslip angle.
Further, in step (2), based on the established virtual atmospheric data system model, considering that the aircraft is affected by the lever arm effect when making angular motion, the velocities above the nose, the right wing and the fuselage are different from the velocity of the aircraft centroid point, and the calculation formula of the fusion attack angle α and the fusion sideslip angle β is obtained:
in the above formula, the first and second carbon atoms are,linear velocity of three axes of mass center under the aircraft system, p, q, r are angular velocities of three axes under the aircraft system, xs,ys,zsThe lever arm length of the lower three shafts of the aircraft system.
Further:
in the above formula, theta, gamma, psi are attitude angles,is the first differential of theta, gamma and psi, i.e. the attitude angular velocity.
Further, in the step (3), selecting a plurality of time periods of aircraft flight, making a difference between the fusion attack angle and the measurement attack angle, and counting data to obtain a mean value and mean square error statistics; different abnormal threshold values are adopted in different time periods, and the maximum value and the minimum value in the difference values of the fusion attack angle and the measurement attack angle are respectively set as xmaxAnd xminLet ε equal max { | xmax|,|xminI.e. the coefficientWhereinIndicating a ceiling operation.
Further, in step (3), Δ T is 2 s.
Adopt the beneficial effect that above-mentioned technical scheme brought:
the virtual atmospheric data system fusing navigation data is adopted to estimate the atmospheric data and monitor the state of the atmospheric data, so that the data abnormality can be timely and effectively monitored, and the safety and the reliability of the aviation aircraft can be enhanced.
Drawings
FIG. 1 is a block diagram of a virtual atmosphere data system of the present invention;
FIG. 2 is a diagram of an anomaly monitoring algorithm of the present invention.
Detailed Description
The technical scheme of the invention is explained in detail in the following with the accompanying drawings.
The invention designs an atmospheric data estimation and state monitoring method fusing navigation data, which comprises the following steps:
step 1: collecting atmospheric data and inertial navigation data, and analyzing the relation between the atmospheric data and the inertial navigation data;
step 2: establishing a virtual atmospheric data system model, and calculating a fusion attack angle and a fusion sideslip angle;
and step 3: when the atmospheric data system works normally, the mean value of the difference between the fusion attack angle and the measurement attack angle is calculatedAnd the mean square error σ; and (3) performing different threshold judgment by adopting a variance v of a difference value between the fused attack angle and the measured attack angle within a period of time delta T, judging that the atmospheric data system is not abnormal when | v | < n σ, and otherwise, judging that the atmospheric data system is abnormal, wherein n is a set coefficient.
In this embodiment, in step 1, the acquired atmospheric data includes an attack angle, a sideslip angle, a mach number, a vacuum speed, and the like, and the acquired inertial navigation data includes an attitude angle, an attitude angle rate, a northeast direction speed, and the like.
The relationship between atmospheric data and inertial navigation data is as follows:
in the above formula, M is the flying Mach number of the aircraft; t issIs the Kelvin temperature; vTIs the speed of the aircraft relative to air, i.e. the vacuum speed;the speed of the air under the machine body system relative to the geographical system, namely the wind speed under the machine body system;a posture transfer matrix for the navigation system to the body system;the speed of the aircraft under the engine system relative to the geographic system, namely the ground speed of the engine system;the machine body is under vacuum speed;the vacuum velocity is the component of the machine system along the transverse axis of the machine body to the right,the component of the vacuum speed of the machine system along the longitudinal axis of the machine body is forward,the component of the vacuum speed along the vertical axis of the machine body under the machine system, α is a fusion attack angle, and β is a fusion sideslip angle.
In this embodiment, in step 2, the virtual atmosphere data system model is established as shown in fig. 1. As shown in FIG. 2, the wind speed is calculated from the vacuum speed and the ground speed at the previous time, and the fused vacuum speed at the current time is obtained from the ground speed at the current time and the wind speed at the previous time according to the principle that the wind speed is not changed for a short timeFurther obtain a fused angle of attackAnd blended sideslip angle
When the airplane does angular motion, due to the influence of the lever arm effect, the speed above the nose, the right wing and the airplane body is different from the speed of the center of mass point of the airplane, so that the influence of the lever arm effect needs to be eliminated by the following formula:
in the above formula, the first and second carbon atoms are,linear velocity of centroid (x, y, z) axis under aircraft system, angular velocity of (x, y, z) axis under aircraft system, xs,ys,zsIs the lever arm length of the lower (x, y, z) axis of the aircraft architecture.
The output of the inertial navigation system is the attitude angular velocity under the geographic systemThe attitude angular velocity under the geographic system needs to be converted into the angular velocity under the mechanical system (p, q, r) by the following formula.
On the basis of the above equation, considering the influence of the lever arm effect, the following equation is obtained:
through the formula, the model of the virtual atmospheric data system is established, and a more accurate fusion attack angle and a fusion sideslip angle can be obtained.
In this embodiment, in step 3, considering that a false alarm and a false alarm may be generated by using a difference at a certain time as a threshold value for determining, and affecting the accuracy of the atmospheric data anomaly monitoring method, the variance of the difference over a period of time is used for determining the anomaly threshold value. Since the system needs to monitor the fault within 2s after the fault occurs, the variance of the difference value within 2s is adopted for threshold judgment.
Selecting a plurality of time periods of aircraft flight, subtracting the fusion attack angle and the measurement attack angle, and carrying out statistics on data to obtain a mean value and mean square error statistics; different abnormal threshold values are adopted in different time periods, and the maximum value and the minimum value in the difference values of the fusion attack angle and the measurement attack angle are respectively set as xmaxAnd xminLet ε equal max { | xmax|,|xminI.e. the coefficientWhereinIndicating a ceiling operation.
The embodiments are only for illustrating the technical idea of the present invention, and the technical idea of the present invention is not limited thereto, and any modifications made on the basis of the technical scheme according to the technical idea of the present invention fall within the scope of the present invention.
Claims (6)
1. An atmospheric data estimation and state monitoring method fusing navigation data is characterized by comprising the following steps:
(1) collecting atmospheric data and inertial navigation data, and analyzing the relation between the atmospheric data and the inertial navigation data;
(2) establishing a virtual atmospheric data system model, and calculating a fusion attack angle and a fusion sideslip angle;
(3) when the atmospheric data system works normally, the mean value of the difference between the fusion attack angle and the measurement attack angle is calculatedAnd the mean square error σ; different threshold judgment is carried out by adopting variance v of the difference value between the fused attack angle and the measured attack angle within a period of time delta T, and when | v | < n σ, atmospheric data is judgedAnd (4) the system is not abnormal, otherwise, the atmospheric data system is abnormal, wherein n is a set coefficient.
2. The method for estimating atmospheric data and monitoring state of fused navigation data according to claim 1, wherein in step (1), the relationship between the atmospheric data and the inertial navigation data is as follows:
in the above formula, M is the flying Mach number of the aircraft; t issIs the Kelvin temperature; vTIs the speed of the aircraft relative to air, i.e. the vacuum speed;the speed of the air under the machine body system relative to the geographical system, namely the wind speed under the machine body system;a posture transfer matrix for the navigation system to the body system;for the speed of the aircraft relative to the geographical system under the system, i.e. the systemTying down the ground speed;the machine body is under vacuum speed;the vacuum velocity is the component of the machine system along the transverse axis of the machine body to the right,the component of the vacuum speed of the machine system along the longitudinal axis of the machine body is forward,the component of the vacuum speed along the vertical axis of the machine body under the machine system, α is a fusion attack angle, and β is a fusion sideslip angle.
3. The method for estimating atmospheric data and monitoring state according to the fused navigation data of claim 2, wherein in the step (2), based on the established virtual atmospheric data system model, the calculation formula of the fused attack angle α and the fused sideslip angle β is obtained by considering that the speed of the aircraft above the nose, the right wing and the fuselage is different from the speed of the aircraft center of mass point when the aircraft is subjected to the lever arm effect during the angular motion:
5. The method for estimating atmospheric data and monitoring state of fusion navigation data according to claim 1, wherein in the step (3), a plurality of time periods of aircraft flight are selected, the fusion angle of attack and the measured angle of attack are differentiated, and data are counted to obtain a mean value and a mean square error statistic; different abnormal threshold values are adopted in different time periods, and the maximum value and the minimum value in the difference values of the fusion attack angle and the measurement attack angle are respectively set as xmaxAnd xminLet ε equal max { | xmax|,|xminI.e. the coefficientWhereinIndicating a ceiling operation.
6. The method for estimating atmospheric data and monitoring the state of the merged navigation data according to claim 1, wherein Δ T is 2s in the step (3).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010090161.XA CN111257593B (en) | 2020-02-13 | 2020-02-13 | Atmospheric data estimation and state monitoring method fusing navigation data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010090161.XA CN111257593B (en) | 2020-02-13 | 2020-02-13 | Atmospheric data estimation and state monitoring method fusing navigation data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111257593A true CN111257593A (en) | 2020-06-09 |
CN111257593B CN111257593B (en) | 2021-05-28 |
Family
ID=70947279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010090161.XA Active CN111257593B (en) | 2020-02-13 | 2020-02-13 | Atmospheric data estimation and state monitoring method fusing navigation data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111257593B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114636842A (en) * | 2022-05-17 | 2022-06-17 | 成都信息工程大学 | Atmospheric data estimation method and device for hypersonic aircraft |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6253166B1 (en) * | 1998-10-05 | 2001-06-26 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Stable algorithm for estimating airdata from flush surface pressure measurements |
US6977608B1 (en) * | 2004-12-15 | 2005-12-20 | Rockwell Collins | Atmospheric data aggregation and forecasting system |
CN104568295A (en) * | 2014-12-08 | 2015-04-29 | 太原航空仪表有限公司 | Monitoring and processing method for pressure faults of embedded air data system |
CN105005099A (en) * | 2015-07-03 | 2015-10-28 | 南京航空航天大学 | Atmospheric parameter calculation method based on strapdown inertial navigation and flight control system |
CN105066994A (en) * | 2015-08-21 | 2015-11-18 | 中国运载火箭技术研究院 | Data fusion method for flush air data system and inertial navigation system |
CN106324643A (en) * | 2016-10-19 | 2017-01-11 | 山东科技大学 | Method for estimating airspeed of unmanned aerial vehicle and detecting failures of airspeed head |
CN106705996A (en) * | 2016-11-25 | 2017-05-24 | 北京航天自动控制研究所 | Aircraft navigation information correcting method based on atmospheric feature parameters |
CN107391891A (en) * | 2017-09-05 | 2017-11-24 | 北京理工大学 | A kind of high aspect ratio wing Optimization Design based on Model Fusion method |
EP3557227A1 (en) * | 2018-04-19 | 2019-10-23 | Honeywell International Inc. | System and method for deriving airspeed from a particle sensor |
CN110455310A (en) * | 2019-05-28 | 2019-11-15 | 中国空气动力研究与发展中心 | The atmosphere parameter measurement method of hypersonic aircraft |
CN110501759A (en) * | 2019-07-03 | 2019-11-26 | 中国商用飞机有限责任公司 | Processing method and processing system for multi-source atmospheric data parameters of airplane |
-
2020
- 2020-02-13 CN CN202010090161.XA patent/CN111257593B/en active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6253166B1 (en) * | 1998-10-05 | 2001-06-26 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Stable algorithm for estimating airdata from flush surface pressure measurements |
US6977608B1 (en) * | 2004-12-15 | 2005-12-20 | Rockwell Collins | Atmospheric data aggregation and forecasting system |
CN104568295A (en) * | 2014-12-08 | 2015-04-29 | 太原航空仪表有限公司 | Monitoring and processing method for pressure faults of embedded air data system |
CN105005099A (en) * | 2015-07-03 | 2015-10-28 | 南京航空航天大学 | Atmospheric parameter calculation method based on strapdown inertial navigation and flight control system |
CN105066994A (en) * | 2015-08-21 | 2015-11-18 | 中国运载火箭技术研究院 | Data fusion method for flush air data system and inertial navigation system |
CN106324643A (en) * | 2016-10-19 | 2017-01-11 | 山东科技大学 | Method for estimating airspeed of unmanned aerial vehicle and detecting failures of airspeed head |
CN106705996A (en) * | 2016-11-25 | 2017-05-24 | 北京航天自动控制研究所 | Aircraft navigation information correcting method based on atmospheric feature parameters |
CN107391891A (en) * | 2017-09-05 | 2017-11-24 | 北京理工大学 | A kind of high aspect ratio wing Optimization Design based on Model Fusion method |
EP3557227A1 (en) * | 2018-04-19 | 2019-10-23 | Honeywell International Inc. | System and method for deriving airspeed from a particle sensor |
CN110455310A (en) * | 2019-05-28 | 2019-11-15 | 中国空气动力研究与发展中心 | The atmosphere parameter measurement method of hypersonic aircraft |
CN110501759A (en) * | 2019-07-03 | 2019-11-26 | 中国商用飞机有限责任公司 | Processing method and processing system for multi-source atmospheric data parameters of airplane |
Non-Patent Citations (5)
Title |
---|
ZHANG,L等: "《Research on Adaptive Multi-Source Information Fault-Tolerant Navigation Method Based on No-Reference System Diagnosis》", 《SENSORS》 * |
叶玮等: "《FADS/INS组合法迎角、侧滑角测量方法研究》", 《飞机设计》 * |
李清东等: "《FADS快速智能故障检测和诊断技术》", 《系统工程与电子技术》 * |
杨胜江等: "《嵌入式大气数据传感与惯性导航信息融合方法研究》", 《战术导弹技术》 * |
陆辰: "《嵌入式大气数据系统算法及大气数据传感信息融合关键技术研究》", 《中国博士学位论文全文数据库》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114636842A (en) * | 2022-05-17 | 2022-06-17 | 成都信息工程大学 | Atmospheric data estimation method and device for hypersonic aircraft |
Also Published As
Publication number | Publication date |
---|---|
CN111257593B (en) | 2021-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111284716B (en) | Flight control system for determining faults based on errors between measured and estimated angles of attack | |
Guo et al. | Multisensor data-fusion-based approach to airspeed measurement fault detection for unmanned aerial vehicles | |
CN111433612A (en) | Neural network system for estimating combined training of aircraft aerial data based on model and flight information | |
EP3663772B1 (en) | Flight control system for determining estimated dynamic pressure based on lift and drag coefficients | |
US20070130096A1 (en) | Fault detection in artificial intelligence based air data systems | |
CN103884339B (en) | The equipment for configuring Navigation for Means of Transportation parameter value | |
EP3663774B1 (en) | Flight control system for determining a common mode pneumatic fault | |
Lyu et al. | A thrust model aided fault diagnosis method for the altitude estimation of a quadrotor | |
Lombaerts et al. | On-line safe flight envelope determination for impaired aircraft | |
CN115167508B (en) | Multi-rotor unmanned aerial vehicle sensor fault safety flight control system and method | |
CN108388229B (en) | Health degree-based four-rotor random hybrid system health assessment method | |
Whitmore et al. | In-flight demonstration of a real-time flush Airdata sensing (RT-FADS) system | |
CN111257593B (en) | Atmospheric data estimation and state monitoring method fusing navigation data | |
CN109612459A (en) | The fault-tolerant air navigation aid of quadrotor inertial sensor based on kinetic model | |
Whitmore et al. | In-flight demonstration of a real-time flush airdata sensing system | |
Deng et al. | An Improved Unscented Kalman Filter for Interrupted and Drift Sensor Faults of Aircrafts | |
Lerro et al. | Test in operative environment of an artificial neural network for aerodynamic angles estimation | |
CN116184986B (en) | Unmanned aerial vehicle fault detection method and system based on flight control log | |
CN112697340A (en) | Fixed-wing aircraft atmospheric data system and fault detection method thereof | |
CN116203931A (en) | Flight fault-tolerant control and reconstruction system under airspeed meter failure | |
Zahed et al. | Virtual sensor development for actual sensor fault detection and flight parameter estimation in real time | |
Raptis et al. | A particle filtering-based framework for real-time fault diagnosis of autonomous vehicles | |
Prabhu et al. | An innovative analytic redundancy approach to air data sensor fault detection | |
Majeed et al. | Multi sensor data fusion based approach for the calibration of airdata systems | |
CN112528492A (en) | Fault detection method and device under wing damage condition |
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 |