飞机结冰是飞行实践中广泛存在的一种物理现象[1]，通常会破坏飞机的动力学特性，影响飞机的操稳性特性，降低飞行品质，严重威胁飞机飞行安全[2]。自20世纪20年代以来，飞机结冰一直是威胁航空安全的重要问题[3]，美国Safety Advisor对1990—2000年间因气象原因造成的飞行事故进行统计，共有3 230起，其中结冰事故占12%[4]，随后在2003—2008年间，又发生了380起因为结冰导致的飞行事故。
针对结冰对飞行安全造成的威胁，国外学者从20世纪30年代开始对结冰现象进行研究。Lynch和Khodadoust深入系统地研究了不同结冰位置、结冰形状对飞机升、阻力系数的影响，并量化分析了结冰引起的气动特性变化[5]。Bragg和Basar在国际结冰会议上提出了进行结冰程度量化的模型[6]，为结冰的量化分析奠定了基础。国内以研究飞机结冰、防除冰机理[7], 操纵保护，风险量化、评估[8]为重点。
为了降低飞机在结冰条件下的飞行风险，必须建立飞机的飞行安全包线，而传统飞行安全包线保护通常只用于防止飞行员操纵飞机进入超出飞机结构和气动极限的状态。1994年，美国鹰航公司一架ATR72-212飞机坠毁于印第安纳州，当时飞机的飞行迎角仅为5°，远低于飞控系统提供的迎角限制值18.1°，可是在机翼结冰的影响下，飞机仍然因失速而导致坠毁[9]。因此，必须获得飞机在复杂结冰环境下的飞行安全包线。目前，国内外确定飞行安全包线的方法很多，其中获得广大研究者公认的方法之一就是可达性分析[10]。可达性分析是指针对微分动力系统，是否可以通过有效控制策略、在特定的时间范围内使初始集范围内的状态达到指定目标集合的过程。该方法应用范围较广，在航空领域安全分析中表现突出。文献[11]利用可达集方法计算了飞机着陆阶段的飞行安全包线，并分析了襟翼切换对于飞机着陆时飞行安全包线的影响，文献[12]以最大可控不变集作为飞机的安全集，并基于安全集研究了机动边界的保护控制，文献[13]基于可达集方法求解飞行器筋斗动作各阶段的操纵边界，为飞行员完成标准动作提供决策支撑。由于可达集理论的实用性，一些学者将可达集理论应于对结冰飞行安全的研究，文献[14-15]将可达集用于结冰飞机着陆阶段的飞行安全研究，并分别得出不同结冰程度下的飞行风险与操纵策略，但是都是在结冰完全已知的情况下得出的结论，具有一定的局限性。
为了分析复杂结冰环境对飞机飞行安全包线的影响，首先利用时间尺度分离的方法将飞机运动方程解耦成两个系统，为方便研究选取快子系统为研究对象，以快子系统输出为虚拟输入，建立考虑飞机横向运动(β、φ)的纵向运动方程。然后根据结冰对飞机气动特性的影响，建立考虑结冰位置、冰型及分布等不确定性的复杂结冰环境下的鲁棒气动导数模型。运用可达集理论，得到飞机在不同结冰程度下的飞行安全包线，对结冰条件下飞机滚转运动进行分析。最后进一步分析考虑不同结冰程度下结冰位置、冰型及分布等不确定性结冰情况时的鲁棒飞行安全包线，由此得出结冰不确定性对飞机飞行安全带来的影响。
1 结冰飞机动力学建模飞机本体运动方程为高维非线性，直接对本体进行研究过于复杂，本文采用时间尺度分离的方法，将飞机本体运动解耦成两个子系统[16-17]，如图 1所示。第1个子系统是以飞机实际舵面控制及推力[δa δe δr T]T为输入，姿态角向量[φ α β]T为输出的快子系统，另一个系统是以姿态角[φ α β]T为虚拟输入、[χ γ V]T为输出的慢子系统。针对不同的控制子系统，飞机的控制方式不同，本文主要研究第2个控制子系统。
|
| 图 1 飞机动力学时间尺度分离 Fig. 1 Time scale separation of aircraft dynamics |
为简化分析过程，以RCAM(Research Civil Aircraft Model)简化模型为研究对象，飞机所有气动参数基于文献[18]。
1.1 飞机非线性动力学模型将运动飞机视为一个质点，其受力如图 2所示。图中：Xb、Yb分别表示机体坐标系下的纵轴及竖轴；Xω、Yω分别表示速度坐标系下的阻力轴及升力轴；α为迎角；θ为俯仰角。
|
| 图 2 飞机受力质点图 Fig. 2 Point mass force diagram for aircraft |
建立飞机非线性动力学方程:
| $ \left\{ \begin{array}{l} {F_X} - mg\sin \gamma = m\dot V\\ {F_Z}\cos \varphi + {F_Y}\sin \varphi + mg\cos \gamma = - mV\dot \gamma \\ - {F_Z}\sin \varphi + {F_Y}\cos \varphi = mV\dot \chi \cos \dot \gamma \end{array} \right. $ | （1） |
式中：m为飞机质量；V为飞行速度；γ为航迹倾角；φ为滚转角；χ为航迹方位角；FX、FY、FZ分别表示沿X、Y、Z轴方向的合力，可表示为
| $ \left\{ \begin{array}{l} {F_X} = T\cos \beta \cos \alpha - D\left( {V, \alpha } \right) \approx T - D\left( {V, \alpha } \right)\\ {F_Y} = - T\sin \alpha - L\left( {V, \alpha } \right) \approx - L\left( {V, \alpha } \right)\\ {F_Z} = - T\sin \beta \cos \alpha + Y\left( {V, \beta } \right) \approx Y\left( {V, \beta } \right) \end{array} \right. $ | （2） |
式中：升力L、阻力D和侧力Y可表示为
| $ \left\{ \begin{array}{l} L\left( {V, \alpha } \right) = 0.5\rho {V^2}S\left( {{C_{{L_0}}} + {C_{{D_\alpha }}}\alpha } \right)\\ D\left( {V, \alpha } \right) = 0.5\rho {V^2}S\left( {{C_{{D_0}}} + {C_{{D_\alpha }}}\alpha + {C_{{D_{{\alpha ^2}}}}}{\alpha ^2}} \right)\\ Y\left( {V, \beta } \right) = 0.5\rho {V^2}S\left( {{C_Y}\beta } \right) \end{array} \right. $ | （3） |
其中：ρ为空气密度；S为机翼参考面积；CL、CD、CY分别为升力、阻力和侧力系数。
选取纵向通道状态为x=[V γ]T，控制输入为u=[T φ α β]T，因此，所研究的飞机动力学方程可表示为
| $ \begin{array}{l} \left[ {\begin{array}{*{20}{c}} {\dot V}\\ {\dot \gamma } \end{array}} \right] = \left[ {\begin{array}{*{20}{l}} { - g\sin \gamma }\\ { - \frac{g}{V}\cos \gamma } \end{array}} \right] + \\ \;\;\;\;\;\;\;\left[ {\begin{array}{*{20}{c}} {\cos \alpha \cos \beta }\\ {\left( {\cos \varphi \sin \alpha \cos \beta - \sin \varphi \sin \beta } \right)\frac{1}{V}} \end{array}} \right]\frac{T}{m} + \\ \;\;\;\;\;\;\;\left[ {\begin{array}{*{20}{c}} { - \frac{{\rho S}}{{2m}}{V^2}\left( {{C_{{D_0}}} + {C_{{D_\alpha }}}\alpha + {C_{{D_{{\alpha ^2}}}}}{\alpha ^2}} \right)}\\ {\frac{{\rho S}}{{2m}}V\left( {{C_{{L_0}}} + {C_{{L_\alpha }}}\alpha } \right)\cos \varphi } \end{array}} \right] + \\ \;\;\;\;\;\;\;\left[ {\begin{array}{*{20}{c}} 0\\ { - \frac{{\rho S}}{{2m}}V{C_Y}\beta \sin \varphi } \end{array}} \right] \end{array} $ | （4） |
飞机结冰是指飞机机体表面的某些部位聚集冰层的现象，主要会导致飞机空气动力学特性变差，升力系数减小、阻力系数增大，同时使得飞机最大飞行迎角减小，如图 3所示。
|
| 图 3 飞机结冰前后升力系数变化示意图 Fig. 3 Diagram of lift coefficient of aircraft changesin clean and iced |
根据Bragg和Basar提出的量化结冰模型，构建结冰程度对飞机气动参数影响的模型[6]：
| $ {C_{\left( {\rm{A}} \right){\rm{iced}}}} = \left( {1 + \eta {k_{{C_{\left( {\rm{A}} \right)}}}}} \right){C_{\left( {\rm{A}} \right)}} $ | （5） |
式中：C(A)和C(A)iced分别为飞机结冰前后的某个气动导数值；kC(A)为飞机结冰因子，反映气动导数受结冰影响程度，通常由飞机本身结构决定, 对于特定飞机为一个常值；η为表征飞机结冰程度的参数，未结冰时η=0，随着结冰程度的增加，η的值不断增大。
Bragg的结冰模型虽然能够定量反映飞机结冰的影响，但是一方面，飞机结冰程度难以准确测量，另外，结冰是一个非常复杂的过程，不同结冰程度、不同冰型、不同结冰位置等，都会对飞机的气动特性造成不同影响，并不能简单用η来描述。本文针对结冰难以准确量化问题，根据结冰对飞机气动特性的影响，提出鲁棒气动导数模型, 可表示为
| $ \left\{ \begin{array}{l} {{\hat C}_L} = {{\hat C}_{L0}} + {{\hat s}_{L0}} \cdot \Delta \\ {{\hat C}_D} = {{\hat C}_{D0}} + {{\hat s}_{D0}} \cdot \Delta \\ {{\hat C}_Y} = {{\hat C}_{Y0}} + {{\hat s}_{Y0}} \cdot \Delta \end{array} \right. $ | （6） |
式中：
基于可达集理论，利用可达性分析获得的反向可达集作为结冰飞机最大飞行安全包线，一旦飞行状态超出反向可达集范围，飞机就有失控的危险。
2.1 可达集理论飞机的动态特性可以表示为一个连续时间控制：
| $ \mathit{\boldsymbol{\dot x}} = f\left( {\mathit{\boldsymbol{x}}, t, \mathit{\boldsymbol{u}}} \right) $ | （7） |
式中：x∈Rn表示n维状态空间；u∈U⊆Rm表示m维控制量且f(·, ·):Rn×U→Rn为有界、连续的Lipschitz函数。
定义状态轨迹为
| $ \mathit{\boldsymbol{\xi }}\left( {{\mathit{\boldsymbol{x}}_0}, {t_0}, \mathit{\boldsymbol{u}}\left( \cdot \right), t} \right):t \to \mathit{\boldsymbol{x}} \in {\mathit{\boldsymbol{R}}^n} $ | （8） |
式中：x0为初始时刻的状态；u(·)为此时的控制输入。则反向可达集RBK, t0(t)可表示为
| $ \begin{array}{l} R_{K, {t_0}}^{\rm{B}}\left( t \right) = \left\{ {\mathit{\boldsymbol{x}} \in {\mathit{\boldsymbol{R}}^n}:\exists \mathit{\boldsymbol{u}} \in \mathit{\boldsymbol{U}}, } \right.\\ \;\;\;\;\;\;\;\left. {\exists {t_{\rm{f}}} \in \left[ {t, T} \right]\left| {\mathit{\boldsymbol{\xi }}\left( {\mathit{\boldsymbol{x}}, {t_0}, \mathit{\boldsymbol{u}}\left( \cdot \right), {t_{\rm{f}}}} \right) \in K} \right.} \right\} \end{array} $ | （9） |
式中：K∈Rn为目标集；tf为达到目标集所需的时间。
本文研究飞机在结冰条件下的飞行安全包线，从理论上确定飞机能够安全飞行的最大可操纵范围，而通过可达性分析得到的反向可达集刚好能够求得飞机可以恢复到安全飞行状态的最大范围，因此，选取反向可达集作为复杂结冰条件下的飞行安全包线。
2.2 可达集计算可达集的求解基本上可以分为两类，包括拉格朗日法和欧拉法[19]。拉格朗日方法利用紧集表示向量场的流向，因此，计算复杂度较高，一般用于计算较高维的可达集；欧拉法主要是以水平集方法为代表，水平集方法用HJB-PDE(Hamilton Jacobi Bellman-Partial Differential Equation)零水平集方程的黏性解表示可达集[20]，该方法将可达性分析与最优控制理论紧密结合起来，能处理复杂的非线性和最优控制策略问题[21]，在航空界应用较为广泛。
2.2.1 水平集方法利用水平集方程求解时间依赖的Hamilton-Jacobi方程零水平集的黏性解，并将Hamilton-Jacobi方程的黏性解作为可达集的隐式表达式。
水平集方程可以表示为
| $ \frac{{\partial \phi \left( {\mathit{\boldsymbol{x}}, t} \right)}}{{\partial t}} + f \cdot \nabla \phi = 0 $ | （10） |
式中：ϕ(x, t)为水平集函数。
目标集K可由水平集函数表示为
| $ K = \left\{ {\mathit{\boldsymbol{x}} \in {\mathit{\boldsymbol{R}}^n}|\phi \left( {\mathit{\boldsymbol{x}}, 0} \right) \le 0} \right\} $ | （11） |
通过计算HJB-PDE方程的黏性解，得到目标集K对应的可达集为
| $ \left\{ \begin{array}{l} \frac{{\partial \phi \left( {\mathit{\boldsymbol{x}}, t} \right)}}{{\partial t}} + \min \left\{ {0, H\left( {\mathit{\boldsymbol{x}}, \mathit{\boldsymbol{p}}} \right)} \right\} = 0\;\;\;\;\mathit{\boldsymbol{x}} \in {\mathit{\boldsymbol{R}}^n}, t < 0\\ \phi \left( {\mathit{\boldsymbol{x}}, t} \right) = \phi \left( \mathit{\boldsymbol{x}} \right)\quad \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\mathit{\boldsymbol{x}} \in {\mathit{\boldsymbol{R}}^n}, t = 0 \end{array} \right. $ | （12） |
式中：
| $ H\left( {\mathit{\boldsymbol{x}}, \mathit{\boldsymbol{p}}} \right) = \mathop {\max }\limits_{\mathit{\boldsymbol{u}} \in \mathit{\boldsymbol{U}}} {\mathit{\boldsymbol{p}}^{\rm{T}}}f\left( {\mathit{\boldsymbol{x}}, t, \mathit{\boldsymbol{u}}} \right), \mathit{\boldsymbol{p}} = \frac{{\partial \phi \left( {\mathit{\boldsymbol{x}}, t} \right)}}{{\partial \mathit{\boldsymbol{x}}}} $ | （13） |
最优控制输入为
| $ {\mathit{\boldsymbol{u}}^ * }\left( {\mathit{\boldsymbol{x}}, \mathit{\boldsymbol{p}}} \right) = \arg \max {\mathit{\boldsymbol{p}}^{\rm{T}}}f\left( {\mathit{\boldsymbol{x}}, t, \mathit{\boldsymbol{u}}} \right) $ | （14） |
表示在状态x下能使哈密尔顿方程H(x, p)达到最大的控制输入。
2.2.2 最优控制输入计算最优控制输入u*(x, p)在飞控系统中的作用是使飞行状态x的轨迹始终保持在反向可达集内部，最终在t=τ时刻达到目标集。
由式(13)可得本文动力学模型的Hamilton函数为
| $ \begin{array}{l} H\left( {\mathit{\boldsymbol{x}}, \mathit{\boldsymbol{p}}, \mathit{\boldsymbol{u}}} \right) = {p_1}\left( { - g\sin \gamma + \cos \alpha \cos \beta \frac{T}{m}} \right) - \\ \;\;\;\;\;\;{p_1}\frac{{\rho S}}{{2m}}{V^2}\left( {{C_{{D_0}}} + {C_{{D_\alpha }\alpha }} + {C_{{D_\alpha }2}}{\alpha ^2}} \right) + \\ \;\;\;\;\;\;{p_2}\left[ { - \frac{g}{V}\cos \gamma + \left( {\cos \varphi \sin \alpha \cos \beta - } \right.} \right.\\ \;\;\;\;\;\;\left. {\left. {\sin \varphi \sin \beta } \right)\frac{T}{{mV}}} \right] + {p_2}\frac{{\rho S}}{{2m}}V\left( {{C_{{L_0}}} + {C_{{L_\alpha }}}\alpha } \right)\cos \varphi - \\ \;\;\;\;\;\;{p_2}\frac{{\rho S}}{{2m}}V{C_Y}\beta \sin \varphi \end{array} $ | （15） |
式中：
1) 若p1 < 0，当
2) 若p1=0，当
3) 若p1>0，当
4) 对于侧滑角β，如果p2sinφ < 0, 则
当考虑复杂结冰条件下飞机气动参数不能完全确定时，可将不确定量Δ视为外界扰动加入最优控制方程，因此，反向可达集RBK, t0(t)可表示为
| $ \begin{array}{l} R_{K, {t_0}}^{\rm{B}}\left( t \right) = \left\{ {\mathit{\boldsymbol{x}} \in {\mathit{\boldsymbol{R}}^n}:\forall \Delta \in D, \exists \mathit{\boldsymbol{u}} \in \mathit{\boldsymbol{U}}, } \right.\\ \;\;\;\;\;\;\;\left. {\exists {t_{\rm{f}}} \in \left[ {t, T} \right]\left| {\xi \left( {x, {t_0}, u\left( \cdot \right), {t_{\rm{f}}}, \Delta } \right) \in K} \right.} \right\} \end{array} $ | （16） |
HJB-PDE方程可以表示为
| $ \left\{ \begin{array}{l} \frac{{\partial \phi \left( {\mathit{\boldsymbol{x}}, t} \right)}}{{\partial t}} + \min \left\{ {0, H\left( {\mathit{\boldsymbol{x}}, \mathit{\boldsymbol{p}}} \right)} \right\} = 0\\ H\left( {\mathit{\boldsymbol{x}}, \mathit{\boldsymbol{p}}} \right) = \mathop {\sup }\limits_{\begin{array}{*{20}{c}} {\mathit{\boldsymbol{u}}\left( \cdot \right) \in \mathit{\boldsymbol{U}}}\\ {\Delta \left( \cdot \right) \in \left[ { - 1, 1} \right]} \end{array}} \frac{{\partial \phi }}{{\partial \mathit{\boldsymbol{x}}}}\left( {\mathit{\boldsymbol{x}}, t} \right)f\left( {\mathit{\boldsymbol{x}}, \mathit{\boldsymbol{u}}, \Delta } \right) \end{array} \right. $ | （17） |
以RCAM简化模型为研究对象，对结冰飞机低空平飞阶段的飞行安全包线变化及运动特性进行分析。正常飞行时飞机及飞行相关参数为：飞机质量m=120 000 kg, 重力加速度g=9.81 m/s2, 飞行高度H=2 000 m, 大气密度ρ=1.006 5 kg/m3，最大推力Tmax=410 920 N，最小推力Tmin=20 546 N，机翼参考面积S=260 m2, 升力系数CL0=1.065 6、CLα=6.072 5，阻力系数CD0=0.159 9、CDα=0.503 5、CDα2=2.117 5，侧力系数CY=-1.6。
3.1 结冰对飞行安全包线的影响首先对确定性结冰进行分析，假设飞机机翼发生结冰，参照文献[22]中结冰程度的区分标准，分别定义3种结冰程度η，即η=0，0.1，0.3，分别对应未结冰、轻度结冰、严重结冰。
计算飞机正常情况飞行时在给定目标集下的反向可达集，并以反向可达集作为此条件下飞机的飞行安全包线，如图 4所示，在可达集范围内，所有飞行状态都可在最优控制操纵下达到目标集。一旦状态超出可达集边界，飞机将无法恢复到正常飞行状态。
|
| 图 4 正常条件下飞行安全包线 Fig. 4 Flight safe envelope in clean condition |
计算飞机未结冰(η=0)发生轻度结冰(η=0.1)、严重结冰(η=0.3)时的飞行安全包线，并进行对比，如图 5所示。图中绿色、粉色和红色包线分别代表未结冰、轻度结冰及严重结冰时飞机飞行安全包线，由仿真图可以看出，随着结冰程度的增加，飞机的可操纵性下降，飞行安全包线不断收缩，在速度轴V方向上收缩最为严重。当飞机发生严重结冰时，飞机的气动参数发生改变，最大失速迎角减小，最小失速速度增大，飞行安全包线向右收缩，飞机发生失速的危险增大。
|
| 图 5 不同结冰程度下的飞行安全包线 Fig. 5 Flight safe envelope under different icingdegrees |
考虑在结冰情形下飞机滚转对于飞行安全包线的影响，分别选滚转角φ=0°, 30°, 60°的情形进行研究，如图 6所示。
|
| 图 6 不同滚转角下的飞行安全包线 Fig. 6 Flight safe envelope at different roll angles |
图中绿色、粉色和红色包线分别代表滚转角为φ=0°，30°，60°时的飞行安全包线，通过对不同滚转角下飞行安全包线的比较可以看出，随着滚转角φ的不断增大，飞行安全包线持续减小，且主要表现在γ轴方向上，飞机的爬升性能降低，这主要因为飞机爬升的力由Lcosφ提供，当滚转角增大时，爬升力减小，爬升性能降低。因此当飞机发生结冰时，应当减少滚转机动，以增大飞机爬升与俯冲的安全操纵范围。
3.2 不确定结冰条件下的飞行安全包线在3.1节中讨论了在确定性结冰条件下飞机飞行安全包线的变化，但在实际飞行中无法完全获取飞机结冰信息，因此，也无法获得飞机准确的飞行安全包线。在此对结冰程度已知，但结冰位置、冰型及分布等存在不确定性情况下的飞机飞行安全包线进行鲁棒性分析。
根据第1节所建模型，选取高斯分布标准差为3σ(3σ置信度为99%)，不确定量Δ分别为0、20%和30%，对轻度结冰(η=0.1)、严重结冰(η=0.3)时的飞行安全包线进行分析，如图 7、图 8所示。
|
| 图 7 轻度结冰条件下飞行安全包线的变化 Fig. 7 Changes of flight safe envelope undermild icing condition |
|
| 图 8 严重结冰条件下飞行安全包线的变化 Fig. 8 Changes of flight safe envelope under severeicing condition |
图 7中蓝色曲面和黑色曲面分别表示飞机在轻度结冰(η=0.1)条件下，不确定量Δ分别为20%和30%时的飞行安全包线。由于结冰位置、冰型及分布等存在不确定性，使得气动导数变化不确定，当考虑Δ=20%时，鲁棒飞行安全包线相较确定性安全包线在φ=0°收缩了30%，当考虑Δ=30%时，鲁棒飞行安全包线相较确定性安全包线在φ=0°收缩了45%。
图 8中蓝色曲面和黑色曲面分别表示飞机在重度结冰(η=0.3)条件下，不确定量Δ分别为20%和30%时的飞行安全包线。当考虑Δ=20%时，鲁棒飞行安全包线相较确定性安全包线在φ=0°收缩了15%，当考虑Δ=30%时，鲁棒飞行安全包线相较确定性安全包线在φ=0°收缩了35%。
通过对比可知，结冰位置、冰型及分布等不确定性因素对轻度结冰时的飞行安全包线的影响大于重度结冰时的影响，这主要是由于确定性轻度结冰时飞机的气动导数变化相对较小，飞机所处于的状态仍然为平衡状态，但是当考虑不确定性时，气动导数变化增大，原平衡点的性质发生了改变，由稳定的平衡点变为不稳定平衡点。因此，当考虑不确定性时，轻度结冰条件下飞行安全包线收缩更加严重。
4 结论建立了考虑横向运动(β、φ)的飞机纵向运动模型，基于可达集理论利用可达性分析研究了飞机在确定性结冰条件下和考虑结冰位置、冰型及分布等不确定性结冰条件下飞机的飞行安全包线，分析了不同结冰程度(η=0, 0.1, 0.3)对于飞行安全的影响，得出以下结论：
1) 以反向可达集作为飞行安全包线相较于传统包线具有明显优势，并且在外部环境发生变化时，仍然能够指导飞行员操纵。
2) 随着结冰程度的增加，飞行安全包线的范围不断收缩，主要体现在速度轴方向上，飞机的气动参数发生改变，最大失速迎角减小，最小失速速度增大, 包线向右收缩。
3) 当考虑滚转角φ对结冰飞机飞行安全包线影响时，随着φ的增大，飞行安全包线主要沿着γ轴方向收缩，飞行安全包线减小。因此，当飞机发生结冰时，应尽量减小滚转角φ同时避免俯冲机动，以提高飞行安全。
4) 当考虑结冰位置、冰型及分布等不确定性结冰情况下的鲁棒飞行安全包线时，不确定量Δ对于飞行安全包线影响较大，且不确定性因素对轻度结冰时的飞行安全包线的影响大于重度结冰时的影响。因此，即使飞机只是发生轻度结冰但是由于结冰位置、冰型及分布等的不确定性，飞行风险仍然较大，飞行员应该保持警惕，适当缩小操纵范围。
| [1] |
屈亮, 李颖晖, 袁国强, 等. 基于相平面法的结冰飞机纵向非线性稳定域分析[J]. 航空学报, 2016, 37(3): 865-872. QU L, LI Y H, YUAN G Q, et al. Longitudinal nonlinear stabilizing region for icing aircraft based on phase-plane method[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(3): 865-872. (in Chinese) |
| Cited By in Cnki (4) | Click to display the text | |
| [2] | REEHORST A, COLANTONIO R O. Examination of icing induced loss of control and its mitigations[C]//AIAA Guidance, Navigation, and Control Conference. Reston, VA: AIAA, 2010. |
| [3] | BRAGGE M B, BROEREN A P. Iced-airfoil aerodynamics[J]. Progress in Aerospace Sciences, 2005, 41(5): 323-362. |
| Click to display the text | |
| [4] | Safety Advisor Aircraft icing[EB/OL]. (2013-05-01)[2019-06-17]. http://www.aopa.org_media/Files/AOPA/Home/Pilot%20Resources/ASI/Safety%20Advisors/sall.pdf. |
| [5] | LYNCH F T, KHODADOUST A. Effects of ice accretions on aircraft aerodynamics[J]. Progress in Aerospace Sciences, 2001, 37(3): 669-767. |
| Click to display the text | |
| [6] | BRAGG M B, BASAR T. Smart icing systems for aircraft icing safety[C]//40th AIAA Aerospace Sciences Meeting Exhibit. Reston, VA: AIAA, 2002: 14-17. |
| [7] |
裘燮纲, 韩凤华. 飞机防冰系统[M]. 北京: 航空专业教材编审组, 1985. QIU X G, HAN F H. Aircraft anti-icing system[M]. Beijing: Aviation Professional Editors Publishing Group, 1985. (in Chinese) |
| [8] |
王健名, 徐浩军, 薛源, 等. 基于极值理论的平尾结冰飞行风险评估[J]. 航空学报, 2016, 37(10): 3011-3022. WANG J M, XU H J, XUE Y, et al. Flight risk evaluation of tailplane icing based on extreme value theory[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(10): 3011-3022. (in Chinese) |
| Cited By in Cnki (2) | |
| [9] | MERRET J M, HOSSAIN K N, BRAGG M B. Envelope protection and atmospheric disturbances in icing encounters: AIAA-2002-0814[R]. Reston, VA: AIAA, 2002. |
| [10] | ZHANG Y, DE VISSER C C, CHU Q P. Online safe flight envelope prediction for damaged aircraft: a database-driven approach[C]//AIAA Modeling and Simulation Technologies Conference. Reston, VA: AIAA, 2016. |
| [11] | ALEXANDRE M B, IAN M, MEEKO M K. Aircraft autolander safety analysis through optimal control-based reach set computation[J]. Journal of Guidance, Control and Dynamics, 2007, 30(1): 68-77. |
| Click to display the text | |
| [11] | MERRET J M, HOSSAIN K N, BRAGG M B. Envelope protection and atmospheric disturbances in icing encounters: AIAA-2002-0814[R]. Reston, VA: AIAA, 2002. |
| [12] |
王爽, 詹浩. 飞行最大可控边界集及其机动边界保护控制[J]. 西北工业大学学报, 2014, 32(4): 523-528. WANG S, ZHAN H. The safe-set aircraft and maneuver-ability envelope protection[J]. Journal of Northwestern Polytechnical University, 2014, 32(4): 523-528. (in Chinese) |
| Cited By in Cnki (0) | Click to display the text | |
| [13] |
刘瑛, 杜光勋, 全权, 等. 基于Hamilton-Jacobi方程的飞行器机动动作可达集分主析[J]. 自动化学报, 2016, 42(3): 347-357. LIU Y, DU G X, QUAN Q, et al. Reachability calculation for aircraft maneuver using Hamilton-Jacobi function[J]. Acta Automatica Sinica, 2016, 42(3): 347-357. (in Chinese) |
| Cited By in Cnki | Click to display the text | |
| [14] |
武朋玮, 李颖晖, 郑无计, 等. 基于可达集方法的结冰飞机着陆阶段安全风险评估[J]. 航空学报, 2018, 39(12): 122139. WU P W, LI Y H, ZHENG W J, et al. Flight risk evaluation based on reachable set method at the phase of icing aircraft landing[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(12): 122139. (in Chinese) |
| Cited By in Cnki | |
| [15] |
周驰, 李颖晖, 郑无计, 等. 结冰飞机着陆阶段飞行安全包线确定及操纵应对策略研究[J]. 航空学报, 2018, 39(12): 122165. ZHOU C, LI Y H, ZHENG W J, et al. Study on flight safety envelope determination and maneuvering coping strategy for icing aircrafts during the landing phase[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(12): 122165. (in Chinese) |
| Cited By in Cnki | |
| [16] |
唐超颖. 应用滑动模态和时间尺度分离的航天器姿态控制[J]. 武汉大学学报(工学版), 2007, 40(6): 117-121. TANG C Y. Sliding mode control of spacecraft attitude based on time scale[J]. Engineering Journal of Wuhan University, 2007, 40(6): 117-121. (in Chinese) |
| Cited By in Cnki) | Click to display the text | |
| [17] | KITSIOS I, LYGEROS J. Launch pad abort flight envelope computation for a personnel launch vehicle using reachability[C]//AIAA Guidance, Control and Navigation Conference. Reston, VA: AIAA, 2005. |
| [18] | LOOYE G, BENNANI S. Description and analysis of the Research Civil Aircraft Model(RCAM): Technical report TP-008-27[R]. Bedfordshire: GARTEUR, 1997. |
| [19] | MITCHELL I. Application of level set methods to control and reachability problems in continuous and hybrid systems[D]. Stanford, CA: Stanford University, 2002. |
| [20] | MITCHELL I, BAYEN A M, TOMLIN C J. Computing reachable sets for continuous dynamic games using level set methods[J]. IEEE Transactions on Automatic Control, 2005, 50(7): 986-1001. |
| [21] | LYGEROS J. On reachability and minimum cost optimal control[J]. Automatica, 2004, 40(6): 917-927. |
| Click to display the text | |
| [22] |
郑无计, 李颖晖, 屈亮, 等. 基于正规形法的结冰飞机着陆阶段非线性稳定域[J]. 航空学报, 2017, 38(2): 567-588. ZHENG W J, LI Y H, QU L, et al. Nonlinear stability region of icing aircraft during landing phase based on normal form method[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(2): 567-588. (in Chinese) |
| Cited By in Cnki |