玻璃纤维增强铝合金层板高速冲击损伤容限

为研究玻璃纤维增强铝合金（GLARE）层板高速冲击损伤容限,对单次、多次冲击载荷下GLARE层板的损伤特性进行了高速冲击试验和数值仿真研究。采用一级气体炮,在GLARE层板靶板的中心位置、边位置、角位置进行弹道冲击试验,获取弹道极限和损伤模式,然后结合数值仿真剖析动态响应特性。结果表明,弹道冲击条件下,GLARE层板主要通过塑性变形、金属层裂纹、脱胶、复合材料纤维脆断等损伤模式吸收冲击能量。边界约束效应对GLARE层板冲击损伤特性具有显著影响,主要表现在：约束效应导致不同位置冲击下GLARE层板损伤模式不同,包括成坑、金属裂纹和冲塞等特征;角位置冲击条件下,GLARE层板的弹道极限速度明显低于中心位置冲击结果;重复冲击条件下,角位置比与边位置和中心位置更容易发生穿透。     

飞行器大迎角下俯仰静、动导数的数值计算

 采用一种适用于非定常流动计算的逐次超松弛（SOR）时间推进方法耦合求解非定常Navier-Stokes方程和强迫运动方程,数值模拟了NACA0012翼型、弹道外形(HBS)和有翼导弹等飞行器在不同初始迎角下做强迫俯仰振荡的黏性动态流场,给出动态气动力及力矩系数的时间历程;采用积分法获得飞行器的静导数和动导数,将计算结果与实验结果及半经验理论预测方法进行比较并分析飞行器在大迎角下动态特性的非线性特征;计算结果与实验值相吻合,表明本文采用的非定常数值模拟方法在大迎角下能够较准确地模拟飞行器的动态特性,具有工程应用价值。     

非线性气动恢复力与飞行器动稳定性

复杂流动中,飞行器非定常气动数值计算后的参数辨识、飞行器模型的风洞定常测量、风洞动态(风洞自由飞)测量均显示了飞行器非线性气动刚度项(气动恢复力项、如俯仰力矩 攻角曲线Cm θ)将出现局部静不稳定( Cm/ θ>0)。　　飞行器气动恢复力对运动稳定性影响以往讨论较少,本文讨论了非线性气动恢复力局部出现静不稳定条件下的动稳定特性,其特点是突跃性、单边性与不确定性。飞行器飞行遥测数据分析支持了上述动稳定特性的描述。     

火星探测器气动热环境和气动力特性的数值模拟研究

针对探测器进入火星大气过程中的高温气体非平衡效应及其气动热环境和气动力特性的精确预测问题,考虑火星大气高温混合气体的化学反应、束缚电子能激发效应、分子振动能激发效应和非平衡效应,以及不同表面催化特性和表面温度条件的影响,建立了数值模拟热化学非平衡流场的计算方法,并通过与平头圆柱试验模型和MESUR火星探测器的试验与计算数据对比分析,验证了计算方法的可行性与有效性,然后对典型火星探测器的气动热环境和气动力特性进行了数值模拟分析。研究结果表明:1)对于本文的计算状态情况,火星探测器流场中化学非平衡效应非常严重,存在一定的热力学非平衡效应;2)热力学模型、化学反应模型和表面催化特性对气动热计算有着重要影响;3)气体模型对气动力特性的预测影响较大,表面温度条件和表面催化特性对气动力系数的计算影响较小。     

高超声速飞行体亚密湍流尾迹RCS特性的相似规律研究

再入或临近空间高超声速飞行器在大气层中飞行时,会形成等离子体绕流和尾迹,将严重影响飞行器的电磁散射特性,从而使得临近空间高超声速飞行器表现出与大气层内常规飞行器和大气层外轨道飞行器在电磁散射特性上的明显不同,对飞行器的探测、跟踪带来新的挑战。弹道靶是高超声速飞行器尾迹电磁散射特性研究的重要设备,为了由弹道靶试验结果外推出真实飞行条件下亚密湍流尾迹的RCS,本文从高超声速流场模拟的双缩尺率和亚密湍流尾迹RCS模拟的Born近似出发,推导出真实飞行条件下和地面弹道靶试验之间亚密湍流尾迹RCS模拟的一种相似规律,推导结果显示当D1/D2=η,V∞1=V∞2,ρ∞1/ρ∞2=1/η,ω1/ω2=λ2/λ1=1/η时,雷达散射截面服从σ1/σ2=η4的相似规律。推导出的相似规律可以为超高速飞行器亚密湍流尾迹RCS研究提供技术支持。     

高超声速全机外形气动加热与结构传热快速计算方法

发展了一种无黏流场解与工程计算方法相结合的高超声速全机外形气动加热与结构传热快速计算方法。该计算方法结合了三维块结构网格无黏流场数值计算技术可处理复杂外形流动的优点与工程计算方法效率高的特点,将气动热的计算简化为绕飞行器的无黏外流(边界层以外)数值解和边界层内热流求解两个部分,同时耦合了防热结构传热计算模型、高温化学非平衡热效应估算方法以及弹道状态动态插值方法,可用于快速计算与分析三维复杂外形高超声速飞行器在弹道飞行状态下全机热环境参数、防热结构内温度场等随飞行时间的变化特性。以RAM-CⅡ、类Ⅹ-37B等典型高超声速飞行器为研究对象,在设定的飞行条件及热防护方案下,进行了气动加热与结构传热问题的求解,给出了全机表面热流密度与防热结构材料温度的时变特性。结果对比表明,所发展的方法具有快速、高效的特点,且计算精度可满足工程设计初期选型需求,可为高超声速飞行器的热防护系统初期设计及热环境特性快速计算分析提供技术支持。     

高超声速飞行器高温流场数值模拟面临的问题

随着高超声速飞行器目标光辐射和电磁散射特性研究的发展和深入,高温流场特性日益引起人们的关注。由于高温流场特性研究中涉及到非常多的复杂气动现象,如气动加热、烧蚀、辐射、燃烧、化学反应以及湍流等,因此其数值模拟面临着诸多挑战。这里基于连续流计算流体力学(CFD)技术和稀薄气体蒙特卡罗直接仿真(DSMC)方法,从化学物理模型建模、方法稳定性与数值求解效率出发,分析了高超声速飞行器外部绕流、尾迹和发动机喷焰三方面的流场特性数值模拟在不同弹道、热防护手段和飞行流域环境下所面临的问题。在此基础上提出了数值求解技术和化学物理模型建模今后需要发展的方向,为有效提高高超声速高温流场特性数值模拟效率、增加流场特性预测精度提供了指导,从而为研究流场对高超声速飞行器目标光辐射和电磁散射特性影响提供有效的基础数据。     

一种临近空间飞行器静/动态气动特性研究

临近空间飞行器典型气动布局给其气动特性尤其是飞行稳定性带来诸多问题,危及飞行安全。为此,针对类HTV-2飞行器布局方案,开展了高超声速飞行时的气动特性研究。通过数值计算和试验结果的对比可见,两者具有较好的一致性,同时看出飞行器具有较高的升阻比,纵横向为静态稳定,并且给出了气动特性在不同飞行速度和高度条件下的变化规律;飞行器单自由度运动具有动态稳定性特性,但是通过风洞试验进行的两自由度耦合运动动态特性试验可以看出,在强迫俯仰运动下滚转出现振荡运动。     

飞行器动态稳定性参数计算方法研究进展

动态稳定性参数（简称动导数）是飞行器控制系统设计、飞行器动不稳定发生边界分析及相应动态稳定性判据研究的关键气动参数。在对飞行稳定性问题进行概述的基础上,介绍飞行器动态稳定性参数数值模拟的国内外研究进展。并按照理论方法、工程近似方法及计算流体力学（CFD）模拟方法的动导数发展方向对近年来主要的动导数计算方法进行了综述,评价了各种动导数预测方法的优缺点,指出了动导数数值模拟在理论基础、非定常气动力建模、预测方法精度和效率等方面存在的问题。最后对动导数数值模拟的发展趋势进行了展望。     

火星探测器气动外形/弹道一体化多目标优化

针对火星探测器概念设计阶段的需求,提出了融合气动外形、弹道和开伞条件的一体化多目标优化设计方法。首先建立了火星探测器进入段三自由度弹道运动方程,基于修正牛顿理论推导了适用于具有较大半锥角球锥外形的气动参数估算模型,采用Sutton-Graves公式计算了驻点热流密度。以开伞高度、总吸热量和容积率为目标函数建立了火星探测器气动外形/弹道一体化多目标优化模型,采用基于分解的多目标进化算法(MOEA/D)进行求解计算并与参考设计进行了对比。数值结果表明:多目标优化方法提供多个三目标均优于参考设计的Pareto最优解,为火星探测器的概念设计提供了一定的参考依据。     

半球形头弹斜撞击下铝合金薄板的失效仿真分析

为分析斜撞击下撞击角度对6061-T651铝合金板失效特性的影响,首先用一级轻气炮发射半球形头弹正撞击铝合金靶板,得到试验数据.其次根据试验工况建立弹靶模型并验证模型的有效性,得到数值模拟结果,不仅能很好地预测靶体的弹道极限速度,还能够较好地分析弹体初始速度—剩余速度变化趋势.最后建立5种不同撞击角度下的弹靶模型,分析...     

面向弹道优化的高超声速再入飞行器模态稳定性分析

为研究高超声速再入飞行器沿弹道的自由扰动运动的稳定性,考虑大气密度随高度的变化和引力梯度,建立了高超声速无动力再入纵向动力学小扰动线性化方程,然后获得转移矩阵和特征方程,在此基础上进行沿弹道的纵向模态分析。利用二次曲线及基于类型函数和形状函数（CST）的方法提出升力式高超声速飞行器气动布局,并采用工程估算方法获得飞行器气动特性数据。针对最大射程、最小射程和跳跃弹道等典型再入弹道进行沿弹道的模态稳定性分析,得到高超声速再入弹道高度模态、沉浮模态和短周期模态稳定性沿弹道的变化特征。从稳定性的角度,对弹道优化提出建议：应避免所设计的弹道产生太大的跳跃,即使是牺牲一些射程上的性能,因为跳跃会使短周期模态和沉浮模态产生更多的不稳定特征根。     

火星进入器高超声速化学非平衡特性数值计算研究

火星探测活动是当前国际深空探测领域的热点,而探测器的进入、下降和着陆阶段的气动特性预测是关系其布局和热防护系统设计的关键问题之一。针对火星进入器高超声速进入过程,发展了火星环境的高超声速化学非平衡流动数值模拟方法,建立了进入段气动特性的有效预测技术。以探路者号火星进入器防热大底为研究对象,采用5组分6化学反应模型开展典型弹道点上的化学反应流动数值模拟。结果表明:CO_2在激波层内发生较大程度的离解反应,火星进入器高超声速化学非平衡特性显著。因化学反应作用强烈,激波层内热传导所致热流量相对较低,而组分扩散造成的热流量占一定比重,计算所得的驻点总加热量与经验公式吻合良好。火星环境的化学非平衡特性的有效预测可为我国未来火星探测任务提供技术支持。     

双弹头牵连体再入动稳定性分析

双弹头牵连体再入过程的动稳定性分析是有关工程设计中的重要环节。依据六自由度弹道方程 ,建立了双弹头牵连体再入的动力学模型 ,利用弹道与气动力的耦合 ,初步分析了双弹头牵连体再入的动稳定特性。计算结果表明 ,绳索牵连减弱了弹头再入的动稳定性 ,由动稳定弹头组成的双弹头牵连体在某些条件下可能出现打转、横飞等动不稳定现象。     

稳定性导数计算方法研究

稳定性导数的准确计算对于飞行器的操稳特性具有重要意义.应用计算流体力学(CFD)方法中的非结构化动网格技术建立能够模拟飞行器做周期性俯仰运动的强迫振荡法,以国际动导数标模Finner导弹为验证算例,获得不同马赫数下俯仰力矩系数的迟滞环曲线,进而计算Finner导弹在不同马赫数下的静稳定性导数和动稳定性导数.结果表明:本文计算得到的俯仰静稳定性导数与试验数据非常接近;在亚音速和超音速范围内,动稳定性导数计算结果与试验值很接近,但在跨音速范围内,本文计算结果与试验曲线的规律性一致,但误差较大.     

动导数数值预测中的相关问题

动导数是飞行器动态稳定性分析、弹道设计和控制系统设计的重要参数,其预测方法主要有工程近似方法、数值模拟和风洞试验。基于动导数、交叉导数的概念,介绍了强迫振荡法、自由振荡法等预测动导数的数值模拟方法,重点就超声速、高超声速动导数数值预测中存在的问题开展讨论,包括时间步长、子迭代步数的选取;振荡频率对辨识结果的影响;强迫振荡法和自由振荡法辨识动导数可能存在差异的原因分析;交叉导数和交叉耦合导数的辨识问题以及飞行器构型对动导数预测的影响等。并结合算例进行了具体分析,在总结现有研究经验的同时,针对当前研究中存在的困惑和难题,提出了相应的改进建议。     

高超飞行器内流道激波振荡问题的数值研究及试验验证

采用数值模拟和风洞试验的方法,对吸气式高超声速飞行器盲腔状态的流动特性进行了研究。采用＂双时间步＂方法进行了内外流一体化的非定常数值模拟,利用彩色纹影系统对高超声速飞行器前体流场进行显示,并采用动态压力传感器测得了飞行器内流道的壁面压力随时间的变化。结果表明：盲腔构型在高超声速飞行中会出现周期性的激波振荡现象。数值模拟所得流场变化特征、内流道壁面压力振荡周期和壁面压力变化趋势与试验结果吻合良好。     

高速高温流场电子能非平衡的数值模拟

采用描述电子能非平衡的三温度模型,结合11组分空气的化学反应模型,对多种高速高温热化学非平衡流场开展数值模拟,并与描述电子能平衡的两温度模型结果进行对比,研究电子能非平衡对高超声速流场特性的影响。圆球弹道靶试验算例表明电子能非平衡不影响激波脱体距离。RAM-C II飞行器的4个飞行工况算例表明,尽管两温度和三温度模型结果存在差异,但二者电子数密度分布的趋势和量级接近,均可与飞行试验数据保持一致,其中三温度模型的预测效果更好。FIRE II飞行器极高温流场模拟结果显示,电子能非平衡几乎不影响飞行器表面的对流传热。     

飞行器动态气动特性研究的问题

本文介绍了飞行器动态气动特性研究的几个问题,包括:大攻角动稳定性的提法;大攻角动稳定性风洞试验技术;大攻角非定常气动力的典型特性;和大攻角非定常气动力风洞试验技术等.这些问题关系到飞行器安全飞行包线和飞行器的机动性和可控制性,因而,在新一代先进飞行器的研制过程中越来越受到重视.     

一种近空间高超声速飞行器滚转稳定性研究

研究近空间高升阻比构型高超声速飞行器面临的复杂稳定性问题,对其机动飞行安全和操控具有重要意义。为此,通过静、动态风洞试验和数值模拟,开展了类HTV-2飞行器布局典型马赫数和迎角下,滚转方向静稳定性、单自由度和两自由度俯仰/滚转耦合运动动稳定性研究。结果表明：飞行器滚转方向为静稳定;自由振动的滚转单自由度运动为衰减振动曲线,动导数为负值,表现为动稳定;强迫俯仰/自由滚转两自由度耦合运动情况下,由于惯性力矩和非对称滚转力矩共同作用,导致滚转出现振幅周期性变化的极限环运动,甚至失稳发散。