CN110989397B - 一种航空器失事搜寻仿真方法及系统 - Google Patents

Abstract

本发明公开了一种航空器失事搜寻仿真方法，涉及航空失事搜寻技术领域，在失事搜寻区域的确定、搜寻路径规划起到重要作用，具体方案为：包括以下步骤：S1：选择基点：视航空器为刚体，以航空器为基点；S2：定义基点的自由度，包括三个移动自由度和三个转动自由度；S3：确定航空器在大气中所受合力和合力矩情况，确定飞行状态；所述S3中飞行状态的确认包括以下步骤：M1：建立坐标系；M2：确定航空器飞行的运动参数；M3：坐标系转换：将运动参数代入坐标系中解算姿态；M4：测算控制量和状态量。通过航空器空气动力学原理，对影响航空器失事的特情进行分类和筛选，确定主要影响航空器飞行而导致失事的特情种类。

Description

一种航空器失事搜寻仿真方法及系统

技术领域

本发明涉及航空失事搜寻技术领域，更具体地说，它涉及一种航空器失事搜寻仿真方法及系统。

背景技术

失事航空器搜寻是民用航空器搜救至关重要的第一步。建立合理的失事搜寻模型，快速确定失事搜寻区域、搜寻路径规划，对失事搜寻数据进行可视化表达对民用航空器搜救具有十分重要得意义。

目前，失事搜寻区域的确定、搜寻路径的确定和失事搜寻数据可视化存在以下问题：失事搜寻区域的确定都是基于对航空器状态和天气情况等客观条件分析的基础上确定的连续空间，所确定的连续空间范围大，不利于快速搜寻到失事航空器，而计算得出的失事搜寻区域的概率分布，都是基于点、线的普通正态分布，未考虑航空器失事特情对概率分布的影响；在航空器遇到特殊情况时，飞行员主观的判断会改变航空器的飞行轨迹、路径，使得实际失事的地点具有不确定性，现有的搜寻路径规划，考虑飞行员人因影响较少，搜寻效率较低；现有的失事搜寻系统大都是失事搜寻决策系统，对航空器失事特情的研究较少，对失事搜寻产生的数据可视化表达较为简单，不利于对失事搜寻数据的理解和把握，可视化表达对搜寻产生积极影响较少。

发明内容

为解决上述技术问题，本发明的目的是一种航空器失事搜寻仿真方法，在失事搜寻区域的确定、搜寻路径规划起到重要作用。

本发明的上述技术目的是通过以下技术方案得以实现的：

一种航空器失事搜寻仿真方法，包括以下步骤：

S1：选择基点：视航空器为刚体，以航空器为基点；

S2：定义基点的自由度，包括三个移动自由度和三个转动自由度；

S3：确定航空器在大气中所受合力和合力矩情况，确定飞行状态；

所述S3中飞行状态的确认包括以下步骤：

M1：建立坐标系；

M2：确定航空器飞行的运动参数；

M3：坐标系转换：将运动参数代入坐标系中解算姿态；

M4：测算控制量和状态量；

M5：气动分析；

M6：结合基点的自由度，得到航空器飞行状态。

作为一种优选方案，S3过程中，航空器飞行状态包括爬升状态、下降状态和滚转状态。

作为一种优选方案，M1过程中坐标系包括：

地面坐标系，确定航空器相对于地面位置；

物体坐标系，描述航空器飞行姿态变化；

空气动力学坐标系，描述航空器轨迹运动；

惯性坐标系，使物体坐标系与地面坐标系重合。

作为一种优选方案，地面坐标系建立过程如下：地面坐标系O_dX_dY_dZ_d。在地面上取一点作为坐标系原点O_d，O_dZ_d轴位于地平面内，指向地球磁北方向，根据情况也可指向真北方向或者规定的任意方向；O_dY_d轴垂直于地平面指向天空；O_dX_d轴垂直于O_dY_dZ_d平面，按照左手系建立坐标系；

物体坐标系建立过程如下：物体坐标系OwXwYwZw，原点Ow位于航空器的重心。OwZw轴沿着机身轴线方向指向航空器头部；OwYw轴与OwZw轴垂直，垂直指向机体正上方；OwXw轴垂直于OwYwZw平面，按照左手系，正方向指向机体的右方；

空气动力学坐标系建立过程如下：空气动力学坐标系O_kX_kY_kZ_k又称作速度坐标系。其原点与航空器物体坐标系重合。O_kZ_k轴与航空器飞行速度矢量V重合；O_kY_k轴与O_kZ_k轴垂直，垂直指向机体上方；O_kX_k轴垂直于O_kY_kZ_k平面，按照左手系，正方向指向机体的右方；

航空器惯性坐标系建立过程如下：惯性坐标系O_gX_gY_gZ_g使得物体坐标系和地面坐标系联系起来，坐标原点与物体坐标系重合，而坐标轴的方向与地面坐标系重合。

作为一种优选方案，M2过程中，运动参数包括：气流角和姿态角，气流角包括迎角α1和侧滑角β1，姿态角包括俯仰角p、偏航角h和滚转角b。

作为一种优选方案，M3过程具体如下：

坐标系绕某一个坐标轴的旋转我们称作一个基旋转，假设坐标系C₁绕X轴旋转α角成坐标系C₂，则可得基旋转矩阵：

同理可得，绕Y轴旋转β角的基旋转矩阵：

同理可得，绕Z轴旋转γ角的基旋转矩阵：

那么绕X、Y、Z轴分别旋转得到的新坐标系的转换矩阵：

L_xyz(α,β,γ)＝L_x(α)L_y(β)L_z(γ) (4)

由式(1)至(3)可得：

惯性坐标系到物体坐标系的转换：由于惯性坐标系与物体坐标系原点重合，可利用姿态角进行旋转获得坐标系转换。由上可得：

可得惯性坐标系坐标到物体坐标系的坐标映射表达式：

同样地，空气动力学坐标系与物体坐标系可以通过迎角α和侧滑角β相联系起来。如下所示，空气动力学坐标系到物体坐标系的变换关系：

可得空气动力学坐标系坐标到物体坐标系坐标映射表达式：

作为一种优选方案，M4过程具体如下：

由上述可得，航空器控制量c和状态量s可表示如下：

其中：升降舵舵偏转角δ_e，方向舵偏转角δ_r副翼偏转角δ_a，发动机油门开合度δ_f，航空器的控制量c，航空器空速V，迎角α₁，侧滑角β₁，俯仰角p，偏航角h，滚转角b，俯仰角速度ω_p，偏航角速度ω_h，滚转角速度ω_b，地面坐标系下的坐标x_d、y_d、z_d，航空器状态量s。

作为一种优选方案，M5气动分析包括以下作用分析：重力作用、发动机推力作用、空气动力作用。

作为一种优选方案，发动力推力作用计算方法如下：

F_Xt＝C_tδ_f (11)

其中：δ_f为开合度；

空气动力作用计算方法如下：

D＝qC_DS (12)

Y＝qC_YS (13)

L＝qC_LS (14)

l＝qC_lSb (15)

n＝qC_nSb (17)

q：空气动压

其中，ρ为空气密度，V_T为航空器空速；C_D是阻力系数，C_Y是侧力系数，C_L是升力系数，C_l是滚转力矩系数，C_m是俯仰力矩系数，C_n是偏航力矩系数；S是机翼面积，b是航空器的翼展长度，

是平均气动弦长；

设航空器机翼完好率Q_jy∈[0,1]，则根据航空器空气动力学原理，航空器在大气中的受力和力矩：

D＝qC_DQ_jyS (19)

Y＝qC_YQ_jyS (20)

L＝qC_LQ_jyS (21)

l＝qC_lQ_jySb (22)

n＝qC_nQ_jySb (24)。

一种航空器失事搜寻仿真系统，仿真系统用于运行上述的仿真方法。

综上所述，本发明具有以下有益效果：

通过航空器空气动力学原理，对影响航空器失事的特情进行分类和筛选，确定主要影响航空器飞行而导致失事的特情种类。利用航空器运动学方程，建立航空器失事运动模型，计算出航空器失事的轨迹及可能的失事地点，而后通过仿真平台，对航空器无控制失事进行模拟仿真。待无控制失事过程完成后，根据失事特情的组合，通过模拟仿真获得航空器可能失事的地点，通过凸包算法和核密度算法对失事区域进行空间分析，计算出最小失事区域范围，生成失事区域热力图和等值线。在失事区域绕等值线由高到低进行连续空间搜寻。

附图说明

图1是本发明实施例的航空器失事搜寻仿真方法的飞行及失事仿真结构示意图；

图2是本发明实施例的航空器迎角、航迹角和俯仰角的结构示意图；

图3是本发明实施例的航空器侧滑角和偏航角的结构示意图；

图4是本发明实施例的航空器滚转角的结构示意图；

图5是本发明实施例的空客A320操作部件的结构示意图。

具体实施方式

以下结合附图对本发明作进一步详细说明。

在飞行仿真和失事仿真中，我们将航空器视为刚体。为了描述航空器在大气中的位置和姿态，通常定义航空器的六个自由度：相对于航空器质心的三个移动自由度、绕质心旋转的三个转动自由度。利用表征了速度矢量与航空器机体之间关系的气流角来确定航空器在大气中所受合力和合力矩等情况。

航空器在大气中运动，机体各个部件受到空气的升力和阻力作用而产生爬升、下降、滚转等动作。如图1所示，飞行及失事仿真结构。

通过对副翼、升降舵、方向舵的舵偏角产生对航空器的气动力和力矩，在大气动压、推力和重力的作用下通过运动方程解算出航空器的六自由度，并输出至应用子模块。

为了准确描述航空器的飞行状态，必须建立合适的坐标系。包括确定航空器相对于地面的位置，需建立地面坐标系；为描述航空器飞行姿态的变化，必须建立航空器物体坐标系；为描述航空器轨迹运动则需建立空气动力学坐标系。为了统一计算，本文均统一建立左手坐标系。

坐标系的分类和定义：

地面坐标系：地面坐标系O_dX_dY_dZ_d。在地面上取一点作为坐标系原点O_d，O_dZ_d轴位于地平面内，指向地球磁北方向，根据情况也可指向真北方向或者规定的任意方向；O_dY_d轴垂直于地平面指向天空；O_dX_d轴垂直于O_dY_dZ_d平面，按照左手系建立坐标系。

航空器物体坐标系：物体坐标系O_wX_wY_wZ_w，原点O_w位于航空器的重心。O_wZ_w轴沿着机身轴线方向指向航空器头部；O_wY_w轴与O_wZ_w轴垂直，垂直指向机体正上方；O_wX_w轴垂直于O_wY_wZ_w平面，按照左手系，正方向指向机体的右方。

航空器惯性坐标系：惯性坐标系O_gX_gY_gZ_g使得物体坐标系和地面坐标系联系起来，坐标原点与物体坐标系重合，而坐标轴的方向与地面坐标系重合。

空气动力学坐标系：空气动力学坐标系O_kX_kY_kZ_k又称作速度坐标系。其原点与航空器物体坐标系重合。O_kZ_k轴与航空器飞行速度矢量V重合；O_kY_k轴与O_kZ_k轴垂直，垂直指向机体上方；O_kX_k轴垂直于O_kY_kZ_k平面，按照左手系，正方向指向机体的右方。

航空器在大气中的运动参数

根据所建立的坐标系，通常定义影响航空器飞行的各个运动参数。

(1)气流角：描述了物体坐标系与航空器速度矢量之间的关系。如图2所示。

迎角α₁：又称作攻角，航空器速度矢量V在航空器对称平面内与物体坐标系O_wZ_w轴的夹角；

侧滑角β₁：航空器速度矢量V与航空器对称面之间的夹角。如图3所示。

(2)姿态角：又称欧拉角，确定了航空器相对于地面的飞行姿态。

俯仰角p(Pitch)：物体坐标系O_wZ_w轴与地面坐标系O_dX_dZ_d平面的夹角，定义航空器抬头为负；

偏航角h(Heading)：物体坐标系O_wZ_w轴在地面坐标系O_dX_dZ_d平面的投影与O_dZ_d轴的夹角，定义顺时针偏转为正；

滚转角b(Bank)：物体坐标系O_wY_w轴与地面坐标系O_dX_dY_d平面的投影与O_dX_dZ_d平面的夹角，定义航空器顺时针偏转为正。如图4所示，如果物体坐标系与地面坐标系OZ轴重合，滚转角时O_dX_d与O_wX_w的夹角。

各坐标系的转换要通过航空器各个方向受到的合力和合力矩，来解算出航空机的姿态，通常都要将各矢量转换到物体坐标系中；利用侧滑角β₁和迎角α₁来描述与气流的相对角度，以此来解算出航空器所受空气动力大小，通过空气动力学方程可以解算出航空器的飞行状态量。通过各坐标的转换可以将航空器飞行参数都统一到某坐标系下进行计算，最后转换到地面坐标系，解算出飞行状态量。

直角坐标系之间的变换，可以通过多次平移、旋转和放缩等操作获得。对于坐标原点相同的坐标系之间可以通过一系列的旋转获得变换矩阵。

通常坐标系绕某一个坐标轴的旋转我们称作一个基旋转。假设坐标系C₁绕X轴旋转α角成坐标系C₂，则可得基旋转矩阵：

同理可得，绕Y轴旋转β角的基旋转矩阵：

同理可得，绕Z轴旋转γ角的基旋转矩阵：

那么绕X、Y、Z轴分别旋转得到的新坐标系的转换矩阵：

L_xyz(α,β,γ)＝L_x(α)L_y(β)L_z(γ) (4)

由式(1)至(3)可得：

惯性坐标系到物体坐标系的转换：由于惯性坐标系与物体坐标系原点重合，可利用姿态角进行旋转获得坐标系转换。由上可得：

可得惯性坐标系坐标到物体坐标系的坐标映射表达式：

同样地，空气动力学坐标系与物体坐标系可以通过迎角α₁和侧滑角β₁相联系起来。如下所示，空气动力学坐标系到物体坐标系的变换关系：

可得空气动力学坐标系坐标到物体坐标系坐标映射表达式：

航空器的控制和状态描述

对于普通的民用航空器，通常通过控制调整升降舵、方向舵、副翼来控制航空器姿态，航空器发动机提供向前的推力。如图5所示：空客A320的操作部件的结构示意图。

升降舵：升降舵位于航空器尾翼，水平安装的升降舵可供上下偏转操作，控制航空器俯仰，舵效产生俯仰力矩。升降舵舵偏转角δ_e向下为正，产生力矩M为负。

方向舵：方向舵安装于航空器尾翼垂直位置，其左右偏转为航空器提供偏航控制，舵效产生偏航力矩。方向舵偏转角δ_r左偏为正，产生偏航力矩N为负。

副翼：副翼安装于机翼两侧翼梢后缘可活动的部分。通过控制左右副翼偏转产生滚转力矩L，副翼偏转角δ_a向下偏为正，产生滚转力矩L为负。

油门：通过控制油门大小来控制航空器推力大小。

因此，利用升降舵偏转角δ_e、方向舵偏转角δ_r、副翼偏转角δ_a和发动机油门开合度δ_f来表征航空器的控制量c。

通常，将航空器空速V，迎角α₁、侧滑角β₁、俯仰角p、偏航角h、滚转角b、俯仰角速度ω_p、偏航角速度ω_h、滚转角速度ω_b以及航空器在地面坐标系下的坐标x_d、y_d、z_d作为航空器状态量s来表征航空器在大气中的状态。

由上述可得，航空器控制量c和状态量s可表示如下：

那么，根据航空器空气动力学原理，利用一定的数学方法，由航空器的控制量可以解算出航空器的状态量。

航空器气动分析

航空器在大气中飞行主要三个外力作用：发动机推力T、空气动力A和航空器自身重力G；以及受到由在空气中运动产生的动力矩M_A和左右发动机产生的发动机力矩M_F。

(1)重力作用

重力G＝mg，m为航空器质量，g为重力加速度。重力方向为铅锤向下。航空器的质量随着燃油消耗等情况发生变化，在失事仿真过程中，燃油的变化量相对于航空器质量变化较小，又在仿真过程中将航空器设为刚体，因此，在进行结算给过程中质量m为一恒定常量。重力加速度也会随着地理位置、飞行高度而变化，但变化很小，在仿真过程中我们将重力加速度g同样设为一恒定常量，通常取为9.81。

(2)发动机推力作用

当左右发动机完好时，左右发动机推力大小相等，产生的力矩相互抵消，M_F＝0，不影响航空器偏航角速度；当航空器发生特情，由于左右发动机发生故障，而产生不同推力。发动机力矩M_F不为0，发动机力矩将影响航空器偏航角速度。

在仿真过程中，我们主要考虑利用开合度δ_f来表征发动机推力：

F_Xt＝C_tδ_f (11)

(3)空气动力作用

空气动力A在航空器物体坐标系下的分量分别为阻力D，升力L和侧力Y。阻力沿着O_wZ_w轴负向方向为正，侧力沿着O_wX_w轴正向方向为正，升力沿着O_wY_w轴正向方向为正。

M_A沿着航空器物体坐标系分量分别为滚转力矩l，俯仰力矩m，偏航力矩n。我们规定航空器向右滚转，滚转力矩为正；航空器向上抬头，其俯仰力矩为正；航空器向右偏航，其偏航力矩为正。

根据航空器空气动力学原理，我们可以计算航空器在大气中的受力和力矩^[33]：

D＝qC_DS (12)

Y＝qC_YS (13)

L＝qC_LS (14)

l＝qC_lSb (15)

n＝qC_nSb (17)

q：空气动压

其中，ρ为空气密度，V_T为航空器空速；C_D是阻力系数，C_Y是侧力系数，C_L是升力系数，C_l是滚转力矩系数，C_m是俯仰力矩系数，C_n是偏航力矩系数；S是机翼面积，b是航空器的翼展长度，

是平均气动弦长。

在考虑航空器发生特情，机翼受损，我们可以通过减少机翼的面积、翼展长度来仿真航空器机翼受损的特情。我们设航空器机翼完好率Q_jy∈[0,1]。则根据航空器空气动力学原理，航空器在大气中的受力和力矩：

D＝qC_DQ_jyS (19)

Y＝qC_YQ_jyS (20)

L＝qC_LQ_jyS (21)

l＝qC_lQ_jySb (22)

n＝qC_nQ_jySb (24)

航空器动力学方程分析

为了便于分析和解算，在失事仿真过程中，我们对航空器和仿真地球做如下假设：

(1)航空器为刚体，质量为一恒定常数；

(2)地球无公转、自转，忽略地球曲率，将地球表面看作一平面；

(3)忽略高度和地理位置对重力加速度的影响，设重力加速度为一恒定常量。

根据上述假设，刚体航空器在大气中运动满足牛顿第二定律，有如下受力分析：

为航空器所受合力矢量，m为航空器质量，

为航空器空速矢量，

为航空器所受合力矩矢量，

为航空器绕中心的动量矩。

由式(25)可得：

根据力学的基本原理，速度矢量对时间的变化率：

表示

上的单位向量，v表示航空器空速的标量，×表示叉乘，

表示航空器的总角速度矢量。

在航空器物体坐标系下可表示为：

其中，

分别表示航空器物体坐标系的O_wX_w、O_wY_w、O_wZ_w轴的单位向量；u、v、w分别表示

在O_wX_w、O_wY_w、O_wZ_w轴的速度分量；ω_b、ω_p、ω_h分别表示滚转角速度、俯仰角速度、偏航角速度。

由式(28)和(29)可得：

又：

其中，

表示u、v、w对时间的一阶微分。

由式(27)、(30)和(31)式可得：

其中，∑F_X、∑F_Y、∑F_Z是

在物体坐标系下O_wX_w、O_wY_w、O_wZ_w轴的分量。那么发动机推力向量

为：

又由矩阵变换可得航空器重力向量

空气动力向量在航空器物体坐标系下的分量：

由式(32)至(35)可得航空器在物体坐标系下的力学方程组：

又根据假设条件，动量矩随时间的变化率为^[34]：

表示

上的单位向量，H表示航空器动量矩的标量，×表示叉乘，

表示航空器的总角速度矢量。

在航空器物体坐标系下

和

可表示为：

其中，

分别表示航空器物体坐标系的O_wX_w、O_wY_w、O_wZ_w轴的单位向量；H_X、H_Y、H_Z分别表示

在O_wX_w、O_wY_w、O_wZ_w轴的速度分量；ω_b、ω_p、ω_h分别表示滚转角速度、俯仰角速度、偏航角速度。

由式(37)和(38)可得：

又：

其中，

表示H_X、H_Y、H_Z对时间的一阶微分。

根据条件假设，由力学原理可得，

I为惯性矩阵^[35]：

由于航空器几何结构在前向和尾部的质量呈对称分布，那么I_XZ＝I_YZ＝0。可得惯性矩阵：

那么，动量矩在物体坐标系下各个分量为：

可得动量矩的微分分量表达式：

由式(37)至(44)可得合力矩在物体坐标系下的分量：

当航空器发动机处于完好状态，左右发动机产生的力矩相互抵消，发动机合力矩M_F＝0。当航空器发生特情，发动机发生故障，左右发动机推力减小，产生力矩差，使航空器发生偏航。本文主要分析适合于国内航线的中小型民用航空器，因此仅考虑双发发动机航空器。设Q_lf和Q_rf分别为左右发动机完好率：

设M_lf为左发动机力矩，M_rf为右发动机力矩，右力学原理可得：

其中，F_lf为左发动机推力，F_rf为右发动机推力，l为发动机间距。规定航空器向左偏转产生负向力矩，可得发动机合力矩M_F：

由于发动机合力矩沿着物体坐标系O_wY_w轴正方向，在其他轴分量为0。那么，航空器合力矩在物体坐标系下的分量：

∑M_X、∑M_Y、∑M_Z分别等于航空器滚转力矩、俯仰力矩和偏航力矩：

航空器运动方程分析

(1)航空器姿态分析

航空器在飞行过程中姿态角p、h、b，它们与角速度的关系为：

(2)航空器位置方程分析

航空器的位置量是建立在地面坐标系下，航空器的速度是位置的一阶微分表达式，利用坐标转换可得位置矢量的一阶微分表达式：

其中，

为物体坐标系到地面坐标系转换矩阵的转置。

(3)航空器气流角和速度方程分析

航空器地速为空速和风速之和，在风速为0，地速和空速相等。通过空气动力学坐标系到物体坐标系的转换，可得空速、迎角和侧滑角的方程：

其中，L_kd为空气动力学坐标系到物体坐标系的转换矩阵。由式(53)可得：

航空器气动力系数和气动力力矩系数分析

一般航空器相对于航空器飞行轴面是基本对称的，包括外形和质量的对称。在考虑航空器基准运动，包括作对称、定常直线飞行，此时运动平面、对称面、铅垂面对称。因此可以有以下推论：

(1)纵向气动力和力矩对于横侧向运动参数在基准运动状态下的导数均为0；

(2)横侧气动力和力矩对于纵向运动参数在基准运动状态下的导数均为0。

航空器的所受空气升力L、阻力D和俯仰力矩m属于纵向力，因此升力系数C_L、阻力系数C_D、俯仰力矩系数C_m均属于纵向运动参数。侧力Y、偏航力矩n和滚转力矩l属于横侧力，因此侧立系数C_Y、偏航力矩系数C_n和滚转力矩系数C_l均属于横侧向运动参数。

因此，利用升降舵偏转角δ_e、方向舵偏转角δ_r、副翼偏转角δ_a和发动机油门开合度δ_f来表征航空器的控制量c。

航空器在大气中飞行，由于对流气流流过航空器上下空气流速不同，而产生升力。通常升力系数与迎角、马赫数、俯仰角速度、升降舵、方向舵和副翼转动等有关。在迎角、侧滑角和马赫数一定范围内，我们可以采取简单的线性模型：

C_L＝C_L0+C_Lαα+C_Lδαδ_α+C_Lδeδ_e (57)

其中，C_L0为航空器基本升力系数；C_Lα为升力系数对迎角α₁的导数；C_Lδα为升力系数对副翼偏转角δ_a的导数；C_Lδe为升力系数对升降舵偏转角δ_e的导数。

当气流流过航空器表面时，由于空气具有黏性，空气微团与航空器表面发生摩擦，阻滞了气流的流动而产生摩擦阻力；在飞行过程中，由于航空器前面压强大，后面压强小，前后形成了巨大的压强差而产生压差阻力；又由于下洗角的存在，升力产生与航空器前进方向相反的水平分力阻止航空器前进而产生诱导阻力；航空器个部件之间，由于气流的相互干扰还会产生部分额外的阻力。通常阻力系数与迎角、马赫数、侧滑角、副翼、升降舵、方向舵、地面效应和起落架等有关，我们在考虑航空器失事过程中除去地面效应和起落架的影响等情况，可以采取简单的线性模型来表示：

C_D＝C_D0+C_Dαα+C_Dδαδ_α+C_Dδeδ_e+C_Dδrδ_r (58)

其中，C_D0为航空器零升阻力系数；C_Dα为阻力系数对迎角α₁的导数；C_Dδα为阻力系数对副翼偏转角δ_a的导数；C_Dδe为阻力系数对升降舵偏转角δ_e的导数；C_Dδr为为阻力系数对方向舵偏转角δ_r的导数。

航空器在不对称气流的作用下会产生侧力。其侧力系数C_Y通常与侧滑角、副翼、升降舵、方向舵偏角等有关，通常表示为：

C_Y＝C_Y0+C_Yββ+C_Yδαδ_α+C_Yδeδ_e+C_Yδrδ_r (59)

其中，C_Y0为航空器基本侧力系数；C_Yβ为侧力系数对侧滑角β的导数；C_Yδα为侧力系数对副翼偏转角δ_a的导数；C_Yδe为侧力系数对升降舵偏转角δ_e的导数；C_Yδr为侧力系数对方向舵偏转角δ_r的导数。

滚转力矩系数C_l通常与迎角、侧滑角、副翼和方向舵的偏转等有关，通常表示为：

其中，C_l0为航空器基本滚转力矩系数；C_lβ为滚转力矩系数对侧滑角β的导数；C_lδe为滚转力矩系数对升降舵偏转角δ_e的导数；C_lδr为滚转力矩系数对方向舵偏转角δ_r的导数；C_lωh为滚转力矩系数对偏航角速度的导数；C_lωb为滚转力矩系数对滚转角速度的导数；b为航空器翼展长度；V为航空器空速；ω_h、ω_b分别为航空器偏航角和滚转角速度。

俯仰力矩系数C_m通常与迎角、马赫数、升降舵偏转和俯仰角度等有关，可表示为：

其中，C_m0为航空器基本俯仰力矩系数；C_mα为俯仰力矩系数对侧滑角α的导数；C_mδα为俯仰力矩系数对副翼偏转角δ_a的导数；C_mδe为俯仰力矩系数对升降舵偏转角δ_e的导数；C_mδr为俯仰力矩系数对方向舵偏转角δ_r的导数；C_mωp为俯仰力矩系数对偏航角速度的导数；

为航空器平均气动弦长；V为航空器空速；ω_p为航空器俯仰角速度。

偏航力矩系数C_n通常与侧滑角、副翼偏转、方向舵偏转、偏航角速度和滚转角速度等因素有关，可是表示为：

C_n＝C_n0+C_nββ+C_nδαδ_α+C_nδeδ_e+C_nδrδ_r (62)

其中，C_n0为航空器基本偏航力矩系数；C_nβ为偏航力矩系数对侧滑角β的导数；C_nδα为偏航力矩系数对副翼偏转角δ_a的导数；C_nδe为偏航力矩系数对升降舵偏转角δ_e的导数；C_nδr为偏航力矩系数对方向舵偏转角δ_r的导数。

当航空器发生特情，升降舵、方向舵、副翼发生故障，航空器的空气动力系数将发生改变，影响航空器的受力和力矩，从而改变航空器的姿态和速度。我们在这里分别设升降舵、方向舵、副翼的完好率为Q_δe、Q_δr、Q_δa。则航空器的空气动力系数表示如下：

C_L＝C_L0+C_Lαα+C_LδαQ_δαδ_α+C_LδeQ_δeδ_e (63)

C_D＝C_D0+C_Dαα+C_DδαQ_δαδ_α+C_DδeQ_δeδ_e+C_DδrQ_δrδ_r (64)

C_Y＝C_Y0+C_Yββ+C_YδαQ_δαδ_α+C_YδeQ_δeδ_e+C_YδrQ_δrδ_r (65)

C_n＝C_n0+C_nββ+C_nδαQ_δαδ_α+C_nδeQ_δeδ_e+C_nδrQ_δrδ_r (68)

采用最小二乘法对QAR数据进行辨识得到A320气动导数如下：

表1 升力系数导数

表2 阻力系数导数

表3 侧力系数导数

表4 滚转力矩系数导数

表5 俯仰力矩系数导数

表6 偏航力矩系数导数

本具体实施例仅仅是对本发明的解释，其并不是对本发明的限制，本领域技术人员在阅读完本说明书后可以根据需要对本实施例做出没有创造性贡献的修改，但只要在本发明的权利要求范围内都受到专利法的保护。

Claims (9)

1.一种航空器失事搜寻仿真方法，其特征在于，包括以下步骤：

S1：选择基点：视航空器为刚体，以航空器为基点；

S2：定义基点的自由度，包括三个移动自由度和三个转动自由度；

S3：确定航空器在大气中所受合力和合力矩情况，确定飞行状态；

所述S3中飞行状态的确认包括以下步骤：

M1：建立坐标系；

M2：确定航空器飞行的运动参数；

M3：坐标系转换：将运动参数代入坐标系中解算姿态；

M4：测算控制量和状态量；

M5：气动分析；

M6：结合基点的自由度，得到航空器飞行状态；

其中，所述M3过程具体如下：

坐标系绕某一个坐标轴的旋转我们称作一个基旋转，假设坐标系C₁绕X轴旋转α角成坐标系C₂，则可得基旋转矩阵：

同理可得，绕Y轴旋转β角的基旋转矩阵：

同理可得，绕Z轴旋转γ角的基旋转矩阵：

那么绕X、Y、Z轴分别旋转得到的新坐标系的转换矩阵：

L_xyz(α,β,γ)＝L_x(α)L_y(β)L_z(γ) (4)

由式(1)至(3)可得：

惯性坐标系到物体坐标系的转换：由于惯性坐标系与物体坐标系原点重合，可利用姿态角进行旋转获得坐标系转换，由上可得：

可得惯性坐标系坐标到物体坐标系的坐标映射表达式：

同样地，空气动力学坐标系与物体坐标系可以通过迎角α1和侧滑角β1相联系起来，如下所示，空气动力学坐标系到物体坐标系的变换关系：

可得空气动力学坐标系坐标到物体坐标系坐标映射表达式：

2.根据权利要求1所述的航空器失事搜寻仿真方法，其特征在于，所述S3过程中，航空器飞行状态包括爬升状态、下降状态和滚转状态。

3.根据权利要求1所述的航空器失事搜寻仿真方法，其特征在于，所述M1过程中坐标系包括：

地面坐标系，确定航空器相对于地面位置；

物体坐标系，描述航空器飞行姿态变化；

空气动力学坐标系，描述航空器轨迹运动；

惯性坐标系，使物体坐标系与地面坐标系重合。

4.根据权利要求3所述的航空器失事搜寻仿真方法，其特征在于，地面坐标系建立过程如下：地面坐标系OdXdYdZd，在地面上取一点作为坐标系原点Od，OdZd轴位于地平面内，指向地球磁北方向，根据情况也可指向真北方向或者规定的任意方向；OdYd轴垂直于地平面指向天空；OdXd轴垂直于OdYdZd平面，按照左手系建立坐标系；

物体坐标系建立过程如下：物体坐标系OwXwYwZw，原点Ow位于航空器的重心，OwZw轴沿着机身轴线方向指向航空器头部；OwYw轴与OwZw轴垂直，垂直指向机体正上方；OwXw轴垂直于OwYwZw平面，按照左手系，正方向指向机体的右方；

空气动力学坐标系建立过程如下：空气动力学坐标系OkXkYkZk又称作速度坐标系，其原点与航空器物体坐标系重合，OkZk轴与航空器飞行速度矢量V重合；OkYk轴与OkZk轴垂直，垂直指向机体上方；OkXk轴垂直于OkYkZk平面，按照左手系，正方向指向机体的右方；

航空器惯性坐标系建立过程如下：惯性坐标系OgXgYgZg使得物体坐标系和地面坐标系联系起来，坐标原点与物体坐标系重合，而坐标轴的方向与地面坐标系重合。

5.根据权利要求1所述的航空器失事搜寻仿真方法，其特征在于，所述M2过程中，运动参数包括：气流角和姿态角，气流角包括迎角α1和侧滑角β1，姿态角包括俯仰角p、偏航角h和滚转角b。

6.根据权利要求1所述的航空器失事搜寻仿真方法，其特征在于，所述M4过程具体如下：

由上述可得，航空器控制量c和状态量s可表示如下：

其中：升降舵舵偏转角δ_e，方向舵偏转角δ_r副翼偏转角δ_a，发动机油门开合度δ_f，航空器的控制量c，航空器空速V，迎角α1，侧滑角β1，俯仰角p，偏航角h，滚转角b，俯仰角速度ω_p，偏航角速度ω_h，滚转角速度ω_b，地面坐标系下的坐标x_d、y_d、z_d，航空器状态量s。

7.根据权利要求1所述的航空器失事搜寻仿真方法，其特征在于，所述M5气动分析包括以下作用分析：重力作用、发动机推力作用、空气动力作用。

8.根据权利要求7所述的航空器失事搜寻仿真方法，其特征在于，所述发动力推力作用计算方法如下：

F_Xt＝C_tδ_f (11)

其中：δ_f为开合度；

所述空气动力作用计算方法如下：

D＝qC_DS (12)

Y＝qC_YS (13)

L＝qC_LS (14)

l＝qC_lSb (15)

n＝qC_nSb (17)

q：空气动压

其中，ρ为空气密度，V_T为航空器空速；C_D是阻力系数，C_Y是侧力系数，C_L是升力系数，C_l是滚转力矩系数，C_m是俯仰力矩系数，C_n是偏航力矩系数；S是机翼面积，b是航空器的翼展长度，

是平均气动弦长；

设航空器机翼完好率Q_jy∈[0,1]，则根据航空器空气动力学原理，航空器在大气中的受力和力矩：

D＝qC_DQ_jyS (19)

Y＝qC_YQ_jyS (20)

L＝qC_LQ_jyS (21)

l＝qC_lQ_jySb (22)

n＝qC_nQ_jySb (24)。

9.一种航空器失事搜寻仿真系统，其特征在于，所述仿真系统用于运行权利要求1至8任一所述的仿真方法。

Images