基于流固耦合的激波/边界层干扰作用下壁板颤振特性

为了解激波/边界层干扰作用下壁板气动弹性及其对流动分离的影响，采用自主开发的双向流固耦合求解器，对不同激波冲击位置下壁板的振动响应和流动特性进行了数值模拟研究。壁板几何非线性运动方程采用有限差分法求解，基于有限体积法求解Navier-Stokes方程组，对流通量采用MUSCL和AUSMPW+格式离散，双向流固耦合采用交错迭代算法。研究结果表明：激波/边界层干扰作用下壁板振动位移先增大后减小，经若干振荡周期后达到稳定颤振状态，呈现二阶振动模态，壁板变形相对于激波冲击位置呈现非对称性，壁板前部分的振幅始终小于壁板后部分；激波冲击位置可显著改变壁板的颤振振幅、频率及分离区长度，当激波冲击位置靠近壁板两端时，壁板振动最终收敛达到静稳定状态；壁板振动响应与流场特征不随激波冲击位置的改变而单调变化，对于激波冲击位置x/a=0.35工况，壁板颤振可有效抑制激波/边界层干扰流动分离。     

冷却塔群塔干扰条件下脉动风压分布模式

冷却塔群塔为常见的群体建筑物组合形式,存在复杂的塔群干扰效应.传统的群塔干扰效应研究通常面向塔筒等效静风荷载,忽略了群塔组合对于脉动风压分布模式干扰效应的量化研究.为了研究群塔组合条件下脉动风压的干扰分布模式,选取具备多种气动干扰特点的八塔组合矩形布置方案为研究对象,对群塔进行刚性模型同步测压风洞试验,获取了不同风向角条件下塔筒表面的脉动风压分布.在"迎风来流"、"侧风来流"和"背风来流"等三种典型来流模式下,结合风向角的变化将脉动风压分布形式总结为"V"型、"M"型、半"V"型+半"M"型三类,探讨了干扰影响脉动风压分布的机理,发现可将脉动荷载的干扰效应分为"弱干扰"、"非对称强干扰"和"对称强干扰"三种类别;在量化表达脉动风压干扰分布方面提出了脉动风压分布的九项式拟合公式,并按照三类分布模式给出了脉动风压的推荐表达式.     

单组元300N发动机低温试验研究

为探究低温环境下单组元300N发动机的工作特性，揭示影响发动机低温性能的主要影响因素，以300N发动机为试验对象，开展了模拟飞行工况的发动机低温试验。给出了低温试验研究方法，分别从温度差异对发动机性能影响、催化剂活性差异对发动机低温启动特性影响和低温对电磁阀响应特性影响等方面获得研究结果。结果表明，低温是影响发动机低温性能的主要影响因素，-48℃条件催化剂无法完成推进剂的催化分解，发动机发生爆炸；-30℃条件下起活时间为80.5～87.5ms，发动机可正常启动，且启动温度与起活时间呈指数关系；催化剂批次差异也对发动机低温工作性能产生一定影响，不同批次催化剂低温起活时间的差异可达91ms；低温试验过程中，电磁阀的关闭受到低温推进剂粘性和背压的影响，产生了明显的迟滞现象，延迟时间约100ms，对发动机在轨的精准控制存在一定影响。     

轻量化红外成像制导镜头设计与仿真

根据红外成像制导系统的轻量化、小型化、低成本发展需求，针对制导系统中的光学镜头进行轻量化设计。通过分析超透镜的波前调控原理和加工工艺，利用传输相位调控原理设计了基于硅和氟化钡材料的中红外波长超透镜，工作波长3 μm,镜头厚度小于1 mm，口径50 mm时重量小于10 g，热膨胀系数较低环境适应性好。利用时域有限差分法对超透镜进行了仿真，通过设置对称边界条件缩短了仿真时间。仿真结果表明，当超透镜数值孔径为0.45时，聚焦效率为85%，焦斑半高宽3.4 μm。     

左行运动激波主导的管内流动特性研究

为了分析左行运动激波主导的管内流动特征，本文采用非定常数值仿真方法，对亚声速进口条件下等直管道内左行运动激波传播与演化特性、左行运动激波/边界层干扰特征开展研究。研究结果表明：在出口周期性强压力脉动干扰下管内存在连续的左行运动激波，该左行运动激波传播特征具有相似性，激波强度、传播速度按幂函数规律衰减。气流经过左行运动激波后总压、总温、静压阶跃式升高，随后受膨胀波影响气流总压、总温、静压下降；左行运动激波/边界层干扰诱发形成翼型回流区，该回流区随运动激波强度衰减逐渐减小。理论与数值分析表明存在左行运动激波后速度为零和运动激波两侧总压相等的两个临界状态。波前马赫数低于临界值或左行运动激波强度高于临界值时，左行运动激波后为倒流、波后总压高于波前。     

风雷软件LES开发设计与验证

基于国家数值风洞风雷软件开源框架，设计开发了LES(Large Eddy Simulation)湍流模型，主要包括Fourier谱/有限差分方法 LES求解器和有限体积/有限差分方法 LES求解器。简要介绍了采用的不可压缩流动求解的投影法、Fourier谱/有限差分混合方法、亚格子模型等理论方法，给出设计的软件框架和计算流程，尤其说明采用的松/紧2种耦合模式。通过数值模拟不可压缩槽道湍流、亚临界雷诺数圆柱绕流、NACA0012临界攻角的低频振荡算例，验证求解器的计算精度和复杂湍流模拟能力。基于风雷开源框架设计的开源LES模型，具备高精度数值格式、亚格子模型、湍流统计等通用模块，可为国内学者提供一个LES湍流模拟研究的开放平台。     

轨控大干扰下的系统角动量管理

摘要： 针对利用角动量交换装置系统实现轨控期间姿态维持的情况，研究当轨控推力器存在大力矩扰动时系统角动量管理方法.结合对地运行航天器轨道运动特性，分析轨控干扰力矩的积累角动量变化规律；考虑角动量管理装置系统的角动量存贮容量，基于角动量积累规律提出了一种结合偏置角动量建立与对称分布轨控位置选择的系统角动量管理方法，实现了对称轨控位置积累角动量互相抵消的自平衡效果，极大程度地提高轨控效率且避免系统角动量饱和现象.所提出的方法适用于喷气欠驱动控制系统，方法的有效性通过了数学仿真验证及在轨型号应用.     

通信拓扑切换下的多飞行器协同拦截方法

针对通信拓扑切换条件下的多飞行器协同拦截问题，提出了一种基于扩张状态观测器的协同制导方法。建立协同制导设计模型，将协同拦截问题转换为视线稳定条件下的剩余飞行时间调节问题。为解决机动目标状态不确定的问题，将目标的状态视作扰动，设计扩张状态观测器来估计机动目标的状态，并在制导律中对目标的机动进行补偿。利用有限时间一致性理论进行一致性控制协议的设计，实现各飞行器剩余飞行时间的有限时间一致，并利用Lyapunov稳定性理论分析通信拓扑切换情况下闭环系统的有限时间稳定性，给出了系统一致收敛时间。仿真结果表明，在通信拓扑变换的情况下，设计的观测器能够有效估计目标状态，且协同制导律能够满足对剩余飞行时间的控制要求，进而实现协同拦截。     

基于CZT的变时间窗时变参数在线辨识方法及应用

针对时变参数的在线辨识问题,提出了基于CZT的变时间窗时变参数在线辨识方法,通过滑动时间窗以及基于方差阈值的重启设计,提高了现有频域辨识方法对时变参数的辨识速度,并且有效避免了数值微分对噪声的放大问题。最后,基于飞机系统的执行机构故障,对该方法进行仿真验证,并与传统辨识方法进行对比,说明了该方法对时变参数有较强的辨识能力,验证了该方法的有效性。     

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

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