考虑地面效应的翼型动态特性数值模拟

为了研究考虑地面干扰的飞行器非定常气动特性,基于滑移网格技术,通过求解二维非定常Euler方程,对不同近地高度下的NACA0012翼型的动态振动特性进行了数值模拟。分析了存在地面干扰时的升力系数以及俯仰力矩系数迟滞环变化,并结合强迫振动法,计算了翼型各个近地高度下的纵向组合动导数。结果表明,地面效应不仅对定常流场产生影响,更显著地影响了非定常气动力及力矩,近地高度越小,这种效应越明显,使得升力系数的迟滞环面积变小,幅值减小,而俯仰力矩系数的迟滞环变化不规律,力矩系数的动导数随着近地高度减小而增大(绝对值减小),纵向的动态稳定性产生损失。地面效应干扰在飞行器非定常气动研究中应该引起重视。     

滚转和锥进运动对弹箭动导数求解的影响

现有的弹箭动稳定性导数计算通常忽略了耦合运动效应的影响,而弹箭在实际飞行中,通常同时存在俯仰、滚转以及锥进运动。为此,本文提出了基于欧拉转动定理和滑移网格技术的复杂角运动模拟方法,该方法利用著名的罗德里格斯转换矩阵插值求得弹箭在每个时间步的角速度修正值,并指定给球形滑移网格区,从而模拟弹箭俯仰滚转锥进相耦合的三自由度耦合角运动,在此基础上,通过对非定常气动参数进行求解辨识,分析了滚转频率和锥进频率对其俯仰组合动导数和升力系数动导数的影响规律。结果表明:本文提出的复杂角运动模拟方法,可有效消除姿态角计算的累积误差,实现对弹箭任意给定角运动的准确模拟;在耦合角运动状态下,弹箭气动力系数迟滞环发生明显的振荡和偏移,俯仰组合动导数和升力系数动导数均随滚转频率和锥进频率增加而显著变化,且变化规律具有非线性特性。     

考虑多种物理效应的钝锥俯仰稳定性参数影响分析

首先采用二阶ROE格式计算了非定常NS方程下高超声速钝锥体的俯仰阻尼导数,通过与文献中内伏牛顿法和实验结果对比验证了本文所采用非定常计算方法的正确性.同时采用二阶ROE格式及双时间步长、刚性旋转动网格技术等,数值分析了NS方程计算条件下高空高马赫数钝锥体模型的强迫振动非定常流场特征,并运用最小二乘法辨识出其俯仰静稳定性导数与俯仰阻尼导数.通过计算结果对比分析了不同高度(含稀薄气体效应)、Maxwell滑移边界条件及五组元化学非平衡模型对钝锥体模型非定常流动特征及静、动导数的影响.结果分析表明,本文所计算的钝锥模型在所假设强迫振动下的俯仰稳定性参数随高度变化剧烈,且考虑滑移边界条件与化学反应会使得俯仰力矩迟滞曲线形态发生较大程度改变.准确计算高空、高马赫数条件下飞行器俯仰阻尼导数应考虑化学非平衡效应和稀薄气体效应的影响,是否考虑这些物理效应会关系到飞行器静、动稳定的程度.     

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

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

基于刚性动网格技术的动导数数值模拟

基于刚性动网格技术,建立了俯仰动导数的非定常数值计算方法.首先使用小幅度强迫俯仰振荡方法求解俯仰组合动导数,然后利用小幅度强迫升沉振荡方法求解洗流时差导数,通过两者相减即可得到俯仰阻尼导数.利用国际动导数标准模型Finner导弹进行算例验证,计算得到的俯仰组合动导数与试验值误差为2.76%,洗流时差导数值约为俯仰阻尼导数的11.5%,与文献的结果一致.结论表明:动导数单独模拟方法具有较好的工程实用价值,且可以推广到横向以及航向的动导数数值模拟.      

翼型纵向俯仰非定常气动力研究

飞行器飞行时受到的非定常气动力是飞行器设计以及机动飞行中须考虑的一个重要因素。采用CFD方法数值模拟了NACA0012翼型在马赫数从0.2~5.0范围变化时做强迫振荡的非定常流场。通过积分法计算了在不同马赫数时翼型的俯仰动导数。计算结果表明,迟滞效应随马赫数的增大先增大后减小,在跨声速区域迟滞效应达到最大;迟滞环的绕向在跨声速区域为顺时针方向;迟滞环顺时针绕向时,动导数大于零,模型动不稳定,逆时针绕向时,动导数小于零,模型动稳定;旋转角速度越大,翼型表面非定常气动效应越明显。     

基于CFD数值仿真技术的飞行器动导数计算

为了研究表征飞行器动态稳定特性的动导数,结合计算流体力学方法,基于滑移网格技术,采用小幅度强迫振动法和差分法两种非定常方法,推导了俯仰组合动导数以及滚转阻尼导数的计算表达式,并结合多参考系模型采用准定常方法计算滚转阻尼导数。建立了基于计算流体力学技术的飞行器动导数计算方法。对国际动导数标模Finner导弹进行了计算验证。两种非定常方法对俯仰组合动导数的计算误差分别为0.65%、6.13%,对滚转阻尼导数的计算误差分别为2.5%、0.06%。求解滚转阻尼的准定常方法对于滚转阻尼导数的计算误差为2.67%。三种方法的计算结果与风洞试验结果吻合的很好。非定常方法计算精度高,准定常方法可以快捷迅速得到结果。本文的数值方法不仅适用于超声速,同样可以用于亚、跨声速时的动导数计算。     

某型导弹动导数风洞试验和数值模拟

动导数的精确预测对于弹箭类武器的成功研制至关重要。以某型导弹动态特性的研究为对象,在FD06风洞中采用自由振动法完成动态特性试验,获取该型导弹俯仰和滚转两个方向的典型振幅衰减曲线和动导数结果;针对该型导弹,基于双时间推进法数值求解非定常N-S方程和强迫运动方程,获取该型导弹俯仰动导数和滚转动导数。结果表明:CFD计算结果与风洞试验结果有较好的一致性;在给定马赫数和攻角下CFD数值模拟结果可以为后续工程型号研制提供参考。     

扑翼非定常运动对平尾影响及俯仰力矩特性研究

为了研究微型扑翼飞行器尾流对其平尾设计、飞行器稳定性以及飞行控制的影响,选取微型扑翼飞行器ASN-211为原模型,采用将其简化为二维的前后串列翼模型进行具体计算和分析。首先以Fluent动网格技术为背景,在用户自定义函数控制扑翼的非定常运动条件下进行不可压、非定常二维流动的计算,并研究在扑翼非定常运动条件下的模型的俯仰力矩特性。然后通过计算不同来流攻角、扑翼扑动频率、扑翼与平尾间距以及不同力矩中心下扑翼、平尾及总的力矩系数,讨论各个参数对力矩特性的影响。最后得出不同的扑翼扑动频率以及扑翼与平尾间距将会对平尾与扑翼的俯仰力矩间的相位差产生影响,所得结论为扑翼飞机的重心布置设计提供参考。     

钝锥在添质作用下的稳定性导数的分析

本文在国外有关工作的基础上进一步分析了烧蚀引起的添质作用的力学模型及其对球钝锥俯仰力矩稳定性导数的影响机理,一起计入了热—化学延迟效应和气流延迟效应,考虑了边界层为层流以及转捩为湍流的各种情况,探讨了稳定性导数与添质率、边界层转捩位置、热—化学延迟相位角、减缩频率以及攻角等诸因素的关系。并且提供了计算具有添质作用的球钝锥在小攻角下以及在不同边界层流态下的俯仰力矩稳定性导数的一套方法,较前扩展了应用范围。分析表明,当添质率较大,而在重心附近又发生边界层转捩时,钝锥的俯仰阻尼力矩导数可以出现严重的负阻尼情况。本文的计算结果与某些弹道靶试验数据以及飞行试验遥测换算数据相当接近。     

翼身融合布局缩比模型飞行试验动力影响分析

动力系统的进排气效应是影响翼身融合布局飞机俯仰力矩特性和全机纵向配平的重要因素。为了研究动力对全机俯仰力矩的影响,基于飞行动力学方程,提出了一种动力影响计算方法。以某翼身融合布局缩比模型为研究对象,结合风洞试验数据,分析了模型飞行试验中爬升段、平飞段和下降段动力对全机俯仰力矩系数的影响,同时计算了迎角测量误差对俯仰力矩计算结果的影响。分析结果表明:动力使该翼身融合布局缩比模型低头力矩增加,会对全机纵向配平产生不利影响;迎角测量精度对俯仰力矩计算结果有较大影响。     

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

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

稳定性导数计算方法研究

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

背撑发动机进排气对翼身融合民机气动特性影响

针对翼身融合(BWB)背撑发动机布局民机的发动机与机体流动干扰问题,以项目团队研究的某BWB布局民机方案为研究对象,采用CFD方法研究了发动机进排气对全机低、高速气动特性的影响规律及其流动机理。结果表明:进排气效应主要影响飞机的高速阻力特性,随着发动机进气流量的减小,机体/短舱之间流动干扰加剧,机体零升阻力增加；进排气效应对低速状态升力和俯仰力矩的影响主要由发动机的抽吸作用对机体表面流动的干扰引起。发动机不同的进气流量会对全机的升阻力特性产生较大影响,在BWB布局机体/发动机一体化设计中,要充分考虑机体与动力装置之间的干扰影响,同时权衡不同进排气状态对机体和发动机当地流动的影响。      

弹性变形对柔性飞艇气动特性的影响

设计加工具有相同几何外形的刚性飞艇模型和可变内压的柔性飞艇模型,并进行了低速风洞实验。研究了内压变化对柔性模型阻力和俯仰力矩特性的影响,比较了柔性模型和刚性模型气动特性的差异。实验结果表明,随着内压的增加柔性模型的阻力呈下降趋势;对于具有相同几何外形的柔性和刚性模型,柔性模型的阻力明显大于刚性模型;刚性模型的尾翼对俯仰力矩的影响较大。     

基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.      

两种有翼飞行器高超声速动态气动特性的对比实验研究

采用模型自由飞技术在脉冲型高超声速风洞中测量了两种类航天飞机外形模型的腐爷阻尼导数，两种模型具有十分接近的外形特征尺寸和投影面积，但机身和机翼的剖面则片此各不相同，实验在名义马赫数M∞＝6.4条件下进行，同一名义实验条件下的重复实验显示一致的运动形态和接近的动导数测量结果。气动参数辨识采用最大似然法，对风沿实验准定常试验时间中模型的平面运动以线性气动参数模型辨识得到它们的俯3仰阻尼导数。结果提产两种外形有差异的模型呈现迥然不同的动态气动特性：带OMS舱的航天飞机仿真模型具有动态稳定性，而简化外形的类航天飞机模型则为动不稳定，虽然对导致这种极大差异的直接物理原因还有待深入研究，但实验揭示了动稳定性对模型外形细节的敏感性。     

双旋翼微型飞行器的非定常气动特性

发展了一种基于计算流体力学(CFD)的非定常气动特性预测方法,计算方法包括了动量源模型、预处理方法、非结构嵌套网格和Spalart-Allmaras(S-A)湍流模型等技术.通过计算悬停Caradonna-Tung算例和俯仰振荡NACA0012算例,验证了计算方法模拟双旋翼微型飞行器动态流场的有效性.数值模拟了双旋翼微型飞行器动态流场,给出了非定常气动系数的迟滞曲线,分析了缩减频率、前飞速度和螺旋桨转速对非定常气动特性的影响.计算结果表明:力矩系数迟滞效应随缩减频率、前飞速度和螺旋桨转速增大而增大,升力系数迟滞效应随缩减频率和前飞速度增大而增大,但基本不随螺旋桨转速变化而变化.      

带连续变弯度后缘操纵面机翼的动态失速减缓

主动偏转后缘操纵面可以减小翼型动态失速对气动特性产生的不利影响。研究连续变弯度后缘操纵面在减缓翼型动态失速方面的性能,利用CFD 结合动网格方法,计算NACA 0012 翼型在大幅度俯仰振荡时的非定常气动力;从减缓效果和能量需求两个方面,对比传统刚性操纵面和两种连续变弯度操纵面的动态失速减缓性能。结果表明：当后缘操纵面按正弦脉冲规律偏转时,可以推迟前缘涡的产生,加速后缘涡的发展,降低压力分布在后缘的峰值,进而减小动态失速时翼型的低头力矩极值;后缘操纵面的弯度构型会影响减缓效果,在相同的偏转策略下,弯度描述函数为2 阶多项式的连续变弯度操纵面的减缓效果最好,且能量需求最小。