基于单纯B样条的航空发动机机载稳态模型研究

为了减小航空发动机稳态建模的模型误差、降低复杂度及提升其实时性,提出了一种基于单纯B样条函数的航空发动机稳态模型建模方法。该函数是局部多项式基函数的线性组合,因此求解该函数为线性回归问题,通过运用广义最小二乘方法来求解B系数,从而提高计算效率和提高模型精度。最后建立了基于该算法的二维和四维涡扇发动机稳态模型,并分别与相同建模样本条件下的多输入多输出约简迭代最小二乘支持向量机稳态模型进行了比较,表明了单纯B样条建模方法不仅继承了B样条的算法复杂度低、存储数据量小和实时性好等优点,同时避免了最小二乘支持向量回归机不能拟合大样本数据的缺点,且拟合效果优于最小二乘支持向量机。     

镁粉对铝/冰燃料燃烧特性的影响

为改善铝/冰燃料的点火燃烧性能,通过添加不同比例的镁粉制备成铝/冰基燃料,研究其对铝/冰燃料燃烧特性的影响。实验通过接触法测得铝/冰基燃料燃烧温度并采用高速摄影系统观测其燃烧现象。研究结果表明,铝/冰燃料中添加镁粉后在100~200℃即出现较为明显的氧化放热反应,镁粉可以明显降低铝/冰基燃料的起始反应温度,加快反应速率。在稳定燃烧阶段,随着镁粉添加量的增加,铝/冰基燃料燃烧剧烈程度降低。通过对燃烧产物的分析发现,添加镁粉后的燃烧产物形貌与未添加镁粉的产物形貌存在明显差别,且未完全反应的铝粉含量较多,燃烧相对不充分。     

地外天体起飞过程真空羽流导流力热效应研究

为研究地外天体起飞真空羽流对探测器分离产生的力热扰动,使用计算流体力学-直接模拟蒙特卡洛(CFD-DSMC)耦合计算模型对锥面导流的真空羽流场进行了计算。采用组分输运模型计算三维连续流场,并获取当地的克努森数作为判断连续流和离散流的依据。使用基于分子动力学的直接模拟蒙特卡洛方法(DSMC)计算离散流场,采用可变软球(VSS)碰撞模型和Larsen and Borgnakke传能模型计算离散流分子间的能量传递,将计算结果与试验进行了对比,验证了计算方法的可靠性。研究结果表明,A器受到的侧向干扰力矩为62N·m,底部受到的最大压力为100Pa,最大热流密度为100k W/m2;B器受到的侧向干扰力矩为558 N·m,表面最大压力为8k Pa,最大热流密度为600k W/m2,喷口与导流装置顶面距离为400mm时,激波已移出喷管内部。     

AZ31B镁合金大应变循环变形行为研究

为了研究AZ31B镁合金在大应变幅条件下的变形机制,开展了该合金在7.嬲应变幅条件下的循环行为研究.结果表明:在拉伸阶段的最大应力值随着循环周次的增加而减小,而在压缩过程中的最大应力值随着循环周次的增加而增大,在整个循环过程中材料呈现循环应变硬化特性,拉应力是导致循环应变硬化的主要原因;随着循环周次增加,滞回曲线的不对称性基本不变.真应力-真应变滞回曲线在卸载和反向拉伸阶段出现3个拐点.在压缩过程中发生{10-12}孪生,反向拉伸过程发生去孪生行为,包申格效应对去孪生行为具有较大影响.研究表明:孪生-去孪生是大应变幅循环变形的主要变形机制;对拉伸、反向压缩过程的变形特征及机制的分析,可为低周疲劳行为的研究提供参考.     

基于流-固耦合传热的热气防冰系统干空气飞行蒙皮温度场计算研究

以机翼热气防冰系统为研究对象,建立了包含热气防冰系统防冰腔内外流场对流换热和固体结构导热的三维稳态流-固耦合传热物理模型,对整个计算区域生成混合网格,边界条件为第三类边界条件,采用计算流体力学方法以 FLUENT 软件为工具,对干空气飞行状态下流-固耦合传热模型进行了求解,获得防冰腔蒙皮内外表面对流换热系数分布和温度场结果,并对计算结果进行了分析。结果表明：防冰腔铝合金蒙皮沿展向和厚度方向导热显著,温度分布较均匀,防冰引气温度为200℃时,防冰腔蒙皮内外表面上最高温度为101℃,最低温度为21℃,3 mm厚的蒙皮同一点处内外表面最大温差仅为4℃,防冰腔排气口处气体的平均温度为63℃。热气防冰系统蒙皮温度场计算方法和计算结果,能够为热气防冰系统干空气飞行试验设计和测试中温度传感器的选型与布置提供依据。     

高度欠膨胀声速射流的自相似特性

构建了高压气体射流模型和高分辨率的计算网格,对喷压比（NPR）为5.60、流动雷诺数Re为105量级的高度欠膨胀射流进行了三维大涡模拟（LES）计算.讨论了时均的射流近场结构,发现大涡模拟成功捕捉到了高度欠膨胀射流近场的典型波系结构,并与文献结果吻合较好.研究着重定量考察流向速度和氮气质量分数的分布规律,以揭示高度欠膨胀射流的自相似特性.结果表明:在射流核心区之后,不同流向位置上流向速度和氮气质量分数沿径向的分布呈现出与亚声速射流类似的自相似特性.但射流流场开始呈现自相似的位置仍然为超声速,这是高度欠膨胀射流自相似特性的独有特征.提出了高斯拟合公式,在靠近和远离射流中心线的区域内均可较好地表征射流的自相似特性.此外,还考察了射流剪切层的发展特征,比较了计算得到的射流剪切层增长率与先前实验测量结果的差异,并分析了其中可能的原因.      

基于操纵面故障影响估计的安全飞行轨迹优化

操纵面故障将改变飞机受到的气动力和力矩,降低机动性能和控制能力,基于传统质点模型的航迹规划方法难以解决此类飞机安全飞行轨迹优化问题.因此,引入由故障引起的侧滑和侧力,估计故障后气动力系数,建立故障飞机模型;针对故障飞机控制和机动能力下降的问题,在原故障模型的基础上引入控制量变化率,构建轨迹优化对象模型,并将控制量变化率最小作为优化指标;采用高斯伪谱法优化计算安全飞行轨迹.仿真实验以副翼和方向舵故障为例进行设计,结果满足运动方程及各控制量、状态量约束,并能实现平稳飞行.     

太阳风——磁层相互作用的磁流体力学数值模拟研究

磁层位于地球空间的最外层, 太阳风与磁层的相互作用是空间天气变化因果链中承上启下的关键环节, 是揭示地球空间天气基本规律的关键科学课题. 地球空间由于时变、多成分、多自由度的关联相互作用,使得传统的理论分析变得非常困难. 数值模拟作为近几十年发展起来的一个新的研究手段,对地球空间的理论和应用研究产生了深刻的影响. 国际上磁层的全球MHD数值模拟工作开始于20世纪70年代末, 最初的研究局限于二维空间. 由于磁层内在的三维特性, 20世纪80年代三维MHD数值模拟工作兴起. 本文简要说明了三维全球磁层MHD (磁流体力学)研究的特点及现状, 给出了三维全球磁层模型的基本框架, 综述了行星际激波与磁层相互作用、大尺度电流体系、重联电压和越极电位、磁层顶K-H不稳定性等方面的太阳风--磁层相互作用的MHD数值模拟的研究进展.      

粒子滤波及其在航天器交会对接相对导航中的应用

在全面分析粒子滤波原理的基础上,提出一种改进高斯粒子滤波方法.该方法利用确定性采样滤波算法进行时间更新,替代高斯粒子滤波算法中的随机采样过程;另外,针对厚尾噪声情况,利用鲁棒统计方法对确定性采样滤波方法进行鲁棒性改进,并将其应用于所提出的改进高斯粒子滤波.将粒子滤波算法应用于交会对接相对导航问题,仿真结果表明,在多种测量噪声情况下,改进高斯粒子滤波较其他粒子滤波,能够在不过多损失估计精度的同时有效降低计算量.文中的研究成果为将粒子滤波应用于航天器导航问题提供了理论参考.     

对接口匹配的航天器双主动最优交会研究

研究了两航天器协同轨道机动（双主动）完成近距离交会任务的最优控制问题。在考虑航天器姿态变化、对接口位置及路径约束的情况下,构建了完整的6自由度、26状态的双主动最优交会数学模型。利用高斯伪谱法分别将燃料总消耗最少和交会时间最短两种最优问题离散为大型非线性规划问题,而后应用SNOPT软件包进行了求解。在此基础上通过大量数值计算分析总结了不同初始参数对最小燃料消耗和最短交会时间的影响规律,并与主被动交会形式进行了对比。结果表明,当两航天器质量接近时,双主动交会通常可明显减少燃料消耗,或缩短交会时间;而当质量差距较大时,双主动最优交会逐渐退化为主被动最优交会。