四元数微分方程最小参数解法研究

对一种求解正交矩阵微分方程的三阶最小参数方法进行改造 ,将其用于四元数微分方程求解 ,并对其在捷联惯性导航系统姿态矩阵实时求解中的应用进行了仿真。结果表明 ,其计算精度比经典的三阶Peano Baker数值解法提高 35倍 ,而浮点运算次数则仅为后者的一半 ,对于提高捷联惯性导航系统姿态计算的精度具有重要意义。     

超声速流动现象的水法流模拟与仿真

根据流体流动控制方程组,研究发现明渠中的浅水波和气体中的扰动波具有相同的传播特征,因而可用直观形象的水流实验来模拟复杂的气体动力学现象.运用收缩扩张形实验明渠得到的水流实验结果与气体流场数值仿真结果具有良好的一致性,验证了明渠流模拟气体动力学现象的有效性.研究结果表明:拉瓦尔喷管中的马赫数M与同尺寸明渠中的弗劳德数Fr具有相同的变化规律;明渠中的超临界流通过楔形体和圆柱体时的流动现象可以比较理想地模拟出超声速气流在相应条件下的膨胀波和激波现象;通过改变出口外部水位的方式在明渠内形成的水跃波,可以非常形象地模拟出拉瓦尔喷管在相应条件下形成的管内激波.     

四元数、Rodrigues参数在卫星姿态解算上的对比研究

详细阐述了罗德里格(Rodrigues)参数及改进型罗德里格(Modified Rodrigues)参数与四元数间的关系;推导了Rodrigues参数及改进型罗德里格参数微分方程和姿态阵,提出了后者避免奇异性的方法。通过仿真对以Rodrigues参数及改进型罗德里格参数描述的卫星姿态算法和以四元数描述的姿态算法在计算量和计算精度上进行了详细比较。仿真结果表明,Rodrigues参数和改进型罗德里格参数计算精度相当且计算量上优于四元数。     

编队飞行航天器平均轨道根数非线性控制研究

基于非线性动力学模型,研究了环绕编队航天器的相对运动控制问题,介绍了一种基于平均轨道根数的航天器编队飞行非线性闭环控制方法。以航天器在惯性坐标系下位置速度向量误差为输入量,给出了一种全状态反馈控制律,利用Lyapunov方法证明该闭环系统是渐进稳定的;最后利用该方法对一个编队进行控制,仿真结果表明了该控制方法的可行性。     

高超声速平板近空间气动特性的计算分析研究

采用DSMC方法对二维小尺度平板在近连续流至自由分子流区域、攻角为0°～30°下的高超声速气动特性进行计算分析,旨在探索高超声速飞行器在近空间空域内气动特性的变化规律和影响因素.结果表明,升阻比随努森数、马赫数的增加单调下降,相对比来说,马赫数对升阻比的影响不是很明显,而努森数的影响比较明显.特别针对摩阻特性进行了分析,在马赫数10、高度55～110km、0°攻角时摩阻在总阻力中所占比重高达60%以上.随着攻角的增加,由于波阻的急剧增加导致摩阻所占比重下降.在10°攻角下,摩阻在总阻力中所占比重约为45%～70%.30°攻角时摩阻在总阻力中所占比重则下降为12%～50%.同时还发现,在近自由分子流区域中,0°攻角时阻力系数随着马赫数的增加而下降,在有攻角时该马赫数效应呈现减弱趋势.     

四元数反馈控制技术研究

弹道导弹、运载火箭在飞行过程中经常会出现三通道大范围调姿的情况,如果用传统的姿态角控制方案,当3个欧拉角中的第2个欧拉角出现过90°的情况时,欧拉角计算出现奇异.本文针对采用捷联惯性测量组合作为敏感器件的弹道导弹、运载火箭提出了四元数控制方案,即采用初始四元数q0袁征初始姿态、终端四元数q1表征终端姿态,实时计算调姿程序四元数qcx以及姿控系统控制量四元数△q(k)的方法,可避免姿态角计算奇异问题,且姿态过程优化.数学仿真证明了该方案的正确性和可行性.     

前体非对称涡流动临界雷诺数效应及分区特性

通过模型表面测压和油流显示,对旋成体于50°迎角在临界雷诺数区域(0.13×106~0.81×106)的压力分布和侧向力特性随雷诺数变化的演化规律进行了研究,结果表明,随着Re数从亚临界增加至临界区域,模型表面的低位涡侧首先出现层流分离气泡成为转捩分离(Tr),而高位涡侧仍处于亚临界层流分离(L),非对称更为显著,侧向力较亚临界区有所增加;随着雷诺数进一步增加,高位涡侧才成为转捩分离,此时非对称流动逐渐演变成对称流动,压力分布呈对称的平台状,侧向力明显减小,因此,通过流动分离前的压力恢复值作为判则,根据旋成体两侧边界层分别处于L/Tr和Tr/Tr状态,可将临界雷诺数区域划分为临界起始发展区和临界区。最后据此判则讨论了旋成体绕流沿轴向多种流态共存的现象。     

PIV技术在前体非对称涡结构研究中的应用

介绍了北航D4风洞PIV系统的布置及具体实验方案,在此基础上实现了PIV技术在前体非对称涡流动结构研究中的应用。在迎角50°、Re=0．14×10°～0．55×106时,对旋成体X／D=2和3．35截面流动结构进行研究。结果表明,随着胁数的增加截面上流动结构存在从非对称二涡向三涡发展的趋势;在亚临界区,旋涡对非对称压力分布的影响起主要作用;在临界起始发展区及临界区,边界层流动状态及其分离形态对非对称压力分布的影响起主要作用;前体非对称涡沿轴向由二涡向三涡的发展状态在临界起始发展区比亚临界区将向更上游的位置发生。     

强激波阵面的非平衡特性研究

利用测量强激波波后N_2~+第一负系(0,0)带和(1,2)带的辐射,对强激波后振动温度历程的进行了测量;利用Langmuir探针技术,在低密度激波管中对强激波后电子数密度历程进行了测量;在弹道靶中进行了激波脱体距离的测量研究.测量和计算结果进行了对比,结果表明:N_2~+B~2∑_u~+态的激发比振动能的激发更快:实验测得的振动温度有明显的周期性振荡;在激波速度7.65km/s～7.85km/s、p_1=1.33Pa、实验段内径0.8m下,实验有效时间只有约6.5μs,实验中的电子数密度不能达到峰值,在约10倍波前自由程的实验有效区域内,电子数密度的测量值与计算值吻合很好.激波脱体距离的测量误差约为5%,是目前国内精度最高的结果.     

旋翼翼型低Ma数动态失速特性计算

低Ma数下,翼型前缘涡强度的增加和移动使法向力系数产生很大的超调量,Beddoes通过增加一项与延迟后后缘分离点有关项来模拟该特性,并改进Leishman-Beddoes二维翼型动态失速模型。在此基础上,本文在非定常法向力系数中引入一阶延迟,推迟失速判断点,得到修正后模型;而后,通过计算NACA0012、OA207翼型在低Ma数下的非定常气动力,并与实验结果进行对比,验证了模型在计算翼型低Ma数下非定常气动力的准确性,并分析了折合频率、迎角平均值、振幅对计算结果的影响。