一种基于多维指标的航空器遭遇尾流风险评估系统

摘要

本发明公开了一种基于多维指标的航空器遭遇尾流风险评估系统。本发明系统用于执行以下步骤：S1，获取飞机遭遇尾涡风险的衡量指标的值，衡量指标包括滚动力矩系数、过载增量、高度下降率和滚转坡度角；S2，建立各个衡量指标的危险等级评估表，并基于衡量指标的值，获取飞机各个衡量指标的危险等级；S3，建立评语集，并基于各个衡量指标的危险等级和评语集，获取各个衡量指标的评分；S4，采用层次分析法对衡量指标进行分析，得到飞机各个衡量指标对应的指标权重；S5，构建飞机的评价函数，并带入飞机各个衡量指标的危险等级评分和指标权重，计算评估值，评估值用于评估飞机遭遇尾流风险。具有能够量化评估后机遭遇尾涡危险风险的技术效果。

说明书

技术领域

本发明涉及尾涡风险评价领域，特别涉及一种基于多维指标的航空器遭遇尾流风险评估系统。

背景技术

尾涡是飞机产生升力的一种副产品，表现为在飞机后方卷起形成一对向外向下旋转的涡流，并在飞机后方持续几分钟时间。尾涡的强度由产生尾涡的飞机重量、飞行速度机翼形状和翼展长度所决定，其中最主要的是飞机的重量。尾涡强度随飞机重量、载荷因数的增加和飞行速度的减小而增大。参照图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，则其危险程度为强烈危险。

可选地，在获取飞机的滚动力矩系数时，将飞机机翼分割为条带状，形成多个微元机翼，再计算滚动力矩系数；所述计算飞机的滚动力矩系数采用下述公式计算，

L

式中，RMC为滚转力矩系数，V为空气的来流速度；S为机翼面积，B为前机翼展；F为机翼上其中一个微元机翼升力的变化量，L

可选地，所述过载增量采用下述公式计算，

式中，F

F

式中，ΔF

可选地，所述高度下降率采用下述公式计算，

其中，G

式中，t

可选地，计算滚转坡度角时，采用滚转角速度计算公式计算滚转角速度，再对所述滚转角速度进行积分，得到滚转坡度角，所述滚转角速度计算公式为：

式中，

可选地，步骤S4具体包括以下步骤，

101，根据所述衡量指标构建四阶判断矩阵；

102，对所述四阶判断矩阵进行一致性检验；当一致性检验的结果时，判断矩阵通过一致性检验，则执行步骤103；

103，基于判断矩阵，采用均方根法计算指标权重。

可选地，所述一致性检验采用下述公式，

式中，CI为判断矩阵的一致性判断指标；λ

可选地，所述均方根法采用下述公式，

式中，a

可选地，所述评价函数采用线性加权法，所述线性加权法采用如下公式，

式中，z为评估值；C(a

与现有技术相比，本发明的有益效果：通过选取滚转力矩系数、过载增量、高度下降率和滚转角速度四个衡量指标来评价飞机遭遇尾涡风险。首先利用层次分析法建立航空器遭遇尾流风险评价指标体系，获取衡量指标的指标权重。获取衡量指标的值，并基于衡量指标和指标权重构建评价函数，对不同翼型的航空器遭遇前机尾流的风险进行评估。该方法可以实时动态地监视不同航段上航空器遭遇尾流风险态势，并对避免尾流危险遭遇，提高进近阶段运行安全，协助管制员实时决策提供了参考。

附图说明：

图1为本发明中飞机尾涡遭遇示意图；

图2为本发明的方法流程图；

图3为本发明中条状带模型示意图；

图4为本发明中条状带模型部分示意图。

具体实施方式

下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例，凡基于本发明内容所实现的技术均属于本发明的范围。

实施例1

一种基于多维指标的航空器遭遇尾流风险评估系统，参照图1，系统用于执行下述步骤，

S1，获取飞机遭遇尾涡风险的衡量指标的值，衡量指标包括滚动力矩系数、过载增量、高度下降率和滚转坡度角。

求取衡量指标的值之前，需要先判断飞机类型，飞机类型包括圆翼飞机和后掠翼飞机；飞机的机型不同，则在评估飞机遭遇尾流风险时的评估过程也不一致，在本实施例中，主要是体现在衡量指标的数据不同。

在获取飞机的滚动力矩系数时，将飞机机翼分割为条带状，形成多个微元机翼，再计算滚动力矩系数；计算飞机的滚动力矩系数采用下述公式，

L

其中，RMC为滚转力矩系数，V为空气的来流速度，与飞机的飞行速度一致；S为机翼面积，B为前机翼展；F为机翼上某个微元机翼升力的变化量，L

在计算滚转力矩系数时，由于机翼各处受力不均，因此可以采用微元思想，将机翼分割成无数条状带，将每条条带状受力当成一个整体，之后按机翼展展长方向进行积分，即可以得到机翼整体的升力变化量。条状带模型图如图3和图4所示，每条条带状上升力变化量乘以其到相应机翼中心的距离，即为后机机翼该条状带区域在前机尾涡诱导作用下的诱导滚转力矩，之后进行累加，即为机翼整体的滚转力矩。

由条状带法思想，我们计算一个条带上由尾涡引起的升力变化量，进而得到此条带上的诱导力矩。

式中，dF(y)表示为局部的升力变化量；ΔL

针对两种飞机类型；椭圆形翼弦弦长cl(y)采用下述公式计算：

后掠翼飞机翼弦弦长cl(y)采用下述公式计算：

式中，c

因此，椭圆形翼弦飞机的滚转力矩系数为：

后掠翼飞机飞机的滚转力矩系数为：

尾涡的强度大小用环量加以衡量，其中，

r

式中，Γ

尾涡耗散速度的快慢与湍流耗散率有关，t

又有，

当ε

当0.0235≥ε

当0.0121≥ε

t

当0.001≥ε

t

尾涡在经过卷起，形成两个方向涡流后，出现线性不稳定，涡量开始径向扩散，之后不稳定性导致两个涡流之间的强烈相互作用，最终崩溃耗散成大气中湍流。把耗散过程分为两个阶段，即扩散阶段和快速衰减阶段；

扩散阶段耗散公式如下：

快速快速衰减阶段公式如下：

且前机产生的尾涡会对后机产生一个诱导速度，根据Burnham-Hallock模型，可得诱导速度采用下式计算：

式中，Γ为尾涡经过一定时间耗散之后的环量，r为该点到涡心处的距离，v为该点的垂直诱导速度。

前机产生的尾涡在后机机翼上每一点的切向速度都不相同，因此每一点的升力变化量也不相同。当后机从前机尾涡场产生的两涡之间任意一点位置进入时，由于左右机翼受力不均，会以两翼中点为中心产生一个诱导滚转力矩，使飞机发生滚转，严重影响飞行安全。

在ICAO标准中，采用滚转力矩系数这一无量纲化参数去比较具有不同物理和空气动力学特性的飞机之间遭遇尾涡的严重性程度,当滚转力矩系数在起飞与进近阶段小于0.04时、平飞阶段小于0.06时，可认为飞机安全。

在求过载增量时，采用下述公式，

式中，F

F

式中，ΔF

式中，ρ为大气密度，V为来流速度，即飞机真空速；cl(y)为升力线斜率；c(y)为飞机翼展长坐标处的翼弦长；V

式中，s(y)为尾翼弦长；C为尾翼展长；

式中，b(y)为发动机翼弦长；D为发动机展长；

ΔF

式中，N为法向力，A为轴向力；α为迎角。

在求高度下降率时，采用下述公式计算；

其中，G

式中，t

计算滚转坡度角时，首先采用下述公式计算滚转角速度，

式中，

I采用下式计算：

式中，

可以得出，椭圆形翼弦的飞机，其滚转角速度计算公式如下：

后掠翼飞机，其滚转角速度计算公式如下：

再对滚转角速度进行积分，得到滚转坡度角。

对于椭圆形翼弦的飞机，其滚转坡度角计算公式如下：

对于椭圆形翼弦的飞机，其滚转坡度角计算公式如下：

S2，建立各个衡量指标的危险等级评估表，并基于衡量指标的值，获取飞机各个衡量指标的危险等级。

其中，滚动力矩系数、过载增量、高度下降率和滚转角速度的危险等级评估表如下所示，危险等级评估表基于衡量指标的干扰强度划分各个衡量指标的危险等级，危险等级包括1、2、3、4、5五个等级。

其中，滚动力矩系数干扰强度对应的危险等级评估如下表所示：

过载增量颠簸强度对应的危险等级评估如下表所示：

高度下降率干扰强度对应的危险等级评估如下表所示：

滚转坡度角干扰强度对应的危险等级评估如下表所示：

S3，建立评语集，并基于各个衡量指标的危险等级和评语集，获取各个衡量指标的评分；评语集包括1、2、3、4、5五个危险等级，分别对应2、4、6、8、10五个评分。

如下表所示：

S4，采用层次分析法对衡量指标进行分析，得到飞机各个衡量指标对应的指标权重。

层次分析法包括下述步骤，

101，根据衡量指标构建四阶判断矩阵；

对航空器遭遇尾流风险指标体系中相关因素指标的相对重要性进行两两比较判断，构造判断矩阵，采用T.L.Saaty建议的1～9标度法。

标度表如下：

准则层包括滚转力矩系数(a

其中，a

102，对四阶判断矩阵进行一致性检验；当一致性检验的结果时，判断矩阵通过一致性检验，则执行步骤103；

在建立判断矩阵时，由于客观事物具有复杂性，主题认识具有主观性和多样性，判断矩阵不可能具有完全一致性，因此需要进行判断矩阵的一致性检验。

一致性检验采用下述公式，

式中，CI为判断矩阵的一致性判断指标；λ

且当判断矩阵的阶数较大时，则引入随机一致性指标RI(Random Index)进行修正，经修正的一致性指标用CR表示。

式中，RI取值参照下表，n为判断矩阵的阶数。

当求出的结果CR＜0.1时，判断矩阵通过一致性检验。

103，基于判断矩阵，采用均方根法计算指标权重。

均方根法计算指标权重采用如下公式，

式中，a

得出四个评价指标所对应的指标权重ω

S5，将衡量指标的评分和各个衡量指标对应的指标权重代入飞机的评价函数，计算出评估值，评估值用于评估飞机遭遇尾流风险。

本实施例中，采用线性加权法构建评价函数，线性加权法如以下公式，

式中，z为评估值；C(a

分别代入两种类型飞机衡量指标的值，得到两种类型飞机的评价函数。

后机与前机的相对位置可由三维坐标系(x,y,z)来表示，后机在前机尾涡影响范围内的任意位置处风险均要考虑滚转力矩系数(RMC)、过载增量(|Δn|)、高度下降率(Gr)和滚转角加速度(Φ)这四种衡量航空器危险程度的指标。

根据民用航空器分为圆翼和后掠翼飞机，可设计评价函数为：

z＝ω

求取评估值后，将评估值用于评估飞机遭遇尾流风险，判断评估值对应的危险程度；

若评估值小于2，则其危险程度为安全；

若评估值大于等于2小于4，则其危险程度为较为安全；

若评估值大于等于4小于6，则其危险程度为轻度危险；

若评估值大于等于6小于8，则其危险程度为中度危险；

若评估值大于等于8小于10，则其危险程度为强烈危险。

如下表所示，

通过查阅上表，工作人员能够根据求取的评估值来实时动态地监视不同航段上航空器遭遇尾流风险态势，并对避免尾流危险遭遇，提高进近阶段运行安全，为协助管制员实时决策提供参考。

两种类型飞机中，四个衡量指标的计算公式不同，则其评价函数也不一致，分别构建评价函数后；能够使用评价函数来将滚转力矩系数、过载增量、高度下降率和滚转坡度角四个衡量指标对航空器遭遇尾流风险进行综合评价，对不同翼型的航空器遭遇前机尾流的风险进行评估。

实施例2

在本实施例中，选取国产ARJ21(中型机)飞机跟随A330-200(重型机)飞机对实施例1中所述的方法进行验证，ICAO规定近地阶段中型机跟随重型机时尾流安全间隔为9.3km。本实施例选取进近阶段时后机穿越前机尾涡流场的三个位置：

1)飞机纵向穿越单涡涡核中心，图1中B位置；

2)飞机纵向穿越距两涡核连线中点处左右各10m处位置；

3)穿越两涡涡核连线中点处。

选取以上三个位置来判断飞机此时遭遇尾流的风险程度，判断ICAO间隔标准是否安全。

其中，ARJ21参数与A330-200参数如下表所示：

尾流及环境参数如下表所示：.

将上述参数带入到实施例1中得到的评价参数公式中，能够计算得到如下表所示的计算结果：

再带入到危险等级划分表中，得到衡量指标的得分结果如下表所示：

再根据，层次分析法得出ω

因此对于纵向穿越单涡涡核中心这种情况：

可以得到，评价函数为：

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处；穿越单涡涡核中心这种情况危险等级高于以上两种情况。根据上述结论，本评估系统能够为管制员实时决策提供参考。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。