CN116380396B - Continuous wind tunnel aircraft atmospheric data system identification test system and method - Google Patents
Continuous wind tunnel aircraft atmospheric data system identification test system and method Download PDFInfo
- Publication number
- CN116380396B CN116380396B CN202310617749.XA CN202310617749A CN116380396B CN 116380396 B CN116380396 B CN 116380396B CN 202310617749 A CN202310617749 A CN 202310617749A CN 116380396 B CN116380396 B CN 116380396B
- Authority
- CN
- China
- Prior art keywords
- wind tunnel
- angle
- atmospheric
- model
- pressure
- 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
Links
- 238000013095 identification testing Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012360 testing method Methods 0.000 claims abstract description 82
- 230000003068 static effect Effects 0.000 claims abstract description 50
- 238000005259 measurement Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims description 36
- 238000009434 installation Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 5
- 238000012797 qualification Methods 0.000 claims description 4
- 238000010998 test method Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims 1
- 238000004891 communication Methods 0.000 abstract description 4
- 238000012827 research and development Methods 0.000 abstract description 4
- 238000012795 verification Methods 0.000 abstract description 2
- 230000008901 benefit Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009736 wetting Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/02—Wind tunnels
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The application provides a continuous wind tunnel aircraft atmospheric data system identification test system and method, and belongs to the technical field of atmospheric data system identification test methods. The system comprises an identification test computer, a wind tunnel measurement and control system, an atmosphere data computer, an atmosphere pressure sensor, a wind tunnel total pressure sensor, a wind tunnel static pressure sensor and an aircraft model; the aircraft model is arranged in the wind tunnel test section and is connected with the atmospheric pressure sensor; the atmospheric pressure sensor is connected with the atmospheric data computer; the atmosphere data computer is connected with the identification computer; the identification computer is connected with the wind tunnel measurement and control system through a communication cable; the wind tunnel measurement and control system is respectively connected with the wind tunnel total pressure sensor and the wind tunnel static pressure sensor through cables; the method solves the technical problems of high dependence degree and high cost on flight tests in the traditional atmospheric data system test in the prior art, long system research and development period and high stability and accuracy verification difficulty on the atmospheric data system in the dangerous envelope range.
Description
Technical Field
The application relates to a system identification test method, in particular to a continuous wind tunnel aircraft atmospheric data system identification test system and method, and belongs to the technical field of atmospheric data system identification test methods.
Background
The air data system is a comprehensive and high-precision air data information system for the aircraft, and the flight control system, the navigation system, the instrument display system and the like of the modern aircraft all need accurate static pressure, dynamic pressure, altitude, indicated airspeed and other information. The atmospheric data system is an important on-board electronic device, and the accuracy of the calculated atmospheric parameters is directly related to the flight safety of the aircraft. The precedent that the output of error signals causes dangerous situations or even accidents of the aircraft is frequent due to the abnormal operation of the atmospheric data system (sensor wetting, icing, program logic problems, etc.).
The existing air data system test mainly uses an air data system component as an independent test piece or uses the appearance of a part of aircraft as a test piece to carry out a wind tunnel test, the test method is used as a conventional pressure measurement wind tunnel test to obtain pressure data in different airspeed, attitude angles and other states, after the test is completed, the air data system correction parameters are obtained through data post-processing, and finally, a flight test is carried out to carry out checking and correction to obtain final air data system parameters.
The existing research method uses a flight test as a checking means to check the stability, the accuracy and the like of an atmospheric data system, and has the defects of higher cost of the flight test, additional research and development period and cost are added after problems are found in the flight test, and the flight test has higher danger to the dangerous envelope checking of the atmospheric data system, influences the safety of an aircraft and is generally difficult to find. With the development of the aircraft technology, the requirements on the accuracy and stability of the air data system are gradually increased, and the requirements on the maneuverability and stall state flight of the aircraft are increased, so that the requirements on the dangerous envelope of the air data system are also increasingly important.
Disclosure of Invention
The following presents a simplified summary of the application in order to provide a basic understanding of some aspects of the application. It should be understood that this summary is not an exhaustive overview of the application. It is not intended to identify key or critical elements of the application or to delineate the scope of the application. Its purpose is to present some concepts in a simplified form as a prelude to the more detailed description that is discussed later.
In view of the above, the application provides a continuous wind tunnel aircraft atmospheric data system identification test system and a continuous wind tunnel aircraft atmospheric data system identification test method, which are used for solving the technical problems of high dependence on flight tests, high cost, long system research and development period and high stability and accuracy verification difficulty of a dangerous envelope range atmospheric data system in the traditional atmospheric data system test in the prior art.
Scheme one: a continuous wind tunnel aircraft atmospheric data system identification test system comprises an identification test computer, a wind tunnel measurement and control system, an atmospheric data computer, an atmospheric pressure sensor, a wind tunnel total pressure sensor, a wind tunnel static pressure sensor and an aircraft model;
the aircraft model is arranged in the wind tunnel test section and is connected with an atmospheric pressure sensor through a gas path hose;
the atmospheric pressure sensor is connected with an atmospheric data computer;
the atmosphere data computer is connected with the identification computer;
the identification computer is connected with the wind tunnel measurement and control system;
the wind tunnel measurement and control system is respectively connected with the wind tunnel total pressure sensor and the wind tunnel static pressure sensor.
Scheme II: a continuous wind tunnel aircraft atmospheric data system identification test method comprises the following steps:
s1, controlling the total pressure of a wind tunnel according to the requirements of an identification test, and controlling Mach number of a test section by adjusting the rotating speed of a wind tunnel compressor to establish the required flow field speed, total pressure and static pressure conditions in an aircraft model installation area of the wind tunnel test section;
s2, transmitting the pressure of the position of the airspeed tube or the pressure measuring hole on the aircraft model to an atmospheric pressure sensor through an air path hose;
s3, transmitting the pressure value acquired in real time to an atmosphere data computer by the atmosphere pressure sensor;
s4, the atmospheric data computer calculates the total pressure, the static pressure, the air pressure height, the indicated airspeed, the attack angle, the sideslip angle and the Mach number parameters of the position of the aircraft model in real time according to the received pressure value, and transmits the parameters to the identification test computer;
s5, the identification test computer carries out an identification test on the air data system.
Preferably, S5 specifically includes the following steps:
s51, calibrating a pressure sensor;
s52, correcting the pressure delay of the air pipeline;
s53, calculating atmospheric parameters of a wind tunnel flow field;
s54, atmosphere parameter identification and comparison.
Preferably, the calibration of the S51 pressure sensor, specifically, taking the height of the position of the test model in the wind tunnel test process as a reference, and maintaining the actual height difference between the total pressure and static pressure sensors for controlling the wind tunnel flow field and the position of the model during the calibrationPerforming calibration, and maintaining the actual height difference of the atmospheric pressure sensor and the position of the model during the calibration>Calibration is carried out to eliminate the difference in height>And->The resulting pressure deviation.
Preferably, the pressure delay correction of the air pipeline in the S52 is specifically to test the pressure delay time generated by the length of the air pipeline for test under different static pressure conditions when the different static pressures are obtained.
Preferably, the calculating of the atmospheric parameters of the S53 wind tunnel flow field specifically comprises the following steps:
s531, directly measuring to obtain total pressure of wind tunnel flow fieldTotal temperature of wind tunnel flow field->Static pressure of wind tunnel residence chamber>The pitch angle theta, the roll angle phi, the pitch installation angle delta theta, the yaw installation angle delta phi, the roll rotation angleA mounting angle delta phi;
s532, resolving to obtain parameters of the Mach number, the total pressure, the static pressure, the air pressure height, the indicated airspeed, the attack angle and the sideslip angle of the model of the wind tunnel at the position of the aircraft model.
Preferably, the wind tunnel flow field Mach number, total pressure, static pressure, air pressure height, indicated airspeed, attack angle and sideslip angle parameters of the model are obtained by calculation:
mach number of wind tunnel flow field
Wherein M is the Mach number of the wind tunnel flow field,for the wind tunnel residence chamber reference point mach number calculated from the residence chamber static pressure,is->And->The correction amount of the test section is a specific calibration relation of each test section of each wind tunnel;
static pressure of wind tunnel flow field is
The air pressure is high
When γ+.0:
when γ=0:
wherein,,is the standard gravity acceleration>The method comprises the steps of carrying out a first treatment on the surface of the R is the gas constant of air, r= 287.05287 (++>/K/>)；/>A lower limit value of the air pressure height of the corresponding height layer; />The lower limit value of the atmospheric static pressure of the corresponding height layer; />An atmospheric temperature lower limit value for the corresponding height layer; gamma is the vertical temperature of the corresponding height layer;
indicating airspeed as
Wherein,,representing sea level standard sound velocity，/>；
The attack angle of the model isThe sideslip angle is +.>
Correcting the installation angle from the angle sensor shafting to the model body shafting meets the following matrix relation:
wherein,,euler rotation matrix of angle sensor relative to ground axis>For the rotation matrix of the angle sensor relative to the model axis, < >>For the rotation matrix of the model axis relative to the ground axis, < >>A rotation matrix of the model shafting relative to the angle sensor;
the rotation matrix from the ground axis to the model body axis is as follows:
wherein,,for Euler matrix rotated around x-axis, < >>For Euler matrix rotated around y-axis, < >>For Euler matrix rotated around z-axis, < >>For the roll angle from the ground axis to the test model body axis, +.>For yaw angle from the ground axis to the test model body axis, +.>The pitch angle is from the ground axis system to the test model body axis system;
the mounting angle between the angle sensor and the model reference plane satisfies the following relation:
wherein the angle sensor is arranged at a pitching installation angle with the modelYaw mounting angle->Roll mounting angle->The measuring method of (2) is that after 0 DEG leveling of the model reference plane, the angle sensor measures +.>And->The model reference plane is 90 DEG vertical to the ground, and the angle sensor measures +.>；
The rotation matrix from the ground axis to the angle sensor is as follows:
wherein,,for the roll angle of the angle sensor output, +.>Yaw angle output by the angle sensor, < +.>，/>Pitch angle output for angle sensor>
And->The matrix multiplication results in the following matrix:
from the following componentsThe method comprises the following steps:
angle of attack and sideslip angle of the model:
。
preferably, S54. The method comprises performing time stamp alignment on an atmospheric data computer, a wind tunnel measurement and control system and an identification test computer before starting a test, wherein the identification test computer receives parameters transmitted by the wind tunnel measurement and control system and parameters transmitted by the atmospheric data computer in real time during the test, performing delay time delta t correction on the received atmospheric data parameters transmitted by the atmospheric data computer, and comparing the corrected atmospheric data parameters with the atmospheric parameters calculated by the parameters transmitted by the wind tunnel measurement and control system, wherein the difference of the atmospheric data parameters is equal to or greater than the difference of the atmospheric data parametersWherein->Is->Respectively representing the calculated values of the same air data parameter at the same moment output by the wind tunnel flow field acquisition system and the air data computer, obtaining the difference of each air data parameter, and judging whether the difference meets the design index.
The third scheme is an electronic device, comprising a memory and a processor, wherein the memory stores a computer program, and the processor realizes the step of the continuous wind tunnel aircraft atmospheric data system identification test method in the second scheme when executing the computer program.
A fourth aspect is a computer readable storage medium having a computer program stored thereon, wherein the computer program when executed by a processor implements a continuous wind tunnel aircraft atmospheric data system identification test method as described in the second aspect.
The beneficial effects of the application are as follows: the application combines the advantages of variable pressure of a continuous transonic wind tunnel, simulates the atmospheric pressure of different heights, airspeeds and the like of the flying vehicle, and obtains Mach number, airspeeds and barometric altitude of the flying vehicle through measurement and calculation of the total pressure of the wind tunnel, a static pressure sensor and a total temperature sensor; by combining the advantage of higher wind tunnel test safety, the dangerous angle envelope of the aircraft is simulated, the attack angle and sideslip angle of the aircraft are calculated through an angle sensor arranged in the model, and the complete envelope range of the atmosphere data system is checked. According to the application, the part identified by the flight test of the aircraft is verified in the ground wind tunnel test, so that problems can be found in the ground test stage and correction can be completed in the research and development stage, the flight test cost is reduced, the flight test safety is improved, and the identification test speed and angle boundary are enlarged.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute a limitation on the application. In the drawings:
FIG. 1 is a schematic diagram of a continuous wind tunnel aircraft atmospheric data system qualification test system;
FIG. 2 is a schematic diagram of a continuous wind tunnel aircraft atmospheric data system identification test method;
FIG. 3 is a highly schematic view of sensor calibration according to the present application;
FIG. 4 is a schematic diagram of pressure delay;
fig. 5 is a schematic diagram of the law of pressure delay with static pressure.
Detailed Description
In order to make the technical solutions and advantages of the embodiments of the present application more apparent, the following detailed description of exemplary embodiments of the present application is provided in conjunction with the accompanying drawings, and it is apparent that the described embodiments are only some embodiments of the present application and not exhaustive of all embodiments. It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other.
Example 1, the description of this embodiment with reference to fig. 1, is a continuous wind tunnel aircraft atmospheric data system identification test system, which is characterized by comprising an identification test computer, a wind tunnel measurement and control system, an atmospheric data computer, an atmospheric pressure sensor, a wind tunnel total pressure sensor, a wind tunnel static pressure sensor and an aircraft model (respectively represented by test models in fig. 1 and 3);
the aircraft model is arranged in the wind tunnel test section and is connected with an atmospheric pressure sensor through a gas path hose;
specifically, the aircraft model comprises an aircraft appearance physical model, an atmospheric data system measuring piece such as a airspeed tube and the like arranged on the aircraft model, a pressure measuring hole for calculating the atmospheric data system applied to the surface of the aircraft model, and an angle sensor arranged in the model;
the atmospheric pressure sensor is connected with the atmospheric data computer through a cable;
the atmosphere data computer is connected with the identification computer through a communication cable;
the identification computer is connected with the wind tunnel measurement and control system through a communication cable;
the wind tunnel measurement and control system is respectively connected with the wind tunnel total pressure sensor and the wind tunnel static pressure sensor through cables;
specifically, the atmospheric pressure sensor and the atmospheric data computer are arranged outside the wind tunnel test section and in the wind tunnel residence chamber.
The operating principle of the continuous wind tunnel aircraft atmospheric data system identification test system is as follows: according to the requirements of the identification test, the required conditions such as flow field speed, total pressure, static pressure and the like are established in the aircraft model installation area of the wind tunnel test section, the airspeed tube or the pressure measuring holes on the surface of the aircraft sense the pressure value of the position in real time, the pressure value is transmitted to the atmospheric pressure sensor in the wind tunnel residence chamber outside the test section through the air path hose, the atmospheric pressure sensor measures the pressure value of each pressure measuring hole in real time and transmits the pressure value to the atmospheric data computer, the atmospheric data computer calculates the airspeed, the pressure and other atmospheric parameters of the position of the aircraft in real time and transmits the atmospheric parameters to the identification test computer through the communication cable, and the identification test computer synchronously, stores and compares the atmospheric data parameters transmitted by the wind tunnel measurement and control system and the atmospheric data computer, so that the identification test of the atmospheric data system is completed.
Example 2, the present embodiment will be described with reference to fig. 2 to 5, which is a continuous wind tunnel aircraft atmospheric data system identification test method, wherein an atmospheric data system and an aircraft model are installed in a wind tunnel, and the wind tunnel flow field is used as a reference to perform an assessment and identification on the atmospheric data system, and the stability and the atmospheric data calculation accuracy of the assessment and identification are examined; the advantages of variable pressure of the continuous transonic wind tunnel are combined, the atmospheric pressure of different heights, airspeeds and the like of the flying of the aircraft can be simulated in the wind tunnel, mach numbers, airspeeds, barometric pressures and the like of the aircraft are obtained through calculation of the pressures measured by the wind tunnel total pressure and the static pressure sensor, meanwhile, the dangerous angle envelope of the aircraft can be simulated by combining the advantages of higher safety of the wind tunnel test, the attack angle and sideslip angle of the aircraft are obtained through calculation by an angle sensor arranged in the model, and the complete envelope range of the atmospheric data system is checked; the method specifically comprises the following steps:
s1, controlling the total pressure of a wind tunnel according to the requirements of an identification test, and controlling Mach number of a test section by adjusting the rotating speed of a wind tunnel compressor to establish the required flow field speed, total pressure and static pressure conditions in an aircraft model installation area of the wind tunnel test section;
s2, transmitting the pressure of the position of the airspeed tube or the pressure measuring hole on the aircraft model to an atmospheric pressure sensor through an air path hose;
s3, transmitting the pressure value acquired in real time to an atmosphere data computer by the atmosphere pressure sensor;
s4, the atmospheric data computer calculates the total pressure, the static pressure, the air pressure height, the indicated airspeed, the attack angle, the sideslip angle and the Mach number parameters of the position of the aircraft model in real time according to the received pressure value, and transmits the parameters to the identification test computer;
s5, an identification test computer carries out an identification test on the air data system, and the method specifically comprises the following steps:
s51, calibrating the pressure sensor, specifically, referring to FIG. 3, by taking the height of the position of the test model in the wind tunnel test process as a reference, performing the following stepsThe actual height difference between the total pressure and the position of the model is maintained during the calibration of the static pressure sensor for controlling the wind tunnel flow fieldPerforming calibration, and maintaining the actual height difference of the atmospheric pressure sensor and the position of the model during the calibration>Calibration is carried out to eliminate the difference in height>And->The resulting pressure deviation was about 12Pa/m under normal pressure conditions.
S52, correcting pressure delay of the air pipeline, specifically referring to fig. 4-5, because the volume of the atmospheric pressure sensor is larger, the space of the inner cavity of the test model is limited, the atmospheric pressure sensor cannot be installed, the atmospheric pressure sensor is installed outside a wind tunnel test section and in a residence chamber, a pressure transmission delay delta t is generated by an air path hose between the atmospheric pressure sensor and a pressure measuring hole of the wind tunnel test model, and the pressure delay time generated by the length of the air pipeline for test under different static pressure conditions is required to be tested before the test.
S53, calculating atmospheric parameters of a wind tunnel flow field, and specifically comprising the following steps of:
s531, directly measuring to obtain total pressure of wind tunnel flow fieldTotal temperature of wind tunnel flow field->Static pressure of wind tunnel residence chamber>The pitch angle theta of the model, the roll angle phi of the model, the pitch installation angle delta theta of the angle sensor and the model, the yaw installation angle delta phi of the angle sensor and the model, and the roll installation of the angle sensor and the modelAngle delta phi;
s532, resolving to obtain parameters of the Mach number, the total pressure, the static pressure, the air pressure height, the indicated airspeed, the attack angle and the sideslip angle of the model of the wind tunnel at the position of the aircraft model.
Mach number of wind tunnel flow field
Wherein M is the Mach number of the wind tunnel flow field,for the wind tunnel residence chamber reference point mach number calculated from the residence chamber static pressure,is->And->The correction amount of the test section is a specific calibration relation of each test section of each wind tunnel;
static pressure of wind tunnel flow field is
The air pressure is high
When γ+.0:
when γ=0:
wherein,,is the standard gravity acceleration>The method comprises the steps of carrying out a first treatment on the surface of the R is the gas constant of air, r= 287.05287 (++>/K/>)；/>A lower limit value of the air pressure height of the corresponding height layer; />The lower limit value of the atmospheric static pressure of the corresponding height layer; />An atmospheric temperature lower limit value for the corresponding height layer; gamma is the vertical temperature of the corresponding height layer; details are shown in Table 1->、A relationship table between the two;
TABLE 1、/>Relationship table between
Indicating airspeed as
Wherein,,represents sea level standard sound velocity,/->。
The attack angle of the model isThe sideslip angle is +.>
Correcting the installation angle from the angle sensor shafting to the model body shafting meets the following matrix relation:
wherein,,euler rotation matrix of angle sensor relative to ground axis>For the rotation matrix of the angle sensor relative to the model axis, < >>For the rotation matrix of the model axis relative to the ground axis, < >>A rotation matrix of the model shafting relative to the angle sensor;
the rotation matrix from the ground axis to the model body axis is as follows:
wherein,,for Euler matrix rotated around x-axis, < >>For Euler matrix rotated around y-axis, < >>For Euler matrix rotated around z-axis, < >>For the roll angle from the ground axis to the test model body axis, +.>For yaw angle from the ground axis to the test model body axis, +.>The pitch angle is from the ground axis system to the test model body axis system;
the mounting angle between the angle sensor and the model reference plane satisfies the following relation:
wherein the angle sensor is arranged at a pitching installation angle with the modelYaw mounting angle->Roll mounting angle->The measuring method of (2) is that after 0 DEG leveling of the model reference plane, the angle sensor measures +.>And->The model reference plane is 90 DEG vertical to the ground, and the angle sensor measures +.>；
The rotation matrix from the ground axis to the angle sensor is as follows:
wherein,,for the roll angle of the angle sensor output, +.>Yaw angle output by the angle sensor, < +.>，/>A pitch angle output by the angle sensor;
and->The matrix multiplication results in the following matrix:
from the following componentsThe method comprises the following steps:
angle of attack and sideslip angle of the model:
。
s54, identifying and comparing the atmospheric parameters, namely, carrying out timestamp alignment on an atmospheric data computer, a wind tunnel measurement and control system and an identification test computer before the test starts, wherein the identification test computer receives the parameters sent by the wind tunnel measurement and control system and the parameters sent by the atmospheric data computer in real time in the test process, carrying out delay time delta t correction on the received atmospheric data parameters sent by the atmospheric data computer, and then comparing the atmospheric data parameters with the atmospheric parameters calculated by the parameters sent by the wind tunnel measurement and control system, wherein the difference of the atmospheric data parameters is calculatedWherein->Is->Respectively representing the calculated values of the same air data parameter at the same moment output by the wind tunnel flow field acquisition system and the air data computer, obtaining the difference of each air data parameter, and judging whether the difference meets the design index.
Different system design indexes and requirements are different, are not important improvement points of the application, and generally require Mach numbers not exceeding +/-0.005, total pressures not exceeding +/-100 Pa, static pressures not exceeding +/-80 Pa, air pressure heights not exceeding +/-20 m, indicated airspeed not exceeding +/-6 km/h, attack angles and sideslip angles not exceeding +/-0.8 degrees.
The application can complete the operation of the identification and correction of the air data system in the flight test in the ground wind tunnel test, thereby reducing the test cost;
according to the application, dangerous flight attitudes which cannot be realized in a flight test can be realized, and more comprehensive atmospheric data system parameters are obtained.
In embodiment 3, the computer device of the present application may be a device including a processor and a memory, for example, a single chip microcomputer including a central processing unit. And the processor is used for executing the computer program stored in the memory to realize the steps of the continuous wind tunnel aircraft atmospheric data system identification test method.
The processor may be a central processing unit (Central Processing Unit, CPU), other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), off-the-shelf programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program (such as a sound playing function, an image playing function, etc.) required for at least one function, and the like; the storage data area may store data (such as audio data, phonebook, etc.) created according to the use of the handset, etc. In addition, the memory may include high-speed random access memory, and may also include non-volatile memory, such as a hard disk, memory, plug-in hard disk, smart Media Card (SMC), secure Digital (SD) Card, flash Card (Flash Card), at least one disk storage device, flash memory device, or other volatile solid-state storage device.
Embodiment 4, computer-readable storage Medium embodiment
The computer readable storage medium of the present application may be any form of storage medium readable by a processor of a computer device, including but not limited to, nonvolatile memory, volatile memory, ferroelectric memory, etc., having a computer program stored thereon, which when read and executed by the processor of the computer device, implements the steps of a continuous wind tunnel aircraft air data system qualification test method as described above.
The computer program comprises computer program code which may be in source code form, object code form, executable file or in some intermediate form, etc. The computer readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth. It should be noted that the computer readable medium contains content that can be appropriately scaled according to the requirements of jurisdictions in which such content is subject to legislation and patent practice, such as in certain jurisdictions in which such content is subject to legislation and patent practice, the computer readable medium does not include electrical carrier signals and telecommunication signals.
While the application has been described with respect to a limited number of embodiments, those skilled in the art, having benefit of the above description, will appreciate that other embodiments are contemplated within the scope of the application as described herein. Furthermore, it should be noted that the language used in the specification has been principally selected for readability and instructional purposes, and may not have been selected to delineate or circumscribe the inventive subject matter. Accordingly, many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the appended claims. The disclosure of the present application is intended to be illustrative, but not limiting, of the scope of the application, which is defined by the appended claims.
Claims (5)
1. The continuous wind tunnel aircraft atmospheric data system identification test system is characterized by comprising an identification test computer, a wind tunnel measurement and control system, an atmospheric data computer, an atmospheric pressure sensor, a wind tunnel total pressure sensor, a wind tunnel static pressure sensor and an aircraft model;
the aircraft model is arranged in the wind tunnel test section and is connected with an atmospheric pressure sensor through a gas path hose;
the atmospheric pressure sensor is connected with an atmospheric data computer;
the atmosphere data computer is connected with the identification computer;
the identification computer is connected with the wind tunnel measurement and control system, and the identification test computer carries out an identification test on the air data system, and specifically comprises the following steps:
s51, calibrating a pressure sensor, namely, taking the height of the position of a test model in the wind tunnel test process as a reference, and maintaining the actual height difference delta H between the total pressure and static pressure sensors for wind tunnel flow field control and the position of the model during calibration 2 Performing calibration, and maintaining the actual height difference delta H between the atmospheric pressure sensor and the position of the model during the calibration of the atmospheric pressure sensor 1 Calibration is performed to eliminate the height difference delta H 1 And DeltaH 2 The resulting pressure deviation;
s52, correcting the pressure delay of the air pipeline, namely testing the pressure delay time generated by the length of the air pipeline for test under different static pressures when different static pressures are obtained;
s53, calculating atmospheric parameters of a wind tunnel flow field, wherein the method comprises the following steps of:
s531, directly measuring to obtain total pressure P of wind tunnel flow field 0 Total temperature T of wind tunnel flow field 0 Static pressure P of wind tunnel residence chamber ct Model pitch angle θ, model roll angleAngle sensor and modelPitch mounting angle Δθ, yaw mounting angle Δψ of angle sensor and model, roll mounting angle +.>
S532, resolving to obtain parameters of the Mach number, total pressure, static pressure, air pressure height, indicated airspeed, attack angle and sideslip angle of the wind tunnel flow field where the aircraft model is located;
s54, identifying and comparing the atmospheric parameters, namely, carrying out timestamp alignment on an atmospheric data computer, a wind tunnel measurement and control system and an identification test computer before the test starts, wherein the identification test computer receives the parameters sent by the wind tunnel measurement and control system and the parameters sent by the atmospheric data computer in real time in the test process, carrying out delay time delta t correction on the received atmospheric data parameters sent by the atmospheric data computer, and then comparing the atmospheric data parameters with the atmospheric parameters calculated by the parameters sent by the wind tunnel measurement and control system, wherein the difference of the atmospheric data parameters is calculatedWherein u is i S i Respectively representing the calculated values of the same air data parameter at the same moment output by the wind tunnel flow field acquisition system and the air data computer, obtaining the difference of each air data parameter, and judging whether the difference meets the design index;
the wind tunnel measurement and control system is respectively connected with the wind tunnel total pressure sensor and the wind tunnel static pressure sensor.
2. The continuous wind tunnel aircraft atmospheric data system identification test method is characterized by comprising the following steps of:
s1, controlling the total pressure of a wind tunnel according to the requirements of an identification test, and controlling Mach number of a test section by adjusting the rotating speed of a wind tunnel compressor to establish the required flow field speed, total pressure and static pressure conditions in an aircraft model installation area of the wind tunnel test section;
s2, transmitting the pressure of the position of the airspeed tube or the pressure measuring hole on the aircraft model to an atmospheric pressure sensor through an air path hose;
s3, transmitting the pressure value acquired in real time to an atmosphere data computer by the atmosphere pressure sensor;
s4, the atmospheric data computer calculates the total pressure, the static pressure, the air pressure height, the indicated airspeed, the attack angle, the sideslip angle and the Mach number parameters of the position of the aircraft model in real time according to the received pressure value, and transmits the parameters to the identification test computer;
s5, an identification test computer carries out an identification test on the air data system, and the method specifically comprises the following steps:
s51, calibrating a pressure sensor, namely, taking the height of the position of a test model in the wind tunnel test process as a reference, and maintaining the actual height difference delta H between the total pressure and static pressure sensors for wind tunnel flow field control and the position of the model during calibration 2 Performing calibration, and maintaining the actual height difference delta H between the atmospheric pressure sensor and the position of the model during the calibration of the atmospheric pressure sensor 1 Calibration is performed to eliminate the height difference delta H 1 And DeltaH 2 The resulting pressure deviation;
s52, correcting the pressure delay of the air pipeline, namely testing the pressure delay time generated by the length of the air pipeline for test under different static pressures when different static pressures are obtained;
s53, calculating atmospheric parameters of a wind tunnel flow field, wherein the method comprises the following steps of:
s531, directly measuring to obtain total pressure P of wind tunnel flow field 0 Total temperature T of wind tunnel flow field 0 Static pressure P of wind tunnel residence chamber ct Model pitch angle θ, model roll anglePitch mounting angle Δθ of angle sensor and model, yaw mounting angle Δψ of angle sensor and model, roll mounting angle +.>
S532, resolving to obtain parameters of the Mach number, total pressure, static pressure, air pressure height, indicated airspeed, attack angle and sideslip angle of the wind tunnel flow field where the aircraft model is located;
s54, identifying and comparing the atmospheric parameters, namely, carrying out timestamp alignment on an atmospheric data computer, a wind tunnel measurement and control system and an identification test computer before the test starts, wherein the identification test computer receives the parameters sent by the wind tunnel measurement and control system and the parameters sent by the atmospheric data computer in real time in the test process, carrying out delay time delta t correction on the received atmospheric data parameters sent by the atmospheric data computer, and then comparing the atmospheric data parameters with the atmospheric parameters calculated by the parameters sent by the wind tunnel measurement and control system, wherein the difference of the atmospheric data parameters is calculatedWherein u is i S i Respectively representing the calculated values of the same air data parameter at the same moment output by the wind tunnel flow field acquisition system and the air data computer, obtaining the difference of each air data parameter, and judging whether the difference meets the design index.
3. The continuous wind tunnel aircraft atmospheric data system identification test method according to claim 2, wherein wind tunnel flow field Mach number, total pressure, static pressure, air pressure height, indicated airspeed, attack angle and sideslip angle parameters of the model are obtained by calculation:
mach number of wind tunnel flow field
M＝M CT +ΔM
Wherein M is the Mach number of the wind tunnel flow field, M CT For wind tunnel residence chamber reference point Mach number calculated by residence chamber static pressure, deltaM is M CT Correction amount of M is a specific calibration relation of each test section of each wind tunnel;
the static pressure of the wind tunnel flow field is P cp
The pressure height is H p
When γ+.0:
when γ=0:
wherein g n G is the standard gravity acceleration n ＝9.80665m/s 2 ) The method comprises the steps of carrying out a first treatment on the surface of the R is the gas constant of air, r= 287.05287 (m 2 /K*s 2 )；H b A lower limit value of the air pressure height of the corresponding height layer; p (P) b The lower limit value of the atmospheric static pressure of the corresponding height layer; t (T) b An atmospheric temperature lower limit value for the corresponding height layer; gamma is the vertical temperature of the corresponding height layer;
indicating airspeed of V
Wherein C is n Represents the standard sound velocity of sea level, C n ＝1225.0584km/h；
The attack angle of the model is alpha, and the sideslip angle is beta
Correcting the installation angle from the angle sensor shafting to the model body shafting meets the following matrix relation:
R ig ＝R im ·R mg
wherein R is ig R is Euler rotation matrix of angle sensor relative to ground shaft system im R is a rotation matrix of the angle sensor relative to a model shaft system mg R is a rotation matrix of a model shafting relative to a ground shafting mi A rotation matrix of the model shafting relative to the angle sensor;
the rotation matrix from the ground axis to the model body axis is as follows:
wherein R is x R is an Euler matrix rotated about the x-axis y R is an Euler matrix rotated about the y-axis z For Euler matrix rotated about the z-axis, phi mg To roll angle from the earth axis to the test model body axis, ψ mg For yaw angle, θ, from ground axis to test model body axis mg The pitch angle is from the ground axis system to the test model body axis system;
the mounting angle between the angle sensor and the model reference plane satisfies the following relation:
the measuring method of the angle sensor and the pitching installation angle delta theta, the yaw installation angle delta phi and the rolling installation angle delta phi of the model comprises the steps that after the model reference plane is leveled by 0 degrees, the angle sensor measures delta theta and delta phi, the model reference plane is vertical to the ground by 90 degrees, and the angle sensor measures delta phi;
the rotation matrix from the ground axis to the angle sensor is as follows:
wherein phi is i Roll angle, ψ, for the output of the angle sensor i For the yaw angle, ψ, output by the angle sensor i ＝0，θ i A pitch angle output by the angle sensor;
R mi and R is ig The matrix multiplication results in the following matrix:
from R mg ＝R mi ·R ig The method comprises the following steps:
angle of attack and sideslip angle of the model:
β＝asin(cosφ m sinψ m cosθ m +sinφ m sinθ m )。
4. an electronic device comprising a memory and a processor, the memory storing a computer program, the processor implementing the steps of a continuous wind tunnel aircraft atmospheric data system qualification test method of claim 2 or 3 when executing the computer program.
5. A computer readable storage medium having stored thereon a computer program, which when executed by a processor implements a continuous wind tunnel aircraft atmospheric data system qualification test method according to claim 2 or 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310617749.XA CN116380396B (en) | 2023-05-30 | 2023-05-30 | Continuous wind tunnel aircraft atmospheric data system identification test system and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310617749.XA CN116380396B (en) | 2023-05-30 | 2023-05-30 | Continuous wind tunnel aircraft atmospheric data system identification test system and method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116380396A CN116380396A (en) | 2023-07-04 |
| CN116380396B true CN116380396B (en) | 2023-08-15 |
Family
ID=86963667
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310617749.XA Active CN116380396B (en) | 2023-05-30 | 2023-05-30 | Continuous wind tunnel aircraft atmospheric data system identification test system and method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116380396B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116499699B (en) * | 2023-06-29 | 2023-08-22 | 中国航空工业集团公司沈阳空气动力研究所 | Continuous wind tunnel pressure measurement test data monitoring and correcting method |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5627311A (en) * | 1994-08-23 | 1997-05-06 | National Aerospace Laboratory Of Science & Technology Agency | Transportable three-dimensional calibration wind tunnel system, verification method of flight control system and flight simulator using same |
| CN106840573A (en) * | 2016-12-19 | 2017-06-13 | 中国航天空气动力技术研究院 | A kind of Flush Airdata Sensing System scaling method |
| CN107091725A (en) * | 2016-12-26 | 2017-08-25 | 中国空气动力研究与发展中心高速空气动力研究所 | A kind of measuring method for improving wind tunnel pressure measuring test pressure scan data precision |
| CN111780949A (en) * | 2020-07-10 | 2020-10-16 | 南京航空航天大学 | CFD analysis-based total pressure correction method for high-speed air inlet channel precursor wind tunnel experiment |
| CN112881004A (en) * | 2021-01-11 | 2021-06-01 | 中国空气动力研究与发展中心高速空气动力研究所 | Airspeed tube wind tunnel check test device and check test method |
| CN113281001A (en) * | 2021-04-15 | 2021-08-20 | 南京航空航天大学 | Full-speed domain atmospheric data resolving method based on integrated micro atmospheric data module |
| CN115758818A (en) * | 2022-11-11 | 2023-03-07 | 成都凯天电子股份有限公司 | Ground verification device for aircraft atmospheric data sensor system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9096330B2 (en) * | 2013-08-02 | 2015-08-04 | Honeywell International Inc. | System and method for computing MACH number and true airspeed |
-
2023
- 2023-05-30 CN CN202310617749.XA patent/CN116380396B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5627311A (en) * | 1994-08-23 | 1997-05-06 | National Aerospace Laboratory Of Science & Technology Agency | Transportable three-dimensional calibration wind tunnel system, verification method of flight control system and flight simulator using same |
| CN106840573A (en) * | 2016-12-19 | 2017-06-13 | 中国航天空气动力技术研究院 | A kind of Flush Airdata Sensing System scaling method |
| CN107091725A (en) * | 2016-12-26 | 2017-08-25 | 中国空气动力研究与发展中心高速空气动力研究所 | A kind of measuring method for improving wind tunnel pressure measuring test pressure scan data precision |
| CN111780949A (en) * | 2020-07-10 | 2020-10-16 | 南京航空航天大学 | CFD analysis-based total pressure correction method for high-speed air inlet channel precursor wind tunnel experiment |
| CN112881004A (en) * | 2021-01-11 | 2021-06-01 | 中国空气动力研究与发展中心高速空气动力研究所 | Airspeed tube wind tunnel check test device and check test method |
| CN113281001A (en) * | 2021-04-15 | 2021-08-20 | 南京航空航天大学 | Full-speed domain atmospheric data resolving method based on integrated micro atmospheric data module |
| CN115758818A (en) * | 2022-11-11 | 2023-03-07 | 成都凯天电子股份有限公司 | Ground verification device for aircraft atmospheric data sensor system |
Non-Patent Citations (1)
| Title |
|---|
| 高超声速飞行器大气数据传感系统关键技术初步分析;杨志红;沈作军;;战术导弹控制技术(第01期);第16-19页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116380396A (en) | 2023-07-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110346605B (en) | Method and system for aircraft airspeed calibration based on static pressure error correction | |
| CN116380396B (en) | Continuous wind tunnel aircraft atmospheric data system identification test system and method | |
| CN108051839B (en) | Vehicle-mounted three-dimensional positioning device and three-dimensional positioning method | |
| WO2015113329A1 (en) | On-board combination navigation system based on mems inertial navigation | |
| EP3425363B1 (en) | Systems and methods for correcting acoustic error in pressure sensors | |
| CN106705996A (en) | Aircraft navigation information correcting method based on atmospheric feature parameters | |
| CN112129322B (en) | Method for detecting and correcting installation error of strapdown inertial measurement unit and three-axis rotary table | |
| KR101107219B1 (en) | Method for navigation of an aircraft, intertial navigation system filter using the same, and navigation system using the same | |
| CN106705995A (en) | Calibration method of MEMS gyroscope g value sensitive coefficient | |
| CN117129179B (en) | Mach number correction method for double-support test under continuous wind tunnel wing | |
| CN108090253B (en) | A kind of helicopter digital air system air speed modification method | |
| CN111854691B (en) | Height calculating method and device based on multi-sensor fusion and unmanned aerial vehicle | |
| CN111322984B (en) | Altitude calculating method and device, unmanned aerial vehicle and storage medium | |
| CN114778887B (en) | Unmanned aerial vehicle wind measuring method and device based on improved triangular vector model | |
| CN110186479B (en) | Inertial device error coefficient determination method | |
| CN114858187A (en) | MEMS inertial navigation system calibration method | |
| CN109211271B (en) | Self-correcting method for magnetic compass | |
| CN111122106A (en) | Correction method for measurement error of full-pressure sensor and computer equipment | |
| CN111649903A (en) | Method for indirectly measuring aerodynamic drag coefficient by using dynamic model experiment | |
| EP3598145A1 (en) | Pitot static systems with angle of sideslip determination and compensation | |
| CN113063444B (en) | Sub-angle second precision star sensor optical axis measurement reference deviation calibration method and system | |
| CN117387630A (en) | Automatic direction correction method and system in strong magnetic interference environment | |
| Jain et al. | Calibration and accuracy determination of airdata system for a modern fighter aircraft | |
| CN114543836B (en) | Calibration-free method for dynamically installing attitude reference of ship | |
| CN114383631B (en) | Real-time calibration method based on least square, taylor expansion and comprehensive residual combination |
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 |