考虑J_2摄动影响的零射程线计算方法研究

考虑到地球扁率J_2摄动的影响,在非正交分解法的基础上,利用自由段弹道地心距迭代解析解,给出了一种弹道飞行器零射程线的计算方法。数值仿真结果表明,该方法具有很高的精度,且计算速度较快,在闭路制导能量耗散和多弹头进攻弹道设计中具有较高的应用价值。     

平／钝头外形俯仰动态特性对落点姿态的影响分析

采用自由振荡法数值模拟了平头、钝头外形的超声速俯仰振荡时间历程,并应用奇异分解线性最小二乘法辨识稳定性导数,得到动导数随马赫数和攻角非线性变化的规律.平头外形在较高马赫数和中小攻角范围内存在动不稳定现象,随着再入时马赫数的降低,动不稳定的攻角范围不断缩小,直至演化为动稳定的;小钝头外形则没有动不稳定现象.蒙特卡洛弹道仿真表明,平头外形初始时刻俯仰角的扰动会经历一个先增长然后衰减的过程,最终落地时俯仰角小于1°,满足设计要求,但其落地速度有所降低,可能影响该外形对靶标的侵彻效果.     

基于样条约束“EMBET”再入轨道测量数据融合方法

针对再入轨道最小二乘交会解算方法数据利用率低、轨道解算精度差的问题,提出并实现了基于样条约束“EMBET”的再入轨道测量数据融合处理方法.该方法认为再入飞行轨道在时序上是相关的,可以利用三次B样条函数精确表示,建立了关于样条函数参数和测量系统误差的测量模型,从而大量压缩了待估参数数量,准确自校准系统误差,提高了轨道估算的精度和稳定性.     

高超声速再入飞行器气动特性的快速预示—局部方法的推广应用

本文将计算高超声速稀薄气流过渡领域中气动特性的局部方法,推广应用到连续介质中弹头型高超声速再入飞行器气动力特性的快速估算。由激波风洞中M_∞=9.9时,一个8°钝锥的气动力测量结果,导出这一实验条件下的领域系数,并以此来估算不同锥角、不同钝度比及不同外形弹头型再入飞行器的气动力和力矩系数,其结果与无粘数值解及实验结果作了比较,在攻角2°～14°范围内吻合得很好。局部方法可用于弹头型高超声速再入飞行器气动特性的快速预示。     

采用几何方式组合壁面压力/马赫数分布的弯曲压缩面

根据给定壁面参数分布规律反设计压缩面的方法,通过几何方式组合研究了壁面采用压力/马赫数复合分布规律的弯曲压缩面,分析了压缩面的参数分布和性能特点.在与参考二维进气道相同约束条件下,将弯曲压缩面应用于二维进气道,分析了其性能并与参考二维进气道进行了对比.数值研究结果表明:几何方式组合能充分发挥壁面按单一参数设计所得弯曲压缩面的优势,使壁面参数分布得到有效改善并趋于更加合理.同参考二维进气道相比,几何方式组合进气道的外压缩面长度有所增加,但在来流马赫数在4~6范围内它的喉道总压恢复系数均得到显著提高.来流马赫数为6时,几何方式组合进气道的喉道总压恢复系数提高10%;来流马赫数为4时其喉道总压恢复系数提高9.6%.     

固体运载器发动机内弹道设计方法

为研究对固体运载器总体性能提升有利的发动机内弹道设计方法,建立了发动机内弹道计算模型,并选取两个燃面设计参数表征内弹道曲线构型。针对某三级固体运载器进行了仿真分析,分别以发动机冲质比最大及运载器射程最大为目标对发动机的内弹道曲线构型进行优化设计,并对这两种设计方法进行对比研究。研究结果表明,发动机性能最优并不一定会实现运载器总体性能最优,以运载器射程最大为目标,对发动机内弹道曲线构型进行优化设计,可显著提高运载器总体性能。研究成果可应用于固体运载器总体方案初步论证阶段,提升固体运载器的性能水平及方案论证效率,使固体运载器总体设计方案更具时效性和竞争性。     

型面旋转变马赫数风洞喷管的优化设计

型面旋转的变马赫数可调风洞喷管因其能够实现出口马赫数的连续快速变化、调节容易和流场品质较好等优点,成为目前地面变马赫数风洞研究的热点。将出口马赫数范围为2.0~4.0的变马赫数风洞喷管作为研究对象,基于iSIGHT和Fluent等软件搭建了适于该风洞喷管的优化设计平台,研究了型线设计中设计点马赫数、初始膨胀线和附面层修正量等因素对变马赫数出口流场的影响,应用NSGAⅡ算法在依据试验设计点所构建的Kriging近似模型上搜索,得到最优的风洞喷管基准型线,并对该基准型线进行二维、三维数值仿真校核,结果显示：通过旋转优选得到的基准型线能够得到马赫数范围为2.0~4.0的变马赫数出口流场,且流场均匀性基本满足固定型面、单设计点风洞喷管的国军标（GJB-1179-91）合格指标,为变马赫数可调风洞喷管在变马赫数风洞中的实际应用奠定了良好的基础。     

磁悬浮助推发射气动力分析及风洞试验

通过分析磁悬浮助推发射装置的空气动力特征, 提出了其气动外形设计的合理原则.分析了磁悬浮橇体和支撑分离机构的气动外形方案, 并确定了优化的气动外形设计参数值.考虑地面效应模拟的同时, 在FD-06风洞中进行了缩比模型试验.结果表明, 升阻力系数随着马赫数或运载器迎角的增大而增大, 装置有上仰的趋势, 迎角为0°时, 阻力系数较小, 各系数的变化也较小.      

吸气式激光推进飞行器发动机特性分析

介绍了吸气式脉冲激光爆轰推进的环聚焦构形飞行器概念,建立了进气道流动模型;利用Chapman-Jouguet爆轰理论和Sedov自相似律,建立了冲量发生与推力传递模型,给出了冲量耦合系数及等效比冲与飞行参数、激光脉冲参数之间关系的解析表达式。计算了发动机冲量耦合系数及等效比冲随飞行高度、马赫数、入射激光强度和脉宽等因素变化的曲线,结果可为相关构形飞行器的弹道与激光器工作模式设计提供参考。     

四段修型弥散反射激波中心体基准流场研究

针对高超声速内收缩进气道宽马赫数范围工作的要求,设计了一种四段修型弥散反射激波中心体基准流场,可明显提高基准流场来流马赫数高于设计点来流马赫数（6.0）时的压缩效率,巡航点（来流马赫数为7.0）出口总压恢复系数较下凹圆弧中心体基准流场提高了2.3%.此外,基于两种基准流场分别设计了圆形进口内收缩进气道并在来流马赫数为5.0~8.0范围内进行数值计算,结果表明:来流马赫数高于设计点来流马赫数时,四段修型进气道的压缩效率更高.有黏条件下,来流马赫数为8.0时二者的增压比近似相等,四段修型进气道喉道截面出口总压恢复系数相对提高了1.1%.      

平头增阻再入体俯仰动态特性计算与流动机理分析

采用自由振荡法数值模拟了一种平头增阻再入体的高超声速俯仰动态特性,沿用稳定性导数模型作为动态特性的指标,通过把稳定性导数视为攻角的函数来考虑非线性影响,并应用奇异分解线性最小二乘法辨识稳定性导数。对于本文研究对象,计算结果表明,平头外形在较高马赫数和小攻角范围内存在动不稳定现象,随着再入时马赫数的降低,动不稳定的攻角范围不断缩小,直至演化为动稳定的。文中对平头外形出现这种动态特性现象的空气动力学原因进行了分析,有必要作进一步的确证与研究。     

高超声速轴对称流道冷流特征及气动力特性研究

对一种轴对称形式的高超声速飞行器全流道开展了风洞实验和数值模拟研究, 分析了不同来流总压、飞行攻角全流道的流场结构和气动力特性.研究结果表明:(1)一定范围内雷诺数的变化对全流道的流动结构和模型的气动力特性无显著影响, 因此所获得的风洞实验结果有望通过某种形式推广到飞行状态下使用;(2)飞行攻角对全流道的流动结构和升力系数有着显著影响, 但阻力系数的影响并不明显;(3)研究范围内来流马赫数的变化对全流道的流动结构有着一定影响, 但研究范围内, 阻力系数随马赫数的变化幅度较小;(4)由于轴对称流道的浸润面积较大, 研究范围内该类飞行器的摩擦阻力在全机阻力中占据了较大的比重, 设计状态下达全机气动力的62%;(5)与实验结果的对照表明, 所采用的数值模拟方法具有较高的精度.      

高超声速锥导乘波体非设计点性能研究

对设计马赫数6的锥导乘波体在马赫数4～7、迎角-6°～+6°的三维流场进行了数值模拟.研究表明:粘性对阻力系数的影响较大;非设计马赫数时,锥导乘波体的升阻比没有明显减小;升力系数与迎角呈线性关系,随着迎角增大而增大;升阻比在+2°迎角下达到最大.     

RBCC飞行器爬升段轨迹设计方法

 火箭基组合循环（RBCC）的发动机推力与飞行轨迹相互影响,导致飞行器轨迹设计与发动机性能分析存在耦合作用。对RBCC飞行器爬升段的轨迹设计方法进行研究,提出了基于马赫数-动压参考曲线的轨迹设计方法。对非均匀有理B样条(NURBS)曲线进行了定义补充,以用于描述具有任意形状的马赫数-动压参考曲线,并对参考曲线的各控制参数选取方法进行了研究;建立了基于二分法求解迎角并实现轨迹方程求解的算法流程。利用提出的轨迹设计方法对空中载机发射的RBCC飞行器进行了爬升段轨迹设计与分析,计算结果表明：(1)马赫数-动压参考曲线法考虑了发动机性能与飞行轨迹的耦合作用,能够适用于RBCC飞行器爬升段的轨迹设计;(2)引射模态低速段（Ma<2.0）消耗的推进剂质量超过爬升段的50%以上,是发动机性能优化的关键;(3)引射模态推进剂流量最大值与最小值之比达到了6.3;(4)当飞行马赫数达到1.7后引射模态下的来流空气冲压作用超过了一次火箭的引射作用,在保证空气捕获方面占主导地位。     

高超声速内收缩进气道轴对称基准流场改进

针对轴对称基准流场中前缘曲激波靠近中心体的部分激波强度过大现象,基于马赫数分布可控反设计方法,将这道前缘曲激波分解为一道较弱弯曲激波和部分等熵压缩波,改进的基准流场存在“四波四区”结构且压缩效率明显提高.基于该改进的基准流场和常规“两波三区”基准流场分别设计了圆形进口的内收缩进气道并对其流场特点和性能进行数值研究.结果表明:改进的进气道的流场能较好保持基准流场的特点;在来流马赫数为4.0~7.0范围内具有较高的压缩效率和良好的流量捕获能力,设计点的出口压比和总压恢复系数分别为17.56和0.540;改进的进气道采用来流马赫数从高到低前缘弯曲激波和汇集的等熵压缩波依次封口的设计概念,在提高流量捕获能力的同时减小了总压损失,总体性能优于常规进气道,来流马赫数为7.0时总压恢复系数相对提高了23.6%,来流马赫数为4.0时流量系数相对提高了5.7%.      

曲面形栅格翼气动特性研究

与常规栅格翼布局相比，曲面形栅格翼布局减小了栅格翼翼元的等固壁通道的长度，从而有效降低了翼元内气流的壅塞，另外，这种栅格翼布局结构简单，易于折叠，减小了飞行器的轮廓尺寸，具有很好的工程应用前景。本文通过数值模拟方法，在亚、跨、超声速条件下，研究了曲面形栅格翼布局的减阻效果和减阻机理，对比了曲面形栅格翼和常规后掠形栅格翼布局的减阻效果，并对不同迎风方式的曲面形栅格翼气动特性进行了分析。     

RBCC进气道双流道变几何方案研究

提出了一种基于双流道变几何的二元进气道设计方案:通过下流道的开关,实现进气道压缩角和收缩比的调节,使得进气道能够适应更宽的工作范围,并能兼顾与飞行器和燃烧室的一体化设计.研究表明:该方案可在保证原有双模态冲压发动机进气道高来流马赫数性能不变的前提下,将自起动马赫数降至引射/亚燃模态过渡点附近,该方案自起动马赫数为2.8;低来流马赫数时进气道具有较高的流量系数（来流马赫数为3.0时达到0.62）、较高的总压恢复系数（来流马赫数为3.0时达到0.79）和较低的阻力（低来流马赫数时降低了30%多）.      

直接的矩阵分裂计算跨音速喷管流场

给出了分离系数矩阵求解欧拉方程组时系数矩阵的直接分离公式，减少和简化了应用该方法时的矩阵运算，参照Ｂｅａｍ－Ｗａｒｍｉｎｇ显式格式构造与ＳＣＭ方法一致的内点差分格式，并分析了它的数值特性，对有斜激波的跨音速内流场进行了数值仿真。计算中采用变系数的当地时间步长，加快了流场达到的收敛速度，计算无需引入任何人工参数，计算值与实验结果吻合。     

槽壁几何参数对跨声速风洞流场品质的影响

采用数值模拟方法研究了再入片开角、槽腔扩开角和开槽形状对跨声速槽壁风洞流场品质的影响。在数值模拟方面,介绍了对应的出入口边界条件设置方法,通过量纲分析得出适用于内流问题的相似参数。研究结果表明:再入片关闭工况下,其前缘处形成了声速喉道,难以建立稳定的试验段超声速流场;再入片打开工况下,试验段流场的不利扰动主要来自于从驻室到试验段的槽壁出流现象,其根源是驻室涡系产生的垂直于槽壁方向的冲击气流,较大的槽腔扩开角对冲击气流有一定的阻挡作用。槽壁各部件对试验段流场品质的影响均是通过增强或抑制槽壁出流来实现的,其中开槽平面形状的影响尤为突出。采用基于自由变形(FFD)的优化设计方法得到了适用于试验段马赫数0.8到1.2状态的开槽形状,优化结果显示初始构型中存在的槽壁出流大部分得到消除,模型区马赫数均方差降低了一个量级,证明了该方法的合理性与实用性。     

双星预警系统弹道参数估计及一些问题分析

主要讨论由2颗地球同步轨道卫星组成的卫星预警系统(简称双星预警系统)中弹道导弹参数的估计问题。阐述了主动段运动的建模问题和弹道数据处理方法,考虑到微分求速的重要性,采用了带光滑约束的样条拟合方法。结合预警系统对弹道导弹预警的具体情况,分别讨论了预警系统下的发射原点估计、射向估计、关机时刻估计、关机点弹道参数估计、落点估计等问题。还讨论了系统中存在的两个问题的分析与解决。模型及分析方法均经过理论和仿真的验证,可应用于双星预警系统下的实际预警的研究。