一种改进的计算压力分布的格林函数方法

对计算压力分布的格林函数方法作了一定的改进,虽然格林函数方法能给出比较准确的速势分布计算,但在压力分布的计算方面却比较欠缺,只给出了一个根据离散点上的速势值进行有限差分计算的简单的公式,采用这个公式无法算出物面准确的压力分布。鉴于此,通过对速势积分方程离散化后所得到的扰动速势的解析表达式进行微分运算,得到一组计算扰动速势导数的解析公式,并由此计算了物面上各离散点上的压力系数,其结果相当准确,与Navier-Stokes方程算法及实验结果相比,均吻合较好。另外,还对库塔条件在数值求解速势积分方程中的应用作了一定的改进,使计算结果更加准确。     

二元全速势跨音速方程的数值解

本报告提供了二元全速势跨音速方程的一种数值解法。本方法的特点是放弃了繁杂的曲线座标系,而在简洁直观的直角座标系中,利用Jameson的旋转差分格式,并在物面边界上采用速势内嵌场,建立全流场的差分方程。二元翼型的样例计算表明本方法和其他全速势方程的数值解法具有同等的精度。对物面外形作一次附面层位移厚度修正后的算例表明计算结果和试验结果符合较好。 本方法比起采用曲线座标系的方法,更易推广到三元流。     

用二维N—S方程计算汽车流场的地面效应

采用有限体积方法求解了Navier-Stokes方程,计算了绕二维类汽车物体的绕流并研究了地面效应。在非正交网格系统上以强守恒形式构造出有限体积方程,采用压力校正算法来求解所得到的方程。对于湍流情况,采用了双方程模型的κ_ε模型。 为了验证采用本程序计算移动壁面效应的正确性,本文给出了Couette流动的计算结果。 本文也给出了汽车绕流的计算结果,它包括层流情况,从汽车头部开始由层流转为湍流的混合流动情况及湍流情况。对于每一种情况都比较了移动壁面的计算结果和固定壁面的计算结果。     

机翼机身跨声速绕流的计算

本文采用变换与数值网格生成方法，对由椭圆抛物而描述的每一个展向曲面进行网格生成，引入一种用于控制网格疏密分布的函数变换方法，其特点是易于对不同族网格线分别控制，根据物面梯度自适应地调整网格的疏密。用有限体积法计算守恒型全速势方程跨声速绕流，对M-100翼身组合体进行计算。理论计算结果与实验结果表明，本文方法有效。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

亚、跨音速机翼携带外挂物气动力干扰的混合差分计算

本文应用跨音速、定常、小扰动势流的混合差分方法,计算了机翼—挂架—外挂物气动力干扰。利用线超松弛改进迭代方法,在物理空间网格点上满足x向大扰动速势方程,在物面上满足精确的边界条件,在涡面上满足库塔条件,远场处取速势方程的线性解。计算网格点上的速势值ψ、下洗速度ψ_y、侧洗速度ψ_z:的分布,以及所有部件、组合体的压强分布、气动力系数,及机翼、外挂物各自所受的气动干扰量。 本程序用BCY语言在上海华东计算技术研究所655机上进行计算。文内三个算例均得到收敛或接近收敛的结果。与可以找到的风洞实验比较尚一致。     

旋转体跨音速零升波阻和压力的差分计算

应用Murman-Cole的有限差分法,求解具有纵向大扰动而横向小扰动的跨音速轴对称速势方程,由此计算旋转体跨音速零升力时的压力分布和波阻力,以及激波位置。物面边界条件被转移到物体轴上。远场边界条件由无穷远处的条件近似代替。计算物面压力系数时,用细长体理论进行物面速势插值。 速势的差分方程用沿半径方向线超松弛改进迭代求解。网格取62×16,迭代初场取零,达到收敛的迭代次数对M_∞<1,M_∞>1以及M_∞≈1分别大约为150,40和300次。松弛因子取为:M_∞<1时,0.9≤ω_b≤1.7,0.9≤ω_p<1.0;M_∞≥1时,0.8≤ω_b≤0.9,0.8≤ω_p≤0.9,这里ω_b,ω_p分别为局部亚音速点和超音速点的松弛因子。 算例为七种不同外形的细长体,计算结果与实验符合尚好。 文中对网格、初场、迭代方法、松弛因子等有关收敛性、收敛速度问题进行了探讨。在局部线化条件下,对定常小扰动轴对称势流的差分方程,进行了线超松弛迭代的稳定性和收敛性分析。数值计算经验与理论分析所得结论相符。     

垂直于物面的横截面上流态的拓扑

对于旋成体有迎角绕流，本文研究了垂直于物面的横截面上的流态结果。利用解析方法，给出了物面轮廓线上、迎风和背风对称线上以及全横截面上截面流线方程奇点的分布规律以及奇点沿流动方向的演变规律，从而可定性地确定截面流态的结构。通过求解ＮＳ方程数值模拟了钝锥有迎角高超声速绕流流场，计算获得的垂直于物面的横截面流线的拓扑结果和理论分析完全一致。     

亚超音速全机纵向气动系数数值计算方法的改善与发展

本报告是在亚超音速位流解基础上改进发展的数值计算方法。用面分布有限基本解计算作用在物面上的压力,用数值积分方法计算作用在物体上的力与力矩。本文采用联合流场、有效剖面吸力概念以及顺流向片条作可压缩紊流附面层的积分解,来改进阻力特性的数值计算。 数值方法已编制成计算程序。 本文提供了用这个计算程序计算的若干算例,并与实验结果作了对比。比较表明,计算结果与实验值有良好的一致性。     

机翼跨音速无粘绕流数值计算

本文介绍机翼跨音速定常无粘绕流的一种数值计算方法。选用精确速势方程作为问题的数学模型。在直角坐标系中,经过适当的坐标变换,先把后掠翼变成矩形翼,再把无限的物理空间变成有限的计算空间。在计算空间中,用有限差分格式(亚音速区用中心差分,超音速区用旋转差分)对精确速势方程离散化。差分方程形成的代数方程组用线松弛迭代求解。     

准静弹性飞机升力分布及气动导数计算

本文介绍一种准静弹性飞机气动特性的数值计算方法。可以用它来计算亚音速时弹性飞机的升力分布及11个主要的纵向气动导数。 文中附有算例,并对准静弹性飞机动态特性方面的一些问题进行了分析和讨论。 本文所得的导数可直接用于弹性飞机操纵性与稳定性的分析。     

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

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

利用气动弹性的柔性振翼模型建立及分析(英文)

Making use of modal characteristics of the natural vibration of flexible structure to design the oscillating wing aircraft is proposed. A series of equations concerning the oscillating wing of flexible structures are derived. The kinetic equation for aerodynamic force coupled with elastic movement is set up, and relevant formulae are derived. The unsteady aerodynamic one in that formulae is revised. The design principle, design process and range of application of such oscillating wing analytical method are elaborated. A flexible structural oscillating wing model is set up, and relevant time response analysis and frequency response analysis are conducted. The analytical results indicate that adopting the new-type driving way for the oscillating wing will not have flutter problems and will be able to produce propulsive force. Furthermore, it will consume much less power than the fixed wing for generating the same lift.     

超临界翼型射流环量控制的数值计算

为提高飞行器的升力特性,改善起飞着陆和机动性能,研究了射流环量控制对超临界翼型气动特性的影响。基于计算流体力学,对影响计算效率和精度的网格参数进行了优化设计,得到一种高效高精度网格的剖分方法;并以此为基础,建立超临界翼型环量控制的计算模型,求解雷诺平均N-S方程,分析了环量控制在不同动量系数和迎角下的作用效果,数值计算结果与实验结果基本一致。结果表明:环量控制大大提高了翼型的气动特性,且在小迎角或负迎角就可以产生较大的升力,但是动量系数并不是越大越好。计算结果可为进一步研究环量控制对飞行器的影响奠定理论基础。     

阻拦索断裂对螺旋桨舰载机着舰安全影响数值分析

喷气动力舰载机着舰拦阻滑跑，如阻拦索断裂，其逃逸复飞的概率极小。螺旋桨动力舰载机零升阻力大，推重比小，其安全复飞的能力值得研究。本文基于建立的螺旋桨舰载机逃逸复飞仿真模型（含阻拦索工作模型、发动机动力响应模型、升降舵操纵模型与气动力的动力影响修正模型），数值模拟了E-2C舰载预警机着舰阻拦索断裂情况下，其逃逸复飞的过程。仿真计算显示对象飞机在不同气动力、离舰速度与舵面操纵逻辑状态下，其纵向动力学方程中敏感参数与航迹下沉量的动态变化，结合视频数据分析其复飞成功的原因。研究表明：动力对螺旋桨舰载机俯仰力矩与升力特性的影响是其逃逸复飞成功的关键。动力影响使对象飞机的俯仰力矩曲线上移0.15，8°迎角下纵向静稳定性减小85%，升力线斜率增大29.7%、最大升力系数增大39%。这显著改善了螺旋桨飞机逃逸复飞状态下俯仰操纵的敏捷性，升降舵操纵效率与失速特性。动力影响使螺旋桨舰载机可在较小的加速度、离舰速度与有限的留空时间情况下，迅速改变其航迹角，减小航迹下沉量，保证逃逸复飞安全。     

战斗机机敏性尺度评估的一种方法

本文提出了一些机敏性度量尺度,研究了它们的分类,并试图对轴向和俯仰机敏性进行量化。用了功率攻击参数和功率损失参数量化轴向机敏性;用了上仰速率和下俯速率量化俯仰机敏性。因此,本文所引入的和所研究的内容可作为进一步研究机敏性时的参考和借鉴.     

CFD/CSD方法分析动力效应对民机气动特性影响

基于Reynolds average Navier-Stockes(RANS)的三维Navier-Stokes流场控制方程耦合结构静力学方程时域分析方法，研究了带有发动机的民用飞机其动力效应对全机气动性能的影响。首先采用数值方法对发动机进排气边界条件进行了模拟，分析了带动力的涡扇发动机模型的流场，并将计算结果与实验进行比较，验证边界条件处理的准确性；以此为基础，考虑结构弹性变形，采用计算流体动力学/计算结构动力学(CFD/CSD)耦合的方法，分别对通气和带动力的翼吊发动机全机的气动性能进行了研究。结果表明:基于通气构型预测的升阻力系数，气动载荷和压心位置与考虑动力效应后的计算结果存在明显不同。弹性变形又会加剧这一差异，使得全机的升阻比下降约12.6%，升力系数下降约8.9%,压心位置后移。数值算例显示，在靠近发动机区域气动载荷受动力效应影响显著，远离该区域，弹性变形效应占主要影响因素，因此在进行带动力效应的民机气动性能分析时，考虑弹性变形的影响是十分必要的。      

C_(nβ.dyn)作为大迎角稳定性判据的分析

 本文通过具体算便分析了C_(n β.d y x)迎角稳定性判据的可靠性,对使用这一判据可能造成失误的原因作了归纳,得出了一些便于使用这一判据的结论。     

亚超音速机翼最佳弯扭综合设计计算方法

本文介绍一种亚超音速机翼最佳弯扭综合设计的计算方法,它应用了有限基本解方法。分别在亚超音速各选取一个设计点(M数和C_L),进行机翼弯扭设计,其目的是减小与升力相关的阻力。在此基础上,顾及亚音速和超音速这两个设计点的气动力特性,还要兼顾到飞机其它性能和结构上实现的可能性,进行机翼的综合设计。本文分别给出了亚音速最佳弯扭设计,超音速最佳弯扭设计和综合设计的计算结果。经过分析表明,计算结果是合理的。     

基于总体布局参数的飞机静稳定性研究

基于飞机初步设计阶段所得到的总体布局参数,建立飞机气动导数数值模型,对亚声速静稳定性进行估算。最后对所建立的数值模型进行算例分析,并与计算流体动力学(CFD)方法所求的气动导数、实际飞行试验数据进行对比,结果表明,基于总体布局参数计算的气动导数方法快速、合理。