超,高超声速机翼俯仰导数的欧拉方程摄动解法

应用片条理论，将文献（１）的尖前缘任意翼型的俯仰稳定性导数的计算方法作、推广、用于计算任意平面形状的三维机翼。在各种特殊的情况下，其计算结果与现有的理论都符合得很好。在激波附体的条件下，本文的计算方法可用于计算任意迎角。任意平面形状机翼的超、高速声速俯仰稳定性导数。     

基于非定常CFD方法跨声速民用飞机标模俯仰气动导数计算研究

采用非定常计算流体力学(CFD)方法和动网格技术,计算了跨声速民用飞机全机构型TCR气动标准模型在0°~30°迎角范围内气动力系数的时均值、俯仰静导数和俯仰动导数;通过与试验结果对比,研究了非定常CFD方法在较大迎角范围内对包含鸭翼、垂尾等部件的复杂外形气动导数的计算精度。与试验结果对比表明,在16°~24°迎角范围内,受鸭翼和机翼上流动分离再附以及内外翼分离涡相互作用等复杂流动结构的影响,俯仰力矩系数的动导数与试验数据有相对较大的差异。     

尖侧缘机身布局的俯仰力矩特性及扰流板控制

针对尖侧缘机身布局在大迎角下存在的正俯仰力矩（抬头力矩）问题,通过风洞试验,首先研究了俯仰力矩的迎角分区特性及流动演化规律：线性增长区（迎角为0°~15°）,俯仰力矩线性增加,全机从附着流到形成进气道前缘涡和机翼涡;非线性增长区（迎角为17.5°~32.5°）,俯仰力矩非线性增加,机头涡出现,机头涡和进气道前缘涡逐渐增强,机翼涡增强后破裂;衰减区（迎角为35°~65°）,俯仰力矩逐渐减小,机头涡增强后破裂,进气道前缘涡破裂发展,机翼涡完全破裂。其次,发现了机身前体是产生正俯仰力矩的主要来源,机头涡是导致大迎角下正俯仰力矩的主控流动。当迎角为40°时,前体各截面正俯仰力矩在进气道前缘处达到最大,主要是由于该处机头涡诱导产生了较强的法向力。最后,提出了大迎角机身扰流板控制技术,产生了较好的控制效果。当迎角为40°时,扰流板可使正俯仰力矩减少62%,其原因是扰流板降低了机头涡涡量及其诱导产生的法向力,减少了机身前体对正俯仰力矩的贡献。该控制技术的缺点是扰流板会带来一些升力损失和附加阻力。基于尖侧缘机身参考宽度的雷诺数为2.59×10⁵。     

牛顿流理论及其在谐振机翼上的应用

根据牛顿流理论,文中讨论了高马赫数、任意迎角、一般剖面在均匀来流中作沉浮和俯仰谐振的翼型和机翼。导得翼型和机翼表面的牛顿撞击压强和安定性导数公式。给出楔、平板三角翼、双抛物线剖面三角翼俯仰刚性导数(-C_(mθ)和阻尼导数(-C_(m))的数值例子。数值结果与现有的牛顿-布兹曼流理论、薄激波层理论、欧拉方程摄动理论的结果作了比较,算例表明本文的结果令人满意。     

翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。     

计算高亚声速情况下三维机翼动稳定性导数的面元法

本文应用我们提出的跨-亚音速三维谐振机翼的局部线化面元法,同时采用了与之配偶的亚临界定常流绕三维厚度机翼的面元法迭代解,对几种不同平面形状的机翼,进行了俯仰和滚转阻尼导数的计算和分析,并和现有的实验数据以及线化理论的计算结果进行了比较,从而表明,本方法可以得到满意的动导数计算结果;还表明,当接近临界马赫数时,动导数计算必须计及跨音速非线性效应。     

非定常超音速侧滑机翼的升力、力矩系数及动导数的数值计算

 本文利用分布在机翼平面上的非定常超音速点源的特征线网格法给出非定常超音速侧滑机翼的数值解。本文给出侧滑机翼作沉浮、俯仰、滚转谐振时的振型以及计算升力、力矩和动导数（包括交叉导数）的计算方法。按本文用非定常理论数值解对典型的无侧滑机翼计算出的动导数与常用的按准定常理论解析公式计算结果进行比较和鉴定,结果表明本文方法是令人满意的。     

典型超声速/高超声速动导数计算方法研究

以当地流活塞理论为基础完成超声速/高超声速动导数计算公式推导,发展了一种高效的动导数计算方法.选取1个超声速、2个高超声速动导数算例标模,计算了超声速/高超声速下飞行器动导数变化规律.研究结果表明:BFM超声速流动算例的俯仰组合动导数绝对值随马赫数的增大而减小,计算重心位置后移亦会导致俯仰组合动导数绝对值减小,纵向动稳定性降低;在0° ～20°迎角范围内,高超声速流动中的尖锥和钝头旋成体俯仰组合动导数绝对值随迎角的增大而增大,纵向动稳定性增强;所提方法计算时间约为双时间动导数计算方法的1/27.     

边条机翼流态和气动力特性实验研究

本文论述了用油流、测压和测力等实验方法研究得出的边条机翼在低速、跨、超音速情况下的绕流情况和气动力特性;分析了低速时绕边条机翼的附体流、前缘分离和附体流并存、前缘涡和气泡涡并存、旋涡破裂等四种主要流态,在跨、超音速时出现的前缘分离和前缘附着两种主要流动类型,以及它们与气动力特性的关系。文中,结合测压和测力等实验结果,分析了边条机翼气动特性随M数和迎角的变化情况,说明了在中等后掠翼上加边条对于改善大迎角下机翼升阻力特性的作用,并对边条引起的俯仰力矩非线性变化现象,提出了可能的改进方法。     

前缘剖面形状对80°／60°后掠双三角翼上涡流运动影响的数值模拟研究

采用上风高分辨率格式———通量差分分裂格式离散求解三维可压非定常薄层N S方程 ,数值模拟前缘剖面形状对 80°/60°后掠双三角翼上涡流运动的影响。计算模型包括尖前缘、菱形前缘、圆形前缘等三种不同前缘形状的机翼。计算结果表明 :三种不同剖面形状的前缘可以诱导产生不同的前缘分离 ,形成的各前缘剪切层的特点也不同。尖前缘机翼的边条翼涡和外翼涡合并点最靠前 ,其次为菱形前缘和圆形前缘的机翼。在大迎角情况下 ,三种机翼上的内翼涡发生破裂。圆形前缘机翼上的内翼涡破裂点比其它两种情况下破裂点的位置较靠前。涡的合并对二次涡的结构和特点也有显著的影响。     

机翼动态气动特性实验研究

介绍了动态特性风洞实验的技术方案、测试系统和测试软件。对前缘后掠６０°的三角翼在迎角０°～９０°范围内快速上仰时的升力进行了测量。结果表明：机翼快速上仰时,随上仰速率增大最大升力系数增高,即动态升力效应更明显;失速迎角也增大。此外,给出了简缩频率、平均迎角、振荡振幅以及俯仰轴位置等对振荡机翼产生的非定常瞬时力的影响     

计算谐振机翼气动力的激波膨胀波级数展开法

本文根据激波-膨胀波理论,忽略了反射波对翼面压强分布的影响,把斜激波之后和通过一族主特征线的压强展成级数,给出了沉浮和俯仰谐振翼型表面压强分布的近似解析解。然后应用片条理论导得三维机翼振动导数的闭式解。文中给出三角形机翼俯仰导数的几个实例,数值结果与现有其它方法的结果有很好的一致性。     

洗流时差阻尼力矩研究

在分析飞机的纵向阻尼力矩时，除考虑俯仰角速度引起的俯仰阻尼力矩外，还应估算出迎角变化率引起的洗流时差阻尼力矩。根据洗流时差阻尼力矩产生的原理，推导并给出了它的计算方法。介绍了各种阻尼导数的表示方法及相互关系。最后较详细地分析了水平尾翼、马赫数、迎角变化率和下洗角随迎角的变化率诸因素对洗流时差阻尼力矩及阻尼导数的影响。     

弹性飞机纵向机动飞行轨迹计算和分析

考虑气动弹性静变形计算某飞机铅垂纵向平面内的机动飞行轨迹。将飞机当作弹性体，在运用飞行动力学方程计算运动轨迹的同时，针对每个时刻，对飞机的迎角升力线斜率、升降舵升力线斜率、俯仰角速度升力线斜率、迎角俯仰力矩导数、升降舵操纵导数以及俯仰角速度俯仰力矩导数进行气动弹性修正，给升降舵一个阶跃输入，计算出弹性飞行情况的飞行轨迹，并与刚体飞行情况对比，从而说明气动弹性对飞机飞行轨迹的影响。     

升力浮力复合型飞艇动导数分析

计算了球体和椭球体的平移、转动加速度导数,数值模拟了NACA0015翼型的俯仰振荡,验证了计算流体动力学(CFD)方法计算动导数的有效性.模拟升力浮力复合型飞艇的强迫振荡运动,得到了全艇以及机翼、艇身、平尾和垂尾的动导数,并与常规飞艇进行了比较分析.机翼对速度导数的贡献导致复合型飞艇升沉速度导数、滚转角速度导数较常规飞艇显著增大.常规飞艇加速度导数以惯性力成分为主,但黏性力成分以及机翼对艇身和尾翼的干扰对复合型飞艇升沉、滚转和俯仰3个方向的加速度导数有重要影响,不能将加速度导数简化为附加质量系数.      

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

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

细长翼涡破裂对气动特性影响的计算

本文在简化的前缘涡模型基础上,沿涡轴分布点源,以模拟细长翼涡破裂对气动特性的影响。对一组不同平面形状细长翼的计算表明,用本方法得到的结果在定性上与机翼在失速迎角前后实际的气动特性变化趋势是一致的。     

机翼大迎角涡破裂时气动特性计算方法的研究

本文计算了小展弦比机翼纵向和横侧气动特性,并计及涡破裂的影响。本方法以Purvis升力面理论和Polhamus吸力比拟为基础。并参照Lan方法进行大迎角特性修正。算例表明,本方法适用于小展弦比机翼气动特性的计算。     

机翼大迎角气动力的快速算法

在定常势流 Green函数方法计算位流气动力的基础上 ,应用前缘吸力比拟法计算机翼大迎角气动力 ;并考虑了涡破裂所产生的影响。计算结果与实验结果吻合良好 ,为高机动飞机初步设计提供了一种快捷、准确的气动力计算方法     

用准涡格法计算前缘分离的细长三角翼

大后掠细长翼在大迎角下产生的前缘分离涡,对气动特性具有重要影响。计算前缘涡气动特性的理论方法已有多种。较早的有Smith的方法,它基于细长体假设,故不能满足后缘条件。前缘吸力比拟法及其推广,对于计算整个机翼的力和力矩来说是一种既简便又比较准确的方法,但给不出详细的载荷分布。文献[3]把涡格法推广引用于有前缘分离的机翼,用若干离散的涡线代表有涡核存在的实际前缘涡,算出的力和力