CN113837639B - 一种基于多维指标的航空器遭遇尾流风险评估系统 - Google Patents
一种基于多维指标的航空器遭遇尾流风险评估系统 Download PDFInfo
- Publication number
- CN113837639B CN113837639B CN202111154614.1A CN202111154614A CN113837639B CN 113837639 B CN113837639 B CN 113837639B CN 202111154614 A CN202111154614 A CN 202111154614A CN 113837639 B CN113837639 B CN 113837639B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- risk
- wing
- wake
- index
- 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
- 238000012502 risk assessment Methods 0.000 title claims abstract description 10
- 238000011156 evaluation Methods 0.000 claims abstract description 62
- 238000005259 measurement Methods 0.000 claims abstract description 58
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000005096 rolling process Methods 0.000 claims abstract description 32
- 230000008569 process Effects 0.000 claims abstract description 13
- 238000005303 weighing Methods 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 32
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 230000006698 induction Effects 0.000 claims description 5
- 230000035484 reaction time Effects 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 2
- 230000009471 action Effects 0.000 claims description 2
- 230000015556 catabolic process Effects 0.000 claims description 2
- 238000006731 degradation reaction Methods 0.000 claims description 2
- 230000007480 spreading Effects 0.000 claims 1
- 238000013459 approach Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- 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
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种基于多维指标的航空器遭遇尾流风险评估系统。本发明系统用于执行以下步骤：S1，获取飞机遭遇尾涡风险的衡量指标的值，衡量指标包括滚动力矩系数、过载增量、高度下降率和滚转坡度角；S2，建立各个衡量指标的危险等级评估表，并基于衡量指标的值，获取飞机各个衡量指标的危险等级；S3，建立评语集，并基于各个衡量指标的危险等级和评语集，获取各个衡量指标的评分；S4，采用层次分析法对衡量指标进行分析，得到飞机各个衡量指标对应的指标权重；S5，构建飞机的评价函数，并带入飞机各个衡量指标的危险等级评分和指标权重，计算评估值，评估值用于评估飞机遭遇尾流风险。具有能够量化评估后机遭遇尾涡危险风险的技术效果。
Description
技术领域
本发明涉及尾涡风险评价领域，特别涉及一种基于多维指标的航空器遭遇尾流风险评估系统。
背景技术
尾涡是飞机产生升力的一种副产品，表现为在飞机后方卷起形成一对向外向下旋转的涡流，并在飞机后方持续几分钟时间。尾涡的强度由产生尾涡的飞机重量、飞行速度机翼形状和翼展长度所决定，其中最主要的是飞机的重量。尾涡强度随飞机重量、载荷因数的增加和飞行速度的减小而增大。参照图1，处在另一架飞机尾涡区域内的飞机可能会造成空气动力学紊乱，可能无法通过自身的控制来恢复姿态，进一步出现下沉、发动机停止运转以及严重的滚转，造成乘客人身安全问题和飞机的损害。
在尾涡遭遇的安全性评估方面，Carsten等对飞机遭遇严重性评估研究中采用了气动响应模型的方法。Gerben等人用尾涡诱导力矩和飞机滚转阻尼力矩的大小来反映飞机遭遇尾涡后的滚转运动。Campos提出通过建立有关力矩模型来模拟飞机在尾涡影响下的滚转状况，并以最大滚转角速度作为评估风险的指标。Visscher提出了诱导力矩系数这一衡量指标，用来评估飞机在遭遇尾涡后的严重性。欧洲空中航行安全组织提出了一项名为“RECAT-EU”的提案，它包括对国际民航组织目前关于进近和起飞的尾流安全间隔最小值规定的修订。冯志勇分析了特定天气条件下、一定时间间隔后残存尾涡的强度和后机能承受的尾涡强度进行比较，引入滚转比例系数对遭遇强度进行评估。
发明内容
目前尚缺乏一种量化评估后机遭遇尾涡危险风险的方法，因此本发明提供一种基于多维指标的航空器遭遇尾流风险评估系统。
为了实现上述发明目的，本发明提供了以下技术方案：
一种基于多维指标的航空器遭遇尾流风险评估系统，所述系统用于执行以下步骤：
S1，获取飞机遭遇尾涡风险的衡量指标的值，所述衡量指标包括滚动力矩系数、过载增量、高度下降率和滚转坡度角；
S2，建立各个衡量指标的危险等级评估表，并基于衡量指标的值，获取飞机各个衡量指标的危险等级；
S3，建立评语集，并基于各个衡量指标的危险等级和评语集，获取各个衡量指标的评分；所述评语集包括1、2、3、4、5五个危险等级，分别对应2、4、6、8、10五个评分；
S4，采用层次分析法对衡量指标进行分析，得到飞机各个衡量指标对应的指标权重；
S5，将所述衡量指标的评分和各个衡量指标对应的指标权重代入飞机的评价函数，计算出评估值，所述评估值用于评估飞机遭遇尾流风险。
通过采用上述技术方案，通过选取滚转力矩系数、过载增量、高度下降率和滚转角速度四个衡量指标来评价飞机遭遇尾涡风险。首先利用层次分析法建立航空器遭遇尾流风险评价指标体系，获取衡量指标的指标权重。获取衡量指标的值，并基于衡量指标和指标权重构建评价函数，对航空器遭遇前机尾流的风险进行评估。该方法可以实时动态地监视不同航段上航空器遭遇尾流风险态势，并对避免尾流危险遭遇，提高进近阶段运行安全，协助管制员实时决策提供了参考。
可选地，步骤S5中，评估值用于评估飞机遭遇尾流风险具体为判断评估值对应的危险程度；
若评估值小于2，则其危险程度为安全；
若评估值大于等于2小于4，则其危险程度为较为安全；
若评估值大于等于4小于6，则其危险程度为轻度危险；
若评估值大于等于6小于8，则其危险程度为中度危险；
若评估值大于等于8小于10，则其危险程度为强烈危险。
可选地，在获取飞机的滚动力矩系数时，将飞机机翼分割为条带状，形成多个微元机翼，再计算滚动力矩系数；所述计算飞机的滚动力矩系数采用下述公式计算，
LR＝L×F；
式中，RMC为滚转力矩系数，V为空气的来流速度；S为机翼面积，B为前机翼展；F为机翼上其中一个微元机翼升力的变化量，LR为该微元机翼受诱导力矩大小，L为该微元机翼到机翼中心处的距离。
可选地，所述过载增量采用下述公式计算，
式中，Fall为飞机所受合力，M为飞机重量，g为重力加速度；其中，
Fall＝ΔFwing+ΔFtail+ΔFengine+ΔFfuselage；
式中，ΔFwing为机翼升力变化量；ΔFtail为尾翼升力变化量；ΔFengine为发动机升力变化量；ΔFfuselage为机身升力变化量。
可选地，所述高度下降率采用下述公式计算，
其中，Gr为高度下降率，ΔH为高度变化量；
式中，tm为飞行员反应时间；tn为航空器响应时间，M为飞机重量，Fall为飞机所受合力。
可选地，计算滚转坡度角时，采用滚转角速度计算公式计算滚转角速度，再对所述滚转角速度进行积分，得到滚转坡度角，所述滚转角速度计算公式为：
式中，为滚转角速度，WRp为滚转阻尼力矩系数；WRc为滚转操纵力矩系数；H(t′-t″)为单位阶跃函数；t′为遭遇尾涡时间，t″为飞行员与飞机反应时间；A为展弦比；I为飞机转动惯量，通过无量纲回转半径来估算飞机相对于中心惯性主轴的惯性矩。
可选地，步骤S4具体包括以下步骤，
101，根据所述衡量指标构建四阶判断矩阵；
102，对所述四阶判断矩阵进行一致性检验；当一致性检验的结果时，判断矩阵通过一致性检验，则执行步骤103；
103，基于判断矩阵，采用均方根法计算指标权重。
可选地，所述一致性检验采用下述公式，
式中，CI为判断矩阵的一致性判断指标；λmax为判断矩阵的最大特征值；n为判断矩阵的阶数。
可选地，所述均方根法采用下述公式，
式中，aij为判断矩阵中的i行、j列的元素。
可选地，所述评价函数采用线性加权法，所述线性加权法采用如下公式，
式中，z为评估值；C(ai)为飞机各衡量指标的评分；ωi为各个衡量指标对应的指标权重。
与现有技术相比，本发明的有益效果：通过选取滚转力矩系数、过载增量、高度下降率和滚转角速度四个衡量指标来评价飞机遭遇尾涡风险。首先利用层次分析法建立航空器遭遇尾流风险评价指标体系，获取衡量指标的指标权重。获取衡量指标的值，并基于衡量指标和指标权重构建评价函数，对不同翼型的航空器遭遇前机尾流的风险进行评估。该方法可以实时动态地监视不同航段上航空器遭遇尾流风险态势，并对避免尾流危险遭遇，提高进近阶段运行安全，协助管制员实时决策提供了参考。
附图说明：
图1为本发明中飞机尾涡遭遇示意图；
图2为本发明的方法流程图；
图3为本发明中条状带模型示意图；
图4为本发明中条状带模型部分示意图。
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例，凡基于本发明内容所实现的技术均属于本发明的范围。
实施例1
一种基于多维指标的航空器遭遇尾流风险评估系统，参照图1，系统用于执行下述步骤，
S1，获取飞机遭遇尾涡风险的衡量指标的值，衡量指标包括滚动力矩系数、过载增量、高度下降率和滚转坡度角。
求取衡量指标的值之前，需要先判断飞机类型，飞机类型包括圆翼飞机和后掠翼飞机；飞机的机型不同，则在评估飞机遭遇尾流风险时的评估过程也不一致，在本实施例中，主要是体现在衡量指标的数据不同。
在获取飞机的滚动力矩系数时，将飞机机翼分割为条带状，形成多个微元机翼，再计算滚动力矩系数；计算飞机的滚动力矩系数采用下述公式，
LR＝L×F；
其中，RMC为滚转力矩系数，V为空气的来流速度，与飞机的飞行速度一致；S为机翼面积，B为前机翼展；F为机翼上某个微元机翼升力的变化量，LR为该微元机翼受诱导力矩大小，L为该微元机翼到机翼中心处的距离。
在计算滚转力矩系数时，由于机翼各处受力不均，因此可以采用微元思想，将机翼分割成无数条状带，将每条条带状受力当成一个整体，之后按机翼展展长方向进行积分，即可以得到机翼整体的升力变化量。条状带模型图如图3和图4所示，每条条带状上升力变化量乘以其到相应机翼中心的距离，即为后机机翼该条状带区域在前机尾涡诱导作用下的诱导滚转力矩，之后进行累加，即为机翼整体的滚转力矩。
由条状带法思想，我们计算一个条带上由尾涡引起的升力变化量，进而得到此条带上的诱导力矩。
式中，dF(y)表示为局部的升力变化量；ΔLR为局部诱导力矩；ρ为尾涡所处位置的空气密度；V为来流速度；cl(y)为飞机处理y位置处的有效升力线系数；c(y)为弦长；Δα(y)表示迎角变化量。
针对两种飞机类型；椭圆形翼弦弦长cl(y)采用下述公式计算：
后掠翼飞机翼弦弦长cl(y)采用下述公式计算：
式中，cr为翼根处的弦长；λ为梢根比。
因此，椭圆形翼弦飞机的滚转力矩系数为：
后掠翼飞机飞机的滚转力矩系数为：
尾涡的强度大小用环量加以衡量，其中，
rc＝0.052b0；
式中，Γ0为尾涡刚从机翼后缘脱落时的环量大小，即初始环量；b0为初始涡核间距；rc为初始涡核半径；M为前机质量大；B为前机翼展；ρ为尾涡所处位置的空气密度；g为重力加速度；V1为前机真空速；ω0为特征速度；t0为参考时间，即尾涡在相互诱导的初始下降速度即以特征速度ω0下降距离b0所需时间。
尾涡耗散速度的快慢与湍流耗散率有关，tc为近场尾涡的持续时间，根据Sarpkaya公式tc与无因次尾涡耗散率ε*之间的关系，如下式：
又有，
当ε*≥0.2535时，
当0.0235≥ε*≥0.0121时，
当0.0121≥ε*≥0.001时，
tc＝9.18t0-180t0ε*；
当0.001≥ε*≥时，
tc＝9t0；
尾涡在经过卷起，形成两个方向涡流后，出现线性不稳定，涡量开始径向扩散，之后不稳定性导致两个涡流之间的强烈相互作用，最终崩溃耗散成大气中湍流。把耗散过程分为两个阶段，即扩散阶段和快速衰减阶段；
扩散阶段耗散公式如下：
快速快速衰减阶段公式如下：
且前机产生的尾涡会对后机产生一个诱导速度，根据Burnham-Hallock模型，可得诱导速度采用下式计算：
式中，Γ为尾涡经过一定时间耗散之后的环量，r为该点到涡心处的距离，v为该点的垂直诱导速度。
前机产生的尾涡在后机机翼上每一点的切向速度都不相同，因此每一点的升力变化量也不相同。当后机从前机尾涡场产生的两涡之间任意一点位置进入时，由于左右机翼受力不均，会以两翼中点为中心产生一个诱导滚转力矩，使飞机发生滚转，严重影响飞行安全。
在ICAO标准中，采用滚转力矩系数这一无量纲化参数去比较具有不同物理和空气动力学特性的飞机之间遭遇尾涡的严重性程度,当滚转力矩系数在起飞与进近阶段小于0.04时、平飞阶段小于0.06时，可认为飞机安全。
在求过载增量时，采用下述公式，
式中，Fall为飞机所受合力，M为飞机重量，g为重力加速度；其中，
Fall＝ΔFwing+ΔFtail+ΔFengine+ΔFfuselage；
式中，ΔFwing为机翼升力变化量；ΔFtail为尾翼升力变化量；ΔFengine为发动机升力变化量；ΔFfuselage为机身升力变化量；其中，
式中，ρ为大气密度，V为来流速度，即飞机真空速；cl(y)为升力线斜率；c(y)为飞机翼展长坐标处的翼弦长；Vz(y)为前机尾涡在后机机翼上诱导产生的上洗或下洗速度；
式中，s(y)为尾翼弦长；C为尾翼展长；
式中，b(y)为发动机翼弦长；D为发动机展长；
ΔFfuselage＝Ncosα-Asinα；
式中，N为法向力，A为轴向力；α为迎角。
在求高度下降率时，采用下述公式计算；
其中，Gr为高度下降率，ΔH为高度变化量；
式中，tm为飞行员反应时间；tn为航空器响应时间。
计算滚转坡度角时，首先采用下述公式计算滚转角速度，
式中，为滚转角速度，WRp为滚转阻尼力矩系数，取值为-0.5；WRc为滚转操纵力矩系数；H(t′-t″)为单位阶跃函数；t′为遭遇尾涡时间，t″为飞行员与飞机反应时间；A为展弦比；I为飞机转动惯量，通过无量纲回转半径来估算飞机相对于中心惯性主轴的惯性矩；
I采用下式计算：
式中，为无量纲回转半径，与发动机的类别和飞机构型有关，从下表获取：
可以得出，椭圆形翼弦的飞机，其滚转角速度计算公式如下：
后掠翼飞机，其滚转角速度计算公式如下：
再对滚转角速度进行积分，得到滚转坡度角。
对于椭圆形翼弦的飞机，其滚转坡度角计算公式如下：
对于椭圆形翼弦的飞机，其滚转坡度角计算公式如下：
S2，建立各个衡量指标的危险等级评估表，并基于衡量指标的值，获取飞机各个衡量指标的危险等级。
其中，滚动力矩系数、过载增量、高度下降率和滚转角速度的危险等级评估表如下所示，危险等级评估表基于衡量指标的干扰强度划分各个衡量指标的危险等级，危险等级包括1、2、3、4、5五个等级。
其中，滚动力矩系数干扰强度对应的危险等级评估如下表所示：
| 干扰强度 | 滚转力矩系数 | 危险等级 |
| 无干扰 | RMC≤0.02 | 1 |
| 轻度干扰 | 0.02≤RMC≤0.04 | 2 |
| 较大干扰 | 0.04≤RMC≤0.08 | 3 |
| 强烈干扰 | 0.08≤RMC≤0.1 | 4 |
| 剧烈干扰 | 0.1≤RMC | 5 |
过载增量颠簸强度对应的危险等级评估如下表所示：
| 颠簸强度 | 过载增量 | 危险等级 |
| 无湍流 | Δn≤0.15 | 1 |
| 轻度颠簸 | 0.15≤Δn≤0.5 | 2 |
| 中度颠簸 | 0.5≤Δn≤0.8 | 3 |
| 较大颠簸 | 0.8≤Δn≤1 | 4 |
| 强烈颠簸 | 1≤Δn | 5 |
高度下降率干扰强度对应的危险等级评估如下表所示：
滚转坡度角干扰强度对应的危险等级评估如下表所示：
| 干扰强度 | 坡度角 | 危险等级 |
| 无安全影响 | Φ＜10° | 1 |
| 轻度影响 | 10°≤Φ＜25° | 2 |
| 较大影响 | 25°≤Φ＜30° | 3 |
| 存在危险 | 30°≤Φ＜45° | 4 |
| 强烈危险 | 45°≤Φ | 5 |
S3，建立评语集，并基于各个衡量指标的危险等级和评语集，获取各个衡量指标的评分；评语集包括1、2、3、4、5五个危险等级，分别对应2、4、6、8、10五个评分。
如下表所示：
| 危险等级 | 对应评分 |
| 1 | 2 |
| 2 | 4 |
| 3 | 6 |
| 4 | 8 |
| 5 | 10 |
S4，采用层次分析法对衡量指标进行分析，得到飞机各个衡量指标对应的指标权重。
层次分析法包括下述步骤，
101，根据衡量指标构建四阶判断矩阵；
对航空器遭遇尾流风险指标体系中相关因素指标的相对重要性进行两两比较判断，构造判断矩阵，采用T.L.Saaty建议的1～9标度法。
标度表如下：
准则层包括滚转力矩系数(a1)、过载增量(a2)、高度下降率(a3)、滚转坡度角(a4)四个子项，得到四阶判断矩阵如下。
其中，aij表示元素ai和aj两两相比得到的数值。
102，对四阶判断矩阵进行一致性检验；当一致性检验的结果时，判断矩阵通过一致性检验，则执行步骤103；
在建立判断矩阵时，由于客观事物具有复杂性，主题认识具有主观性和多样性，判断矩阵不可能具有完全一致性，因此需要进行判断矩阵的一致性检验。
一致性检验采用下述公式，
式中，CI为判断矩阵的一致性判断指标；λmax为判断矩阵的最大特征值；n为判断矩阵的阶数。
且当判断矩阵的阶数较大时，则引入随机一致性指标RI(Random Index)进行修正，经修正的一致性指标用CR表示。
式中，RI取值参照下表，n为判断矩阵的阶数。
| n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
| RI | 0 | 0 | 0.58 | 0.90 | 1.12 | 1.24 | 1.32 | 1.41 | 1.45 | 1.49 |
当求出的结果CR＜0.1时，判断矩阵通过一致性检验。
103，基于判断矩阵，采用均方根法计算指标权重。
均方根法计算指标权重采用如下公式，
/>
式中，aij为判断矩阵中的i行、j列的元素。
得出四个评价指标所对应的指标权重ω1、ω2、ω3、ω4。
S5，将衡量指标的评分和各个衡量指标对应的指标权重代入飞机的评价函数，计算出评估值，评估值用于评估飞机遭遇尾流风险。
本实施例中，采用线性加权法构建评价函数，线性加权法如以下公式，
式中，z为评估值；C(ai)为衡量指标的危险等级评分；ωi为衡量指标的指标权重。
分别代入两种类型飞机衡量指标的值，得到两种类型飞机的评价函数。
后机与前机的相对位置可由三维坐标系(x,y,z)来表示，后机在前机尾涡影响范围内的任意位置处风险均要考虑滚转力矩系数(RMC)、过载增量(|Δn|)、高度下降率(Gr)和滚转角加速度(Φ)这四种衡量航空器危险程度的指标。
根据民用航空器分为圆翼和后掠翼飞机，可设计评价函数为：
z＝ω1C(RMC)+ω2C(Δn)+ω3C(Gr)+ω4C(Φ)。
求取评估值后，将评估值用于评估飞机遭遇尾流风险，判断评估值对应的危险程度；
若评估值小于2，则其危险程度为安全；
若评估值大于等于2小于4，则其危险程度为较为安全；
若评估值大于等于4小于6，则其危险程度为轻度危险；
若评估值大于等于6小于8，则其危险程度为中度危险；
若评估值大于等于8小于10，则其危险程度为强烈危险。
如下表所示，
| 危险程度 | 评价函数分值 | 危险等级 |
| 安全 | z＜2 | 1 |
| 较为安全 | 2≤z＜4 | 2 |
| 轻度危险 | 4≤z＜6 | 3 |
| 中度危险 | 6≤z＜8 | 4 |
| 强烈危险 | 8≤z＜10 | 5 |
通过查阅上表，工作人员能够根据求取的评估值来实时动态地监视不同航段上航空器遭遇尾流风险态势，并对避免尾流危险遭遇，提高进近阶段运行安全，为协助管制员实时决策提供参考。
两种类型飞机中，四个衡量指标的计算公式不同，则其评价函数也不一致，分别构建评价函数后；能够使用评价函数来将滚转力矩系数、过载增量、高度下降率和滚转坡度角四个衡量指标对航空器遭遇尾流风险进行综合评价，对不同翼型的航空器遭遇前机尾流的风险进行评估。
实施例2
在本实施例中，选取国产ARJ21(中型机)飞机跟随A330-200(重型机)飞机对实施例1中所述的方法进行验证，ICAO规定近地阶段中型机跟随重型机时尾流安全间隔为9.3km。本实施例选取进近阶段时后机穿越前机尾涡流场的三个位置：
1)飞机纵向穿越单涡涡核中心，图1中B位置；
2)飞机纵向穿越距两涡核连线中点处左右各10m处位置；
3)穿越两涡涡核连线中点处。
选取以上三个位置来判断飞机此时遭遇尾流的风险程度，判断ICAO间隔标准是否安全。
其中，ARJ21参数与A330-200参数如下表所示：
尾流及环境参数如下表所示：.
将上述参数带入到实施例1中得到的评价参数公式中，能够计算得到如下表所示的计算结果：
再带入到危险等级划分表中，得到衡量指标的得分结果如下表所示：
再根据，层次分析法得出ω1＝ω2＝ω3＝ω4，四种情况同样重要。
因此对于纵向穿越单涡涡核中心这种情况：
可以得到，评价函数为：
z＝0.25×1.642+0.25×0.06+0.25×0.742+0.25×1.42；
z＝0.966；
对于纵向穿越两涡涡核连线中心这种情况：
z＝0.25×0+0.25×0.248+0.25×0.309+0.25×0；
z＝0.140；
对于纵向穿越两涡连线中心左右各10m处的这种情况：z＝0.25×0.295+0.25×0.32+0.25×0.394+0.25×0.54；
z＝0.387；
根据上述计算结果，可见在ICAO间隔为9.3km时，三种不同的后机穿越前机尾涡流场的穿越方式模糊综合评价得分均小于2，评估为安全。三种方式比较来看，穿越两涡涡核连线中心这种方式最安全，其次是穿越两涡连线中心左右各10m处；穿越单涡涡核中心这种情况危险等级高于以上两种情况。根据上述结论，本评估系统能够为管制员实时决策提供参考。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于多维指标的航空器遭遇尾流风险评估系统，其特征在于，所述系统用于执行以下步骤：
S1，获取飞机遭遇尾涡风险的衡量指标的值，所述衡量指标包括滚动力矩系数、过载增量、高度下降率和滚转坡度角；
S2，建立各个衡量指标的危险等级评估表，并基于衡量指标的值，获取飞机各个衡量指标的危险等级；
其中，滚动力矩系数干扰强度对应的危险等级评估如下表所示：
；
过载增量颠簸强度对应的危险等级评估如下表所示：
；
高度下降率干扰强度对应的危险等级评估如下表所示：
滚转坡度角干扰强度对应的危险等级评估如下表所示：
S3，建立评语集，并基于各个衡量指标的危险等级和评语集，获取各个衡量指标的评分；所述评语集包括1、2、3、4、5五个危险等级，分别对应2、4、6、8、10五个评分；
S4，采用层次分析法对衡量指标进行分析，得到飞机各个衡量指标对应的指标权重；
S5，将所述衡量指标的评分和各个衡量指标对应的指标权重代入飞机的评价函数，计算出评估值，所述评估值用于评估飞机遭遇尾流风险；
在获取飞机的滚动力矩系数时，将飞机机翼分割为条带状，形成多个微元机翼，再计算滚动力矩系数；所述计算滚动力矩系数采用下述公式计算，
LR＝L×F；
式中，RMC为滚转力矩系数，V为空气的来流速度；S为机翼面积，B为前机翼展；F为机翼上其中一个微元机翼升力的变化量，LR为该微元机翼受诱导力矩大小，L为该微元机翼到机翼中心处的距离；ρ为尾涡所处位置的空气密度；
将机翼分割成无数条状带，将每条条带状受力当成一个整体，之后按机翼展展长方向进行积分，即可以得到机翼整体的升力变化量；每条条带状上升力变化量乘以其到相应机翼中心的距离，即为后机机翼该条状带区域在前机尾涡诱导作用下的诱导滚转力矩，之后进行累加，即为机翼整体的滚转力矩；
一个条带上由尾涡引起的升力变化量计算公式为：
式中，dF(y)表示为局部的升力变化量；ΔLR为局部诱导力矩；；V为来流速度；cl(y)为飞机所处y位置处的有效升力线系数；c(y)为弦长；Δα(y)表示迎角变化量；
所述过载增量采用下述公式计算，
式中，Fall为飞机所受合力，M为飞机重量，g为重力加速度；其中，
Fall＝ΔFwing+ΔFtail+ΔFengine+ΔFfuselage；
式中，ΔFwing为机翼升力变化量；ΔFtail为尾翼升力变化量；ΔFengine为发动机升力变化量；ΔFfuselage为机身升力变化量；
所述高度下降率采用下述公式计算，
其中，Gr为高度下降率，ΔH为高度变化量；
式中，tm为飞行员反应时间；tn为航空器响应时间，M为飞机重量，Fall为飞机所受合力；
对于椭圆形翼弦的飞机，其滚转坡度角计算公式如下：
对于后掠翼飞机，其滚转坡度角计算公式如下：
式中，p为滚转角速度，WRp为滚转阻尼力矩系数，取值为-0.5；WRc为滚转操纵力矩系数；H(t′-t″)为单位阶跃函数；t′为遭遇尾涡时间，t″为飞行员与飞机反应时间；Vz(y)为前机尾涡在后机机翼上诱导产生的上洗或下洗速度；A为展弦比；I为飞机转动惯量，通过无量纲回转半径来估算飞机，为无量纲回转半径。
2.根据权利要求1所述的一种基于多维指标的航空器遭遇尾流风险评估系统，其特征在于，步骤S5中，评估值用于评估飞机遭遇尾流风险具体为判断评估值对应的危险程度；若评估值小于2，则其危险程度为安全；若评估值大于等于2小于4，则其危险程度为较为安全；若评估值大于等于4小于6，则其危险程度为轻度危险；若评估值大于等于6小于8，则其危险程度为中度危险；若评估值大于等于8小于10，则其危险程度为强烈危险。
3.根据权利要求1-2任意一项所述的一种基于多维指标的航空器遭遇尾流风险评估系统，其特征在于，步骤S4具体包括以下步骤，
101，根据所述衡量指标构建四阶判断矩阵；
102，对所述四阶判断矩阵进行一致性检验；当一致性检验的结果时，判断矩阵通过一致性检验，则执行步骤103；
103，基于判断矩阵，采用均方根法计算指标权重。
4.根据权利要求3所述的一种基于多维指标的航空器遭遇尾流风险评估系统，其特征在于，所述一致性检验采用下述公式，
式中，CI为判断矩阵的一致性判断指标；λmax为判断矩阵的最大特征值；n为判断矩阵的阶数，ωi为各个衡量指标对应的指标权重。
5.根据权利要求3所述的一种基于多维指标的航空器遭遇尾流风险评估系统，其特征在于，所述均方根法采用下述公式，
式中，aij为判断矩阵中的i行、j列的元素。
6.根据权利要求3所述的一种基于多维指标的航空器遭遇尾流风险评估系统，其特征在于，所述评价函数采用线性加权法，所述线性加权法采用如下公式，
式中，z为评估值；C(ai)为飞机各衡量指标的评分；ωi为各个衡量指标对应的指标权重。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111154614.1A CN113837639B (zh) | 2021-09-29 | 2021-09-29 | 一种基于多维指标的航空器遭遇尾流风险评估系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111154614.1A CN113837639B (zh) | 2021-09-29 | 2021-09-29 | 一种基于多维指标的航空器遭遇尾流风险评估系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113837639A CN113837639A (zh) | 2021-12-24 |
| CN113837639B true CN113837639B (zh) | 2023-12-19 |
Family
ID=78967567
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111154614.1A Active CN113837639B (zh) | 2021-09-29 | 2021-09-29 | 一种基于多维指标的航空器遭遇尾流风险评估系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113837639B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114065407B (zh) * | 2022-01-18 | 2022-05-17 | 中国民航大学 | 一种基于机场风速多元感知的尾流间隔动态缩减方法 |
| CN114611945B (zh) * | 2022-03-15 | 2024-02-02 | 中国民用航空飞行学院 | 一种基于多维数据的飞行学员胜任力画像方法 |
| CN115526279B (zh) * | 2022-11-29 | 2023-04-28 | 中国人民解放军国防科技大学 | 基于深度学习的飞机尾流识别及强度分级方法及装置 |
| CN116186906B (zh) * | 2023-04-27 | 2023-06-30 | 中国民航大学 | 基于顶风的尾流间隔动态缩减方法及系统 |
| CN117521425B (zh) * | 2024-01-05 | 2024-03-29 | 中国民航大学 | 基于无人机与有人机混合运行的尾流间隔确定方法及系统 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103810332A (zh) * | 2014-01-27 | 2014-05-21 | 中国人民解放军空军航空大学军事仿真技术研究所 | 直升机舰面尾流影响仿真方法 |
| CN104376400A (zh) * | 2014-10-27 | 2015-02-25 | 广州市中南民航空管通信网络科技有限公司 | 基于模糊矩阵层次分析法的风险评估方法 |
| CN106650077A (zh) * | 2016-12-15 | 2017-05-10 | 中国航空工业集团公司西安飞机设计研究所 | 一种弹性飞机尾涡遭遇动响应分析方法 |
| CN109918764A (zh) * | 2019-03-01 | 2019-06-21 | 中国民用航空飞行学院 | 一种衡量飞机遭遇尾涡后的滚转风险程度的方法 |
| CN110084464A (zh) * | 2019-03-06 | 2019-08-02 | 中国矿业大学 | 一种基于物联网的智慧城市火灾危险性实时评估方法 |
| CN110705850A (zh) * | 2019-09-19 | 2020-01-17 | 安徽百诚慧通科技有限公司 | 一种基于层次分析法的道路风险评估方法 |
| CN110717675A (zh) * | 2019-10-10 | 2020-01-21 | 广西电网有限责任公司 | 一种电网建设风险评估方法 |
| CN110728449A (zh) * | 2019-10-10 | 2020-01-24 | 广西电网有限责任公司 | 一种风险管控评价方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6963291B2 (en) * | 2002-05-17 | 2005-11-08 | The Board Of Trustees Of The Leland Stanford Junior University | Dynamic wake prediction and visualization with uncertainty analysis |
| DE50312100D1 (de) * | 2003-07-25 | 2009-12-17 | Nikolai Alekseevich Baranov | Verfahren und system zur verhinderung, dass ein fl eines wirbelerzeugers eindringt |
-
2021
- 2021-09-29 CN CN202111154614.1A patent/CN113837639B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103810332A (zh) * | 2014-01-27 | 2014-05-21 | 中国人民解放军空军航空大学军事仿真技术研究所 | 直升机舰面尾流影响仿真方法 |
| CN104376400A (zh) * | 2014-10-27 | 2015-02-25 | 广州市中南民航空管通信网络科技有限公司 | 基于模糊矩阵层次分析法的风险评估方法 |
| CN106650077A (zh) * | 2016-12-15 | 2017-05-10 | 中国航空工业集团公司西安飞机设计研究所 | 一种弹性飞机尾涡遭遇动响应分析方法 |
| CN109918764A (zh) * | 2019-03-01 | 2019-06-21 | 中国民用航空飞行学院 | 一种衡量飞机遭遇尾涡后的滚转风险程度的方法 |
| CN110084464A (zh) * | 2019-03-06 | 2019-08-02 | 中国矿业大学 | 一种基于物联网的智慧城市火灾危险性实时评估方法 |
| CN110705850A (zh) * | 2019-09-19 | 2020-01-17 | 安徽百诚慧通科技有限公司 | 一种基于层次分析法的道路风险评估方法 |
| CN110717675A (zh) * | 2019-10-10 | 2020-01-21 | 广西电网有限责任公司 | 一种电网建设风险评估方法 |
| CN110728449A (zh) * | 2019-10-10 | 2020-01-24 | 广西电网有限责任公司 | 一种风险管控评价方法 |
Non-Patent Citations (7)
| Title |
|---|
| Air traffic management process quality assessment model based on improved fuzzy matter element analysis;Pan weijun 等;Journal of Physics: Conference Series;第1187卷(第5期);1-9 * |
| 基于安全阈值的尾流危险区评估;谷润平 等;安全与环境学报;第20卷(第05期);1677-1682 * |
| 梁延安.天府机场交叉跑道航空器尾流安全间隔研究.《中国优秀硕士学位论文全文数据库工程科技I辑》.2020,(第09期),1-2,27,32-38. * |
| 飞机尾流遭遇响应与危险性评估方法研究;胡皓;中国优秀硕士学位论文全文数据库 工程科技I辑(第02期);B026-103 * |
| 飞机遭遇尾涡的安全性分析;韩红蓉 等;交通运输工程学报;第12卷(第1期);45-49 * |
| 魏志强.尾涡流场及安全间隔的建模与仿真计算研究.《中国优秀硕士学位论文全文数据库工程科技II辑》.2009,(第03期),32,34. * |
| 魏志强.尾涡流场及安全间隔的建模与仿真计算研究.中国优秀硕士学位论文全文数据库工程科技II辑.2009,(第03期),31. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113837639A (zh) | 2021-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113837639B (zh) | 一种基于多维指标的航空器遭遇尾流风险评估系统 | |
| CN111382522B (zh) | 一种基于起飞滑跑数据的航空发动机安装推力评估方法 | |
| Hossain et al. | Drag analysis of an aircraft wing model with and without bird feather like winglet | |
| Lan et al. | Flight characteristics analysis based on QAR data of a jet transport during landing at a high-altitude airport | |
| US7382283B2 (en) | Method and device for detecting an overstepping of design loads of the tailplane of an aircraft | |
| US7271741B2 (en) | Method and device for detecting an overstepping of design loads of the fin of an aircraft | |
| CN109918764A (zh) | 一种衡量飞机遭遇尾涡后的滚转风险程度的方法 | |
| EP2645055A1 (en) | Method for evaluating the structural compatibility of an aircraft for use on rough runways | |
| CN112149299B (zh) | 用于飞行模拟机数据包开发的试飞任务的优化分析方法 | |
| Chang et al. | Flying qualities for a twin-jet transport in severe atmospheric turbulence | |
| De Visscher et al. | A simple wake vortex encounter severity metric | |
| CA2810801A1 (en) | Method for evaluating the structural compatibility of an aircraft for use on rough runways | |
| Robinson et al. | A flight evaluation of methods for predicting vortex wake effects on trailing aircraft | |
| Vicroy et al. | A numerical simulation study to develop an acceptable wake encounter boundary for a B737-100 airplane | |
| Dębski et al. | The Conception of The Continuous Fatigue Estimation of The Airplanes. | |
| Valenzuela et al. | Wing Design Strategies for Vehicles Designed to Operate in Ground Effect | |
| Castellani et al. | Flight flutter test of an aircraft equipped with active winglets | |
| Persoon et al. | Measurement of transonic dips in the flutter boundaries of a supercritical wing in a wind tunnel | |
| Chang et al. | Hazardous levels of commercial aircraft response to atmospheric turbulence | |
| Chang et al. | Evaluations of Flying Qualities for a New and an Aging Transport Aircraft in Severe Atmospheric Turbulence | |
| Papadakis et al. | Aerodynamic performance of a swept wing with simulated ice shapes | |
| Andrikaitis et al. | 863. Development of a finite element model of the sailplane fuselage | |
| Jiang et al. | Severity assessment of sudden plunging motion for jet transport aircraft in response to severe atmospheric turbulence | |
| Pitt et al. | Uncertainties in vibratory mode shapes and their effect on flutter speeds | |
| Stuff | The Relationship Between Near-and-Far-Field of Vortex Wakes from Aircraft with High Aspect Ratio Wings |
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 |