尖前缘任意翼型的超音速和高超音速俯仰安定性导数

本文推广了文[1]的欧拉方程摄动法,讨论了尖前缘任意翼型俯仰安定性导数的计算方法,并应用该方法计算了双圆弧翼型的俯仰安定性导数,其数值结果与文[2]的超音速线化位流方程有限差分法的数值结果比较,更接近于实验值。在激波附体的条件下,本文的方法可用于计算大迎角尖前缘任意翼型的超音速和高超音速俯仰安定性导数。     

某改型飞机时差动导数对尾旋时间历程估算的影响

提供了利用欧拉方法计算的某改型飞机尾旋时间历程结果,文中比较了使用和没使用洗流时差动导数数据对计算结果的影响,说明洗流时差动导数数据在这种模拟计算中是非常重要的,它使计算结果更合理。     

针对旋翼流场模拟的自适应网格算法

针对基于格心格式求解器的旋翼流场模拟,提出了相应的自适应笛卡尔网格的数据存储结构及自适应算法。给出了相应的单元处理策略,简化了对自适应笛卡尔网格的处理;对于频繁的自适应加密过程中产生的大量重复点,采用了高效的交替数字树算法(Alternating digital tree,ADT)予以删除;对于自适应疏化过程中产生的大量无用点,提出了标记-删除-移动(Mark,delete,move,MDM)算法予以快速地删除,减少了不必要的计算资源消耗。对CaradonnaTung旋翼在不同悬停状态下进行了模拟验证,对比了压力分布系数与桨尖涡位置。之后对HELISHAPE 7A旋翼在前飞情况下进行了模拟验证,计算值与实验值吻合。此外,求解器对桨尖涡的捕捉效果得到了明显的提高,表明本文方法具有良好的有效性与鲁棒性。     

任意迎角二维平板的超声速和高超声速俯仰振动导数

应用频率的摄动法,本文得到了任意迎角二维平板机翼的超声速和高超声速俯仰振动导数。该导数是由包括直到缩减频率的五阶项所组成的。数值结果表明高阶项的贡献在较低马赫数或大迎角情况下是很重要的。作为平均迎角趋于零的特殊情况,本文的数值结果与超声速线化的位流理论是一致的,在高超声速情况与活塞理论一致。     

有限元计算模型的试验修正

本文阐述了有限元计算模型试验修正的一种方法,推导了求特征值和特征向量对结构参数导数的二阶矩阵摄动近似解,并在直升机桨叶的有限元计算模型的试验修正中得到成功应用。     

旋转流场下的振荡动导数试验技术研究

为研究飞机在旋转流场下的非定常气动特性,中国空气动力研究与发展中心低速所在Φ5m 立式风洞开展了旋转流场下的振荡动导数试验技术研究。本文推导了在旋转流场下识别组合动导数的方法,介绍了试验设备,获得了在旋转的同时,由振荡产生的3个组合动导数,并对试验结果进行了分析与讨论。将单自由度动导数结果与Φ3.2m 风洞试验结果进行了对比,旋转/振荡耦合试验结果表明:旋转运动使得俯仰组合动导数变得不稳定,而对于横向组合动导数,大转速则会显著增大非线性。该试验技术能够为研究旋转流场下的非定常气动特性提供一个有效的试验平台。     

结构固有振动特征值问题矩阵摄动法的一个拓广

本文探讨振动结构的物理、几何参数与固有模态之间的关系。处理这类问题的常见方法是特征值问题摄动法。但该方法得到的解只有当参数作小扰动时才有效。为克服这一缺陷,本文将前人已得到的模态关于参数的导数表达式看成是一组常微分方程式,从而运用常微分方程的定性研究方法和数值计算方法来研究这类问题。文中用这一拓广的方法探讨了模态在参数空间上的连续性及重频点对模态演变的影响,并对一个轮盘叶片模型进行了数值计算和实验分析。     

重根特征向量导数计算的特殊迭代法

在重根特征向量导数计算方法的发展中，为了从无限个特解中确定出唯一的通解，提出过不同定解条件，但正确定解条件是唯一的。员后公认的定解条件为Z^TMZ Z^TMZ=-Z^TMZ,这里认为Z^TMZ≠Z^TMZ，由此便引出物理短阵(刚度阵和质量阵)的二阶导数。为了避免物理阵二阶导数的计算，本文在满足前述定解条件的前提下，利用一种特殊的迭代格式回避了物理阵二阶导数的引入。同时该迭代过程可直接获得通解，不同于目前流行的做法：先由一个不定支配方程求其特解，然后才由定解条件确定通解。另外，数值表明，动柔度迭代式的精度可与直接迭代式相匹敌，可是动柔度迭代式用于许多特征向量导数计算时，却省时得多。     

基于非定常气动力建模的动导数仿真实验

应用系统辨识理论,利用纵向大振幅非定常气动力实验数据,在频率域内建立了基于Ｆｏｕｒｉｅｒ 变换法的非定常气动力数学模型,并发展了一种获取常规动导数的仿真实验方法,分析了迟滞环的物理意义。对６０°三角翼和Ｆ－１８ 飞机模型进行动导数仿真计算,研究了不同振幅、频率及迎角变化对动导数的影响。结果表明,Ｆｏｕｒｉｅｒ 变换模型可以用于动导数仿真计算,使大振幅非定常实验与小振幅常规动导数实验合二为一。迟滞环的方向反映了飞机运动的阻尼特性。     

短钝外形飞行器自由振动动导数试验技术

为研究短钝外形飞行器的动稳定特性,基于自由振动动导数试验方法在1.2 m量级亚跨超声速风洞中建立了动导数测量试验技术.通过新设计的弹性铰链和轴承铰链解决了短钝外形飞行器弹性支撑和低频振动模拟问题.利用新建立的试验装置研究了马赫数、迎角、减缩频率对动稳定特性的影响.在短钝外形飞行器气动力特点下,新设计的弹性铰链能够满足模...     

超声速化学反应流场迎风格式数值模拟及并行计算

基于有限体积迎风格式对超声速燃烧流场进行了的数值模拟.由于超声速燃烧流场绕流的复杂性,要求对多组分Euler/N-S方程求解的数值模拟方法应具有较高的计算精度及效率.本文引用辅助点方法建立了具有空间二阶精度的van Leer迎风矢通量分裂格式,并应用于超声速燃烧流场绕流的数值模拟.化学反应为氢气/空气十反应模型,采用考虑了化学反应特征时间的当地时间步长显式Runge-Kutta时间推进格式.对钝头体模型爆轰现象、后向台阶氢气喷射及二维内外流超声速燃烧流场模型进行了区域分裂技术的并行计算.计算结果与参考文献作了对比,得到了满意的结果.     

圆转方非轴对称喷管无粘流动的多层网格法求解

本文利用多层网格法求解三维无粘亚音和跨音速圆转方非轴对称喷管的流动,基本格式是把Ni[1]的二维有限体积积分格式推广到三维流场。边界条件采用文[2]中的“预测一修正”处理方法。为了进一步提高格式稳定性,在计算流动变量的一阶时间变化量时,采用单元体积平均法而不是文[1]中的算术平均法。引入当地时间步及多层网格方法,以加快格式的收敛速度。文中还对多层网格构造形式及其收敛效果作了讨论。通过对两种喷管模型流动的计算和与试验数据的对比,验证了该算法的精度。     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

超声速和高超声速粘性绕流QPNS方程的GLS有限元解

从完全非定常Ｎ－Ｓ方程（ＦＮＳ）中略去流向粘性导数导出拟抛物Ｎ－Ｓ方程（ＱＰＮＳ），ＱＰＮＳ技术结合经典抛物推进和拟时松弛，拟时松弛能有效地抑制解的发散。本文导出的熵变量形式ＱＰＮＳ方程具有对称性和自动满足热力第二定律，这将提高解的稳定性，伽辽金／最小二乘法（ＧＬＳ）被用来构造ＱＰＮＳ的弱解式。对于钝前缘物体，ＱＰＮＳ在前缘区并不适用，本文代之以用ＦＮＳ求解该区，由离散解得到的非对称线性方程组，对     

Euler方程无网格算法及布点技术

研究了Bitina的显式无网格算法。在该算法的基础上,给出了Euler方程无网格离散形式,运用RungeKutta显式时间推进格式推进求解。文中的耗散项采用了与非结构网格上类似的守恒型耗散算子,边界条件的处理借鉴了结构化网格处理技术。此外,本文还描述了一种区域离散布点方法,研究了点云生成的选点准则,并成功地数值模拟了二维翼型的典型绕流。     

基于SNGR方法的二维空腔噪声数值模拟

基于SNGR方法,采用有限差分法求解带源项的线化欧拉方程。数值方法采用色散关系保持(DRP)格式。低耗散色散龙格-库塔格式显式时间推进,采用无反射远场边界条件。对二维空腔,长深比为9∶1,在亚声速、跨声速和超声速多个马赫数下进行噪声预测,并与实验结果对比,两者吻合较好。     

强迫振动法提取桥梁气动导数研究

利用节段模型的自由振动法来提取气动导数，这种方法其装置简单、易于实现，但在试验中结果并不完美。原因之一是模型可能会受到涡的影响，并且在高风速时，信号受干扰严重，信噪比较低。另外激振方式的不同也会造成重复性差。强迫振动法一直被认为是提取气动导数最精确的方法，但由于装置较为复杂，研究和实施此法者目前较少。笔者给出了针对两种桥梁断面进行的一次强迫振动法实验检验，表明效果好。     

旋转导弹风洞动态测力试验技术研究

常规风洞静态气动力测量技术无法得到旋转导弹的非定常气动特性数据,需要研究在风洞中模拟旋转导弹运动特征以及对气动力实现动态测量的试验技术。在1.2 m量级超声速风洞中,研究了大长细比导弹模型旋转运动主动控制技术以及与旋转运动对应的动态测量试验技术。采用旋转导弹模型（长细比为20）对建立的试验技术进行了风洞试验验证。结果表明:采用微型驱动系统并对旋转组件与导弹模型进行一体化设计,可以对大长细比导弹模型转速进行稳定控制;建立的风洞动态测力试验技术可以对导弹模型旋转运动下的动态数据进行测量,试验数据重复性精度良好。