三角翼机翼摇滚主动控制多学科耦合数值模拟

细长机身和大后掠机翼气动构型的飞行器大攻角飞行时,由于缺少横向阻尼,易发生以绕体轴滚转振动为主的摇滚运动,飞行安全受到严重威胁。针对三角翼摇滚问题,采用动网格技术,建立了气动、运动和控制多学科耦合的数值模拟方法。通过耦合非定常Navier-Stokes方程、刚体运动方程和经典控制律,采用控制面差动偏转的方式对三角翼摇滚主动控制过程进行了数值模拟,并分析了不同控制状态下三角翼受控滚转的运动特性。在来流马赫数为0.3的条件下,实现了80°后掠三角翼摇滚现象的有效控制。     

战斗机机翼摇滚特性研究

为了研究战斗机的机翼徭滚特性，运用风洞试验和数值模拟手段，对一典型三角翼布局开展了研究工作。风洞试验研究探讨了不同攻角和初始角位移因素对机翼摇滚特性的影响；运用非定常建模技术建立了机翼摇滚过程中的滚动力矩系数的表达式并进行了机翼摇滚的数值模拟，预测了发生机翼摇滚的临界攻角和轴承阻尼系数对摇滚特性的影响。最后对机翼摇滚的发展、稳定阶段的能量转换进行了讨论。研究结果表明机翼摇滚的数值模拟与试验结果具有较好的一致性。     

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

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

F12全机动态特性数值模拟

采用非定常正弦振荡和准定常旋转坐标系两种数值模拟方法,对F12全机模型三个攻角进行动态特性数值模拟,得到飞行器的纵向组合动导数和阻尼动导数,计算结果与试验数据对比具有良好的一致性,同时用无粘和粘性两种计算方程进行准定常旋转坐标系模拟,计算结果一致,表明两种数值模拟方法准确性高,计算速度快,可为设计提供快速可信的数据支持。     

多体分离非定常气动特性数值模拟

发展建立了用于研究多体分离的非定常流场数值模拟方法,该方法采用网格弹性变形的结构网格动网格技术,采用有限体积法求解ALE(Arbitrary Lagrangian-Eulerian)描述形式的三维非定常Euler方程,与刚体运动方程耦合,数值模拟了子弹穿越母弹头部激波的子母弹分离过程,并与采用攻角修正的准定常计算结果进行了对比分析,表明多数情况下两种数值方法得到的结果比较接近。     

弹丸旋转空气动力效应非定常数值模拟

准确计算马格努斯力和力矩对旋转弹箭设计、弹道计算和稳定性研究都至关重要。采用非定常雷诺平均Navier-Stokes(RANS)方程对高速旋转的SOCBT弹丸进行了数值模拟,纵向气动特性在整个攻角范围内都较好地与试验值保持一致,而侧向气动特性马格努斯力和力矩在小攻角范围与试验数据吻合较好,在大攻角范围却存在一定的差异。采用延迟分离涡模拟(DDES)方法的计算结果有较为明显的改善,对比研究表明分离点位置对马格努斯效应有着显著影响。表明DDES方法对于提高旋转弹箭马格努斯效应的数值模拟精度有较大的潜力。     

乘波体荷兰滚模态特性研究

乘波体非轴对称、扁平、大长细比的几何外形特点决定了其存在严重的横航向耦合动稳定性问题。目前对乘波体横航向稳定性的研究还相对较少,且一般只针对单一外形,对飞行器设计未能获得有指导意义的定性、定量结论。以幂次乘波体为研究对象,首先引入设计参数kw和φ描述其外形特点,然后结合CFD数值模拟和Kriging代理模型,获得了整个设计参数空间内乘波体的静/动导数,进而通过求解飞行动力学耦合方程特征根获得乘波体的横航向耦合动稳定荷兰滚模态特性。定义了荷兰滚动稳定性导数的概念,推导了荷兰滚阻尼近似表达式,解释了不同攻角下荷兰滚模态发散/收敛的成因,获得了荷兰滚阻尼随设计参数和攻角的分布规律。推导了荷兰滚频率近似表达式,获得了荷兰滚频率随设计参数和攻角的分布规律。     

有侧滑角椭圆弹身大攻角状态绕流的数值模拟

用数值模拟方法研究了椭圆细长体在大侧滑角、大攻角状态下的流动特性。数值模拟的出发方程为三维可压缩流的全Ｎ－Ｓ方程，数值格式为ＴＶＤ格式。通过数值模拟，研究了椭圆截面细长体在Ｍ∞＝２．５、侧滑角β＝４５°状态下，攻角从１０°￣３５°范围内流场中的分离特征。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

三角翼大攻角绕流场数值模拟研究

本文基于N-S方程采用Jameson的显式中心有限体积格式和Roc的上风格式对三角翼大攻角绕流场进行初步数值模拟，空间网格为H-0型网格。计算结果表明两种计算格式均可较好的模拟大攻角绕流场。本文计算结果为今后开展全机大攻角绕流场数值模拟提供了计算经验。     

80°／65°双三角翼滚转稳定特性预测研究

在FL-23风洞中开展了80°／65°双三角翼大迎角下的滚转特性研究．包括静态测力试验．动导数试验和自由滚转试验,通过静态测力试验及动导数试验获得了双三角翼模型在大迎角条件下的滚转力矩特性以及动导数特性,从而对双三角翼大迎角条件下的滚转运动特性进行了预测．最后通过自由摇滚试验对预测结果进行了验证。研究结果表明随着模型迎角的增加,双三角翼呈现不同的滚转运动形态,包括静态稳定、舣周期震荡、准极限环摇滚．通过静态气动力及动导数可以较准确地对模型的运动形态及对应的迎角范围进行预测。     

一种近空间高超声速飞行器滚转稳定性研究

研究近空间高升阻比构型高超声速飞行器面临的复杂稳定性问题,对其机动飞行安全和操控具有重要意义。为此,通过静、动态风洞试验和数值模拟,开展了类HTV-2飞行器布局典型马赫数和迎角下,滚转方向静稳定性、单自由度和两自由度俯仰/滚转耦合运动动稳定性研究。结果表明：飞行器滚转方向为静稳定;自由振动的滚转单自由度运动为衰减振动曲线,动导数为负值,表现为动稳定;强迫俯仰/自由滚转两自由度耦合运动情况下,由于惯性力矩和非对称滚转力矩共同作用,导致滚转出现振幅周期性变化的极限环运动,甚至失稳发散。     

小迎角高超声速非旋转钝锥非线性运动分析

采用分段拟合技术对小迎角高超声速(马赫数Ma=6)非旋转钝锥双平面拍摄风洞自由飞试验结果进行了分析,获得了非旋转钝锥在小迎角高超声速下气动导数非线性的具体形式,分析了非旋转钝锥在小迎角下的非线性运动特点。研究发现在小迎角范围内钝锥的动导数系数均呈现明显的非线性,而静导数系数非线性较弱,可近似为线性。各组试验的静、动导数系数的非线性形式表明,除模型Ⅱ外的其余模型均为两方向振幅不同的极限圆锥运动,模型Ⅱ为极限平面运动。每组试验两个平面的静导数系数在小迎角范围内保持基本相等。而无论出现极限圆锥运动的四组试验还是出现极限平面运动的一组试验,模型在俯仰和偏航两个方向的动导数系数均存在较大差异。     

迎角动态变化对二元高超声速进气道气动特性的影响

针对高超声速飞行器受到扰动后迎角可能会发生瞬时大幅度改变或振荡的问题,对来流马赫数为4.03、迎角动态变化的二元高超声速进气道流场进行了非定常数值模拟,分别在进气道等速上仰和按正弦波振荡条件下研究了迎角动态变化对进气道气动特性的影响。计算结果表明:迎角动态变化时,进气道流场特征和性能参数的变化规律和稳态时情况基本保持一致,但是存在明显的迟滞现象;二元高超声速进气道动态上仰的迎角速率越大,进气道发生不起动时的迎角值越大;在进气道下壁面附近的低速区,非定常效应的影响显著,受进气道固壁面运动的影响,低速区的产生、发展及消失影响进气道气动特性。     

SDM 标模大攻角动导数试验

介绍了ＣＡＲＤＣ高速所１．２ｍ风洞大攻角的动导数试验装置与测量系统以及在１．２ｍ风洞中对标准动态模型（ＳＤＭ标模）所作的一系列风洞试验结果。试验结果与国外文献数据具有较好的一致性。试验结果随马赫数和攻角的变化表现出明显的非线性,而减缩频率的影响并不显著。天平与测试系统的重复性精度较高。     

用极大似然法辨识轴对称飞行器的气动参数

采用极大似然法完成了轴对称飞行器的非线性气动参数辨识。根据飞行弹道的特性,采用考虑了横侧向运动耦合的纵向三自由度飞行动力学模型;为了反映气动参数随马赫数的变化,给出了一种考虑马赫数变化的气动参数数学模型,并采用了分段辨识的技巧;通过采用建模前估计法的辨识结果作为极大似然法的初值,保证了极大似然法的收敛,计算结果表明：这种算法对大攻角下的非线性气动参数能获得满意的辨识结果。     

方形截面飞行器上仰机动对滚转特性影响的数值模拟

飞行器从中小迎角至大迎角范围内,由于背风区流动分离形态的演化,静态气动特性特别是横侧向气动特性也随迎角显著变化,可能诱发复杂的滚转运动。但飞行器一般是上仰机动时,才从平飞状态快速拉起至大迎角,此机动过程对横侧向气动特性和滚转运动可能产生较大影响。本文发展了刚体动力学方程和Navier-Stokes方程的松耦合求解技术,并通过数值模拟航天飞机脱落碎片的六自由度运动轨迹进行了验证。针对背风区涡流形态及横侧向气动特性复杂的方形截面飞行器,数值模拟研究了其不同迎角下的静态滚转气动特性、自由滚转运动特性,以及上仰机动时不同拉起速率对滚转运动特性的影响。结果表明,对于此飞行器,静态时存在临界迎角约为13°,当迎角小于临界迎角时,滚转方向是静不稳定的,诱发快速滚转运动;当迎角大于临界迎角时,滚转方向是静稳定的,其滚转运动是收敛的。但上仰机动时,滚转运动的形态还与拉起速率相关,即使拉起的终止迎角大于临界迎角,如果拉起速率较慢,也可能出现快速滚转运动。     

Aerodynamic investigation of twist angle variation based on wing smarting for a flying wing

In this paper, the effects of twist angle variation on aerodynamic coefficients and flow field on the wing with wing smarting approach are studied using numerical simulation. The simulation was performed using incompressible Reynolds-Averaged Navier-Stokes (RANS) equations based on the two-equation k-ω Shear Stress Transport (SST) turbulent model for flow speed 30 m/s and a Reynolds number of 69000. Investigations have been carried out for several twist angles and at a specific range of angles of attack. The twist applied is the type of geometric twist (wash-out), which is linearly distributed along the span. The test case is a lambda-shaped tailless aircraft with a wing fracture on the trailing edge, and a sweep angle 56°. The results show that with increasing twist angle, the aerodynamic efficiency improves over a wide range of angles of attack, but at 0° angle of attack it will decrease significantly. By increasing the angle of attack, the effect of twist on the flow field and aerodynamic coefficients will gradually decrease; hence, at a certain amount of angle of attack, the effect of twist will stop, that angle is called the neutral brink angle. Longitudinal stability analysis shows that by growing the twist angle, the conditions required for longitudinal stability are satisfied, and the pitch-up phenomenon will be delayed.     

微型飞行器低雷诺数矩形扑翼非定常气动特性的数值模拟

采用双时间步方法求解三维可压缩非定常N-S方程,数值模拟了微型飞行器低雷诺数矩形扑翼的非定常绕流,首先将得到的结果与文献进行了对比,数据间具有较好的一致性.然后针对不同的展弦比、减缩频率及初始攻角,计算了矩形扑翼的非定常气动特性及表面流态和动态压力分布,并分析了翼尖涡对扑翼非定常气动特性的影响.