CN113609771B - 一种基于mbse的飞机大气参数模型计算方法 - Google Patents
一种基于mbse的飞机大气参数模型计算方法 Download PDFInfo
- Publication number
- CN113609771B CN113609771B CN202110900253.4A CN202110900253A CN113609771B CN 113609771 B CN113609771 B CN 113609771B CN 202110900253 A CN202110900253 A CN 202110900253A CN 113609771 B CN113609771 B CN 113609771B
- Authority
- CN
- China
- Prior art keywords
- pressure
- aircraft
- static pressure
- altitude
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000012821 model calculation Methods 0.000 title claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000004088 simulation Methods 0.000 claims abstract description 10
- 230000003068 static effect Effects 0.000 claims description 38
- 230000000712 assembly Effects 0.000 claims description 11
- 238000000429 assembly Methods 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 5
- 238000005094 computer simulation Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000004891 communication Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 7
- 238000009530 blood pressure measurement Methods 0.000 abstract description 4
- 238000011161 development Methods 0.000 description 3
- 230000002706 hydrostatic effect Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013213 extrapolation Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 208000018634 fetal akinesia deformation sequence Diseases 0.000 description 1
- 208000012165 fetal akinesia deformation sequence syndrome Diseases 0.000 description 1
- 238000002421 fluorescence-activated droplet sorting Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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)
- Evolutionary Computation (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种基于MBSE的飞机大气参数模型计算方法,采用每一个压力测量产品的二组压力进行对比计算出二组参数,再根据大气数据标准文件通过Simulink建立大气参数计算模型计算出大气参数,通过模型仿真和嵌入式测压孔的优势可以计算高高空低压力的大气参数测量和计算。可以测量高度范围达0~4.5km;Vi数范围:0~2000km/h;迎角范围:‑10°~75°;侧滑角范围:‑45°~45°。
Description
技术领域
本发明涉及飞机大气数据参数计算方法,特别涉及一种基于MBSE的飞机大气参数模型计算方法。
背景技术
大气数据系统是飞机上重要的机载电子子系统,其提供的总压、静压、攻角、侧滑角、空速、马赫数和总温等大气参数对于飞行器的飞行安全操纵至关重要。随着现代飞行器对于飞行速度、机动性、舒适性以及隐身性等要求的不断提高,传统大气数据系统所采用的探针式测量方式无法满足新一代飞行器的性能需求,寻求一种新型的大气数据测量方式和解算方法,成为新一代飞行器发展的重要部分。
当飞机大迎角高机动飞行时,传统的安装在飞机头部的外插式大气数据感受器(空速管、风标等),容易诱发飞机结构振动,使飞行员无法工作和搭载的飞行测量、操控设备失效;还可能加剧飞机绕流的非对称特征,引起附加的强烈的侧向力和偏航力矩,可以引起飞机的操控能力和飞行品质降低,甚至导致飞机失控而引发灾难。
对于高超声速飞机而言,外插式大气数据系统产生的激波可能对发动机的进气道的气流产生不利的屏蔽作用,进而导致进气道效率下降甚至丧失;同时由于高超声速飞行,其强烈的气动热效应使外插式大气数据系统的感受部件烧蚀而失效。
这就需要发展一种非插入式即嵌入式大气数据系统,具有实时准确的测量飞机飞行速度、高度、升降速度以及流向角的能力,可以获取多余度的大气数据用于飞行后系统结构重建。
嵌入式大气数据系统(简称FADS)研制是基于在飞机表面适当的位置布置阵列式测压孔提取大气数据的概念。具有在高速度(可达M9)、大迎角(可达75°)和从超低空到高高空条件下实时感知大气数据的突出优点,完全可以满足新一代飞机在高机动、高速度和隐身飞行条件下准确感受大气数据的要求。
以美国为代表的西方航空航天技术发达国家自上个世纪70年代,以天地往返航天工程为背景,开始着手进行系列的理论分析、风洞试验和飞行验证研究,开发出较为成熟的嵌入式大气数据系统产品,被认为是应用空气动力学的8个重大进展之一。
发明内容
为了解决上述技术问题,本发明中披露了一种基于MBSE的飞机大气参数模型计算方法,本发明的技术方案是这样实施的:
一种基于MBSE的飞机大气参数模型计算方法,基于SysML建立飞机系统模型;使用模型仿真方法测试验证大气参数模型;对大气系统的每一个LRU(一种算法)功能进行建模;对数据总线进行数据通信带宽(数据流量)进行仿真。
包括步骤如下,
S1:对每个总压测压孔和静压测压孔建立SysML模型;
S2:使用simulink对大气参数计算方法进行仿真建模;
S3:测量总压信号和静压信号;
S4:将S3中获取的信号转化成数字信号;
S5:计算高度和空速的输入总压和静压;
S6:计算出指示空速;
S7:测量静压测压孔组件的圆周角;
S8:计算迎角和侧滑角;
S9:将计算出的结果传输给其他子系统。
优选地,S4步骤中,通过安装在飞机一个横向切面的四套静压测压孔组件和机头3组总压测压孔组件形成大气数据系统的气压采集通道感受计算所需的6组总压信号和8组静压信号Ps1,Ps2,Ps3,Ps4,Ps5,Ps6,Ps7,Ps8。
优选地,所述静压测压孔组件包括第一组测压孔和第二组测压孔,所述静压测压孔组件包括4个内部蒙皮安装孔。
优选地,所述静压测压孔组件和所述总压测压孔组件结构相同。
优选地,S5步骤中,
解算高度和空速的输入总压
静压
其中,Pt1,Pt2,Pt3,Pt4,Pt5,Pt6为6组总压信号,Ps1,Ps2,Ps3,Ps4,Ps5,Ps6,Ps7,Ps8为8组静压信号。
优选地,S6步骤中,气压高度Hp=Hr+(Pb-Ps/4)*Tt,其中,Hr是高度分层的高度均值,Pb是飞机相对机场高度的高度H相对的静压压力值,Tt是不同高度层的修正系数,P0是海平面大气压力标准值,Vn2为海平面标准声速,Pt为解算高度和空速的输入总压,
中间变量
根据公式计算出指示空速Vi。
优选地,S8步骤中侧滑角用Beta标记,迎角Aoa标记;
Aoa=(sinPr1+cosPr3)-(sinPr2+cosPr4);Pr1,Pr2,Pr3和Pr4为S7步骤中测量出的四组圆周角。
本发明主要是采用多压力源,采用建模仿真,进行压力前期修正,采用压差原理进行攻角和侧滑角计算。多压力源就是利用神经网络的方法对大气压力源进行修正,通过SysML建模理论进行建模,对各产品进行建模,对大气数据进行仿真计算提供压力输入,采用非线性映射方法模拟大气压力对计算函数的重复验证,等出大气参数就算公式。
采用每一个压力测量产品的二组压力进行对比计算出二组参数,再根据大气数据标准文件通过Simulink建立大气参数计算模型计算出大气参数,通过模型仿真和嵌入式测压孔的优势可以计算高高空低压力的大气参数测量和计算。可以测量高度范围达0~4.5km;Vi数范围:0~2000km/h;迎角范围:-10°~75°;侧滑角范围:-45°~45°。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一种实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
其中相同的零部件用相同的附图标记表示。需要说明的是,下面描述中使用的词语“前”、“后”、“左”、“右”、“上”和“下”指的是附图中的方向,词语“底面”和“顶面”、“内”和“外”分别指的是朝向或远离特定部件几何中心的方向。
图1为本发明的模型解算图;
图2为静压测压孔组件和总压测压孔组件的测压孔结构图;
图3为大气参数仿真建模算法流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例
在一种具体的实施例中,如图1、图2和图3所示,一种基于MBSE的飞机大气参数模型计算方法,基于SysML建立飞机系统模型;使用模型仿真方法测试验证大气参数模型;对大气系统的每一个LRU(一个独立的产品单元)进行功能建模;对数据总线进行数据通信带宽(数据流量)进行仿真。
本实施例中,在仿真飞机模型的前机身同一切面的上下左右布置4组静压测压孔组件,机头侧前方左右和右上方各布置一组总压测压孔组件。
由于仿真模型的飞机型号不同,具体位置根据不同飞机型号安装的具体位置也不同,且该部分为现有技术。
所述静压测压孔组件包括第一组测压孔和第二组测压孔,所述静压测压孔组、件包括4个内部蒙皮安装孔。所述静压测压孔组件和所述总压测压孔组件结构相同。即每个静压测压孔组件和所述总压测压孔组件可测出两组数据,本实施例中静压测压孔组件和总压测压孔组件采用嵌入式安装。
采用每个静压测压孔组件和所述总压测压孔组件的二组压力进行对比计算出总压信号和静压信号的参数,再根据大气数据标准文件通过Simulink建立大气参数计算模型计算出大气参数,通过模型仿真和嵌入式测压孔的优势可以计算高高空低压力的大气参数测量和计算。可以测量高度范围达0~4.5km;Vi数范围:0~2000km/h;迎角范围:-10°~75°;侧滑角范围:-45°~45°。迎角代表飞机的飞行方向的地平线的夹角,侧滑角就是飞机的机头方向与飞机实际前进方向的夹角。
本实施例主要是采用多压力源,采用建模仿真,进行压力前期修正,采用压差原理进行迎角和侧滑角计算。多压力源就是利用神经网络的方法对大气压力源进行修正,通过SysML建模理论进行建模,对各产品进行建模,对大气数据进行仿真计算提供压力输入,采用非线性映射方法模拟大气压力对计算函数的重复验证,得出大气参数就算公式。
本实施例包括步骤如下,
S1:对安装在飞机一个横向切面的四套静压测压孔组件和机头3组总压测压孔组件建立SysML模型;测压孔如图2所示,
S2:使用simulink对大气参数计算方法进行仿真建模;
S3:通过安装在飞机一个横向切面的四套静压测压孔组件和机头3组总压测压孔组件形成大气数据系统的气压采集通道感受计算所需的6组总压信号和8组静压信号;每个静压测压孔组件和所述总压测压孔组件各测得二组压力值;于是测得数据:Pt1,Pt2,Pt3,Pt4,Pt5,Pt6为6组总压信号,Ps1,Ps2,Ps3,Ps4,Ps5,Ps6,Ps7,Ps8为8组静压信号。
S4:将S3中获取的信号转化成数字信号;通过测压孔集成在一起的解算器把压力信号转换成数字信号。
S5:通过驻留在综合集成航电平台的大气数据解算驻留应用计算器飞机需要的高度和空速的输入总压和静压参数;
解算高度和空速的输入总压
静压
S6:计算出指示空速;
气压高度Hp=Hr+(Pb-Ps/4)*Tt,其中,Hr是高度分层的高度均值,Pb是飞机相对机场高度的高度H相对的静压压力值,Tt是不同高度层的修正系数,P0是海平面大气压力标准值,Vn2为海平面标准声速,Pt为解算高度和空速的输入总压,
中间变量
根据公式计算出指示空速Vi。
S7:根据安装位置的不同,测量出四套静压测压孔组件中压力传感器的圆周角Pr1,Pr2,Pr3和Pr4;
S8:根据飞机的地坐标系和测压孔相对飞机坐标系的值修正根据压差计算出来的侧滑角和迎角。侧滑角用Beta标记,迎角Aoa标记;
S9:将计算出的结果传输给其他子系统。
如图3所示,场压数据就是机场的PS值,也就是机场的静压压力值,M是马赫数。是一个比值,无量纲参数。VHn是第n次的H相对的速度值,V为飞机的实时空速,基于本实施例得到的大气数据系统参数进行二次计算可得出图3所示的数据,这些计算内容可采用其他标准公式计算,为现有技术,因此本发明不做赘述。
需要指出的是,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于MBSE的飞机大气参数模型计算方法,其特征在于:基于SysML建立飞机系统模型;使用模型仿真方法测试验证大气参数模型;对大气系统的每一个LRU进行功能建模;对数据总线进行数据通信带宽仿真;
所述方法包括步骤如下,
S1:对每个总压测压孔和静压测压孔建立SysML模型;
S2:使用simulink对大气参数的飞机相关系统进行仿真建模;
S3:测量总压信号和静压信号;
S4:将S3中获取的信号转化成数字信号;过安装在飞机一个横向切面的四套静压测压孔组件和机头3组总压测压孔组件形成大气数据系统的气压采集通道感受计算所需的6组总压信号和8组静压信号Ps1,Ps2,Ps3,Ps4,Ps5,Ps6,Ps7,Ps8;
S5:计算高度和空速的输入总压和静压;
解算高度和空速的输入总压
静压
其中,Pt1,Pt2,Pt3,Pt4,Pt5,Pt6为6组总压信号,Ps1,Ps2,Ps3,Ps4,Ps5,Ps6,Ps7,Ps8为8组静压信号
S6:计算出指示空速;
气压高度Hp=Hr+(Pb-Ps/4)*Tt,
其中,Hr是高度分层的高度均值,Pb是飞机相对机场高度的高度H相对的静压压力值,Tt是不同高度层的修正系数,P0是海平面大气压力标准值,Vn2为海平面标准声速,Pt为解算高度和空速的输入总压,
中间变量
根据公式计算出指示空速Vi;
S7:测量静压测压孔组件的圆周角;
S8:计算迎角和侧滑角;
侧滑角用Beta标记,迎角Aoa标记;
Aoa=(sinPr1+cosPr3)-(sinPr2+cosPr4);
Pr1,Pr2,Pr3和Pr4为S7步骤中测量出的四组圆周角;
S9:将计算出的结果传输给其他子系统。
2.根据权利要求1所述的一种基于MBSE的飞机大气参数模型计算方法,其特征在于:所述静压测压孔组件包括第一组测压孔和第二组测压孔,所述静压测压孔组件包括4个内部蒙皮安装孔。
3.根据权利要求2所述的一种基于MBSE的飞机大气参数模型计算方法,其特征在于:所述静压测压孔组件和所述总压测压孔组件结构相同。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110900253.4A CN113609771B (zh) | 2021-08-06 | 2021-08-06 | 一种基于mbse的飞机大气参数模型计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110900253.4A CN113609771B (zh) | 2021-08-06 | 2021-08-06 | 一种基于mbse的飞机大气参数模型计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113609771A CN113609771A (zh) | 2021-11-05 |
CN113609771B true CN113609771B (zh) | 2024-04-12 |
Family
ID=78339584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110900253.4A Active CN113609771B (zh) | 2021-08-06 | 2021-08-06 | 一种基于mbse的飞机大气参数模型计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113609771B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112685972A (zh) * | 2020-12-29 | 2021-04-20 | 南京航空航天大学 | 一种高效率的高超声速飞行器压力传感器分布优化方法 |
-
2021
- 2021-08-06 CN CN202110900253.4A patent/CN113609771B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112685972A (zh) * | 2020-12-29 | 2021-04-20 | 南京航空航天大学 | 一种高效率的高超声速飞行器压力传感器分布优化方法 |
Non-Patent Citations (1)
Title |
---|
先进战斗机过失速机动大气数据融合估计方法;杨朝旭等;《航空学报》;正文第1-3章 * |
Also Published As
Publication number | Publication date |
---|---|
CN113609771A (zh) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108152529A (zh) | 一种基于飞行参数计算风速及风向的方法 | |
US6604029B2 (en) | Multi-function air data probes using neural network for sideslip compensation | |
Cumming et al. | Aerodynamic flight test results for the adaptive compliant trailing edge | |
AU2014339084A1 (en) | Air data sensor for an aircraft | |
Raab et al. | In-flight testing of MEMS pressure sensors for flight loads determination | |
RU2341775C1 (ru) | Способ определения аэродинамического угла летательного аппарата | |
CN113609771B (zh) | 一种基于mbse的飞机大气参数模型计算方法 | |
CN113525711B (zh) | 通过飞行试验辨识飞机气动焦点的方法 | |
CN112764424B (zh) | 一种飞机飞行控制系统关键传感器故障重构方法 | |
Li et al. | A compact embedded flight parameter detection system for small soaring UAVs | |
US10296013B2 (en) | Vehicle guidance system and method that uses air data from surface-mounted pressure sensors for vehicle orientation control | |
Nebula et al. | Virtual air data: A fault-tolerant approach against ADS failures | |
Wan et al. | Fault-tolerant FADS system development for a hypersonic vehicle via neural network algorithms | |
Howe et al. | Quiet Spike Near Field Flight Test Pressure Measurements with CFD Comparisons | |
CN112762960A (zh) | 一种飞行器所处风场的在线计算方法 | |
Jain et al. | Calibration and accuracy determination of airdata system for a modern fighter aircraft | |
Jain et al. | A novel methodology to update table in air data system of a high performance fighter aircraft | |
CN118579273B (zh) | 适用于旋翼机的大气数据探头安装位置选定方法 | |
Weishaeupl et al. | Flight Test Driven Development of Low Cost UAVs-Pitfalls and Opportunities | |
Corrêa et al. | Helicopter air data systems calibration using DGPS | |
Soldatkin | Design Features and Efficiency Analysis of the Aircraft Complex Air Data System with a Fixed Flush-Mounted Flow Receiver | |
Raab | Practical examples for the flight data compatibility check | |
Ellingson et al. | A combined experimental and numerical analysis of UAV Pitot-static system error at low Reynolds number | |
Kelbe et al. | Loads Certification of the Cirrus Aircraft SF50 Vision Jet® | |
CN117332711A (zh) | 一种飞翼布局的大气数据压力建模方法 |
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 |