超声速飞行器侧向喷流干扰流场传统数值模拟方法的误差分析

使用CFD方法,分别就真实喷管边界和简化喷口边界,计算超声速飞行器侧向喷流干扰流场,研究边界条件对干扰流场及气动力的影响.使用k-ε湍流模型封闭雷诺平均N-S方程,利用非结构网格对流场进行空间离散.通过对比,计算结果与实验值吻合良好,证明该方法具有一定可靠性.进一步研究表明喷流边界条件对喷流干扰流场具有一定影响:相对于简化喷口边界,真实喷管边界喷流出口的非均匀性导致喷口上游分离涡和激波位置较为靠前,从而引起附加气动力和力矩的变化;由于摩擦阻力的作用,真实喷管静推力存在损失;喷流压比为500时,总法向力和总俯仰力矩在两种边界条件之间的误差分别为8.21%和22.4%,误差较大.在进行侧向喷流干扰流场的精确计算时,需要考虑边界条件的影响.      

全太阳系深空探测轨迹全局优化(英文)

针对太阳系中全部的248997颗行星的探测问题,给出了一种关于探测飞行器的深空探测全局四维轨迹(t,x,y,z)优化方案,即飞行器从地球发射进入太阳系并采用小推力控制,优化方案的性能指标为飞行器与太阳系中全部行星中相遇和交会的星的数量最多并且燃料消耗最少。本方案给出了四维飞行轨迹进行全局优化的一套算法,该算法由搜索算法和四维轨迹优化算法组成。此搜索算法从太阳系的248997颗行星中寻找获得尽可能多的经过近地球3维走廊内的行星;而四维轨迹优化算法由改进的动态规划算法、基于最优控制理论的共轭梯度算法和静态参数优化算法组成,其中静态参数优化算法用于搜索最优发射时间窗口。基于该组合算法,通过长时间的大规模的飞行数字仿真,最终计算出探测器的四维最优飞行轨迹,在一年内路过了太阳系中全部行星中的12颗行星。     

考虑弹性转轴的扑翼三维气动数值计算

通过求解雷诺平均Navier-Stokes方程,研究了同时具有扑动和俯仰运动的三维扑翼气动特性.为了模拟扑翼的弹性特性,在扑翼俯仰轴处设计了扭簧.研究了扑翼平面形状、俯仰弹性轴与质心的相对位置、俯仰刚度变化等参数的影响.结果表明:俯仰角相对扑动角有一个大的相位超前量,且刚度系数越大该超前量越大;考虑了俯仰弹性特征后,计...     

弹性飞行器大机动飞行气动特性分析

针对两种不同气动布局的前、后掠翼无人机,采用CFD/CSD松耦合方法,对弹性飞行器承受过载为10和20情况下的气动特性与刚性飞行器进行了对比分析。计算结果表明,与刚性飞行器假设情况相比,考虑过载作用下的弹性前掠翼飞机升力减小、阻力减小;而弹性后掠翼飞机的升力增大、阻力增大;两种气动布局的弹性飞行器与刚性情况相比,升阻比都略微增大,而且俯仰力矩系数明显增大。这对于飞行器纵向气动操纵特性会产生较大影响,同时也说明在大机动飞行过程中计入飞行器的弹性变形是非常必要的。     

一种新型垂尾抖振抑制方法实验研究

现代高性能三角翼/双垂尾布局战斗机的垂尾结构普遍受到严重的非定常抖振载荷的困扰。根据自诱导理论提出了一种新型的垂尾抖振抑制方法,利用机头处的静态或振动式硬质鼓包,使三角翼前缘涡涡核弯曲、扭转,从而改变前缘涡的轨迹,延缓涡的破裂,减弱前缘涡破裂尾迹在垂尾周围流场处的脉动强度,以达到抑制垂尾抖振的目的。在西北工业大学低湍流度风洞实验室进行了风洞实验,实验所用模型为一个铝制的全机模型,该模型由一个70°大后掠的三角翼,以及两个31°后掠的垂尾组成。风洞内实验段的风速为10m/s 以及20m/s,迎角范围为20°~50°。实验目的是测量机头处的静态或振动式球形鼓包对垂尾抖振的抑制效果。在尾翼根部两侧粘贴有半桥连接的应变片,用以测量尾翼根部的应变,以此应变作为尾翼抖振强度的衡量标准。实验结果表明,不论是静态的还是振动式的鼓包都不同程度地减缓垂尾的抖振响应,振动式鼓包对垂尾的抖振抑制效果与鼓包的振动频率有关。某一侧的鼓包仅对该侧的垂尾抖振有抑制效果,它不影响另一侧垂尾的抖振响应。频谱分析的结果表明,鼓包在抑制垂尾抖振的同时并没有改变垂尾振动的主频。     

大迎角下翼型振动的气动性能研究

众所周知,风洞试验中的飞机模型,尤其是带有大展弦比机翼的模型有时会出现翼梢振动现象,振动模式主要表现为翼梢沉浮和俯仰形式,以致影响实验结果的精度和可靠性.选取相对厚度较小的NACA0008翼型,在求解非定常Navier-Stokes方程的基础上,采用改进的无限插值理论和绕翼型的C型网格,模拟风洞实验中模型振动条件下的流场,研究振动模式及其不同耦合对流场、尤其是大迎角流场的影响,并考虑了模型弹性轴不同位置对结果的影响.研究结果表明:在临近传统定常失速迎角的大迎角条件下,翼型的振动可以引起翼型大尺度的分离,导致翼型失速的提前发生;振动在不同的相位滞后条件下,对翼型流场的分离程度不尽相同.     

半模机翼振动对气动性能影响的风洞试验研究

为了验证翼型振动测压试验结果和有关振动数值模拟文献的计算结果,设计完成了对称翼型NA-CA0012和层流翼型NACA64-210两种半模机翼的低速风洞振动试验.通过专门研制的半模振动模型,选取5种激振方式,用直接测力法得到了模型在静态和不同激振方式下升力特性变化曲线,分别研究了振动频率、雷诺数、自然转捩和固定转捩、数据采集方式等参数对机翼气动特性的影响.试验结果表明:半模机翼振动的气动特性与二元翼型和数值计算情况有所不同,其影响因不同翼型构型、不同采集方式和翼面不同流动模式等会产生不同效果.并对产生原因进行初步探讨.     

一种新型颤振激励系统的数值计算研究

针对一种新型颤振激励系统的特点,采用非结构旋转运动网格技术,通过数值求解非定常欧拉方程,计算了该激励系统工作时的非定常气动力特性,并与准定常计算结果进行了比较。研究结果表明:采用非定常方法所得结果与准定常计算结果有明显的差别;气动力的变化与激励系统中圆柱的旋转方向有关,旋转角速度越大,气动力变化的幅值越小;气动力的幅值与圆柱缝道宽度也有一定关系。     

风洞实验中机翼振动对气动特性影响的计算分析

通过求解非定常N-S方程研究了弹性机翼在低速气流中振动时的气动特性,并与刚性机翼气动特性做了比较。通过结构有限元方法计算出模型机翼的结构模态和相应振动频率,假设机翼以最小刚度的弯曲模态振动,计算了不同振幅下不同迎角、尤其是失速迎角附近的升力特性。计算结果表明,在小振幅下机翼的平均升力特性和刚性机翼相差很小,但当振幅达到一定幅度后,振动机翼下平均升力特性与刚性机翼的结果有较明显的差别,主要表现在最大升力系数的降低和失速迎角的减小。此结果可以解释实际风洞实验中最大升力系数往往小于CFD结果、实验中最大升力系数对应的迎角也比CFD结果小的现象。     

三维非结构自适应多重网格技术

发展了一种基于几何外形的非结构自适应多重网格技术。根据流场参数的变化梯度确定加密边,由加密准则对粗网格进行自适应剖分得到分布合理的较密网格。通过预先生成的初始极密表面网格来将边界网格的加密点投影到边界上,使得边界保持初始外形。为了提高流场求解速度,将自适应加密后的一系列网格作为多重网格,采用了多重网格算法进行计算。M6机翼的流场计算表明该方法可以准确捕捉激波位置,并且表面的压力系数分布与实验值相吻合。