结冰对民用飞机飞行性能的影响分析

在获取民用飞机适航合格审定证明前,必须通过飞行试验,并在试验中演示验证民用飞机能在CCAR25部附录C确定的连续和间断的最大结冰状态下安全运行。基于此背景,计算并得到某型机两种构型在45min待机冰型下的飞行性能,与无冰状态下对应构型的飞行性能进行对比分析。结果表明,结冰对飞机的飞行性能具有较大的影响。     

一种超声速状态下机载推进系统模型的建模方法

为了满足机载航空推进系统模型对精度、实时性及数据存储量的苛刻要求,提出了一种机载模型建模方法.对传统的机载稳态模型建立过程中使用的相似准则进行了充分讨论,得出了一种具备更高建模精度的相似准则.在新的相似准则基础上,建立了包含泰勒展开式中非线性2阶项的高精度稳态变量模型,有效解决了分段线性稳态变量模型在大包线、变状态条件下精度难以保证的问题;为了解决超声速状态下安装推力与发动机净推力差别较大的问题,建立了考虑外流特性的简化进气道模型,给出了一种计算安装推力的方法,即由简化发动机模型计算安装推力中符合相似准则的部分,由进气道模型计算溢流阻力、放气阻力及发动机推力中不符合相似准则的部分.仿真结果表明:基于新的相似准则所建立的机载模型具有大包线、变状态适应性,输出参数误差均在0.5%以内.      

考虑陀螺效应的转子动力学相似准则

针对考虑陀螺力矩的转子系统缩尺相似模型与原型是否满足动力学相似的问题,采用量纲分析法推导了考虑陀螺力矩的转子-简单支撑系统相似准则,得到了模型与原型各物理量相似比.通过求解不同相似比下Jeffcott转子系统临界转速的解析解,和利用有限元软件ANSYS对某转子系统临界转速和振型相似性进行了数值模拟,验证了所推导的考虑陀螺力矩转子系统相似准则的正确性.研究表明:考虑陀螺力矩、简单支撑和Rayleigh比例阻尼的转子系统,在相似转速下缩尺模型与原型的固有频率具有相似性,并且临界转速相似.      

民机飞行结冰参数的模拟与校准

结冰问题的研究对于民用飞机设计十分重要,采用结冰风洞试验是主要的方法之一,欧美国家有很多成熟的结冰风洞用于型号研究。结冰风洞中需要对很多参数进行控制,其中水雾参数精确控制的难度很大,而且一般通过间接方法来进行测试校准,往往误差较大,导致不同结冰风洞试验结果往往有所差异。文章介绍了在意大利宇航中心的IWT 结冰风洞完成的冰风洞校准试验,对结冰风洞水雾参数校准方法进行了阐述,为民机结冰风洞研究提供参考。     

旋翼桨叶结冰对直升机飞行性能的影响

采用由翼型冰风洞试验推导而来的工程公式,将旋翼结冰模型嵌入离散格式的旋翼气动模型,分别计算结冰后每个桨盘网格区域的翼型升、阻力系数的增量和挥舞系数的增量,建立了结冰后的UH-60A直升机飞行力学模型.结合发动机功率数据,通过改变结冰条件,研究了不同结冰参数对最大平飞速度、最大垂直爬升率和最大爬升率的影响.结果表明:在环境温度为0~-10℃时,降低环境温度会严重降低飞行性能;液态水含量对所选性能指标有一定影响;水滴直径对最大平飞速度的影响较小,对爬升性能的影响可忽略不计;飞行性能指标随结冰时间呈线性下降趋势.      

基于欧拉法模拟旋转帽罩水滴撞击特性

基于欧拉两相流理论,对旋转帽罩水滴撞击特性进行了数值模拟。通过引入耗散函数求解水滴相无黏欧拉方程,实现了欧拉法求解三维旋转帽罩的水滴撞击特性。分析了旋转帽罩转速对水滴撞击特性的影响。发现旋转帽罩转速越大,其表面水滴撞击极限稍有减小。在典型飞行条件与气象条件下,离心力对水滴运动轨迹的影响远小于惯性力对其的影响,旋转帽罩转速对水滴撞击特性影响较小。     

结冰飞机非线性稳定域确定及安全操纵方法

结冰会恶化飞机的动力学特性,造成飞行包线收缩,威胁飞行安全,研究结冰后飞机的非线性稳定域变化对于驾驶员操纵应对策略设计以及飞行安全的提高具有重要意义。以NASA的GTM为案例飞机,首先对飞机气动参数进行多项式拟合,同时结合结冰因子模型,建立了飞机在结冰条件下的纵向通道动力学模型;然后通过分岔分析方法对飞机在不同程度结冰条件和操纵指令下的飞行状态变化进行了研究,并将其用于指导驾驶员操纵,同时考虑到分岔分析方法的局限性,利用微分流形理论确定了飞行系统的非线性稳定域,并将其作为飞行安全边界;最后针对结冰情形,提出将分岔分析方法与微分流形理论相结合共同用于操纵指导,并进行了操纵时域验证。研究结果表明,结冰会使安全边界收缩,在小扰动的作用下都可能使飞行状态超出安全边界。随着结冰程度增加,飞机的稳定性质甚至会发生变化,此时飞行状态将很难维持在原有的安全边界以内,提出了通过指导驾驶员操纵指令变化使飞行状态到达新的安全边界。研究结果对于飞行安全操纵及边界保护都具有一定的指导意义。      

一种起落架载荷谱相似性判别方法

起落架结构的寿命监控对保障其安全性与经济性具有重要作用,但由于地面谱复杂的高低载荷非线性交互作用,难以准确计算单机谱的损伤。以当量损伤计算方法为基础,通过对寿命监控中的单机谱进行初步筛选,判别单机谱与基准谱的相似程度,进而分析损伤计算方法的适用性。提出了一种基于时间序列分析中的动态时间弯曲方法的起落架载荷谱相似性判别方法,进行了基准谱与4个单机谱下的疲劳实验,通过分析损伤计算误差与起落架载荷谱相似距离的关系,验证所提载荷谱相似性判别方法的合理性。      

基于FMECA信息的测试性验证试验样本分配方法

针对目前国内外测试性验证试验中普遍采用的基于故障率的分层抽样方法考虑因素单一,可能导致样本分配不尽合理,而现有考虑多因素的样本分配方法代表性不足,难以实现工程上的应用的问题,对影响样本分配的因素与故障模式影响及危害性分析(FMECA)信息的关系进行了详细分析,提出了基于FMECA信息的样本分配方法。首先定义了单元影响系数和影响因子向量,提出了向量各元素赋值规则,解决了影响因素选取的问题;其次通过基于逼近理想解排序法(TOPSIS)优化的层次分析法确定了权值矩阵,实现了样本分配。最后通过实例验证表明该方法考虑因素全面,运用灵活,更具代表性和工程应用价值。     

Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。