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

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

二元非线性机翼随机动力学行为

 为考查在高斯白噪声作用下二元机翼随机颤振动力学行为,采用数值仿真的方法对其进行了研究。首先把机翼简化成一个具有扭转和上下自由度的平板,其中扭转弹簧具有3次非线性刚度,同时假设空气是线性不可压的而且气动力是定常的,再对任意运动和高斯白噪声作用下的机翼进行建模。其次采用了蒙特卡罗仿真来求解随机微分方程的数值解,根据数值解结果进行统计分析,计算出最大Lyapunov指数。最后得出：在随机激励作用下机翼的动力学行为与确定性颤振相比有很大差异,同时随机颤振点提前于确定性颤振点。     

大攻角机翼的气动弹性计算方法

 本文利用大攻角带分离流机翼的完全非定常非线性气动力计算方法,通过与机翼运动方程的同时求解,在时间域内实现了大攻角机翼非线性气动弹性的数值模拟,根据弹性机翼各种状态下的运动过程,可以得到大攻角机翼的颤振速度等重要参数以及亚临界、超临界等飞行状态的运动规律。算例结果表明,大攻角下机翼的气动弹性问题需引起足够的重视。     

一种无人机变后掠翼机构的设计与仿真研究

提出了一种变后掠翼机构的设计方法,能根据机翼平面形状参数,合理布置变掠机构、确定机构部分主要参数,并能得到机翼的运动规律。首先提出了变后掠翼机构方案,将运动机构抽象为数学模型;然后对各参数进行详尽的设计和计算,并对机翼的运动进行了数值分析,在此基础上建立了一种变后掠翼机构设计方法。运用方法对具体的算例,建立机构三维模型,进行运动模拟仿真,CATIA运动仿真结果与MATLAB数值计算结果一致,证明了方法的有效性和合理性。     

基于 Kriging模型机翼平面外形气动优化设计

使用基于Kriging模型的优化设计方法,进行了非常规布局机翼的平面外形多目标优化设计。利用CFD技术进行机翼升力系数和阻力系数的气动计算,通过拉丁超立方试验设计生成样本点,建立了Kriging代理模型,结合多目标遗传算法对机翼平面外形进行多点多目标优化设计,最终得到了Pareto最优解集。根据设计需求,从Pa-ret0前沿选取一个非劣解作为优化结果。结果表明：陆ging模型与cFD计算误差很小,可信度高;在不问设计状态下,机翼气动性能都得到了提高,表明优化设计方法具有可行性和高效性。     

后掠机翼跨音速绕流的二级近似方法

本文采用TSDH方程计算三维后掠机翼的跨音速绕流,考虑了适用于机翼钝前缘的前缘边界条件和前缘速势方程。采用Jameson格式在不等步长格网中的推广形式,把TSDH方程离散化为差分方程组。然后,在整个计算空间内布置稀网,在机翼附近再布置密网,进行稀密网的交替迭代,以加速收敛和提高计算精度。对ONERA M6机翼的超临界无激波和有激波情况的计算表明,TSDH解与FVP解和风洞试验符合良好。     

绕三角翼纵向俯仰大迎角气动特性计算研究

 采用数值计算方法, 对三角翼从0°上仰至90°的动态流场结构进行了计算, 在此基础上, 对三角翼在上仰过程中受到横侧小扰动情况下的流场结构和气动力特性进行了计算研究。给出了三角翼纵向动态情况下的气动力系数变化, 特别是大迎角横侧力矩系数的变化特征, 并对受到横侧小扰动后横侧运动的稳定性进行了计算与分析。结果表明, 机翼的上仰运动延迟了机翼上翼面旋涡的破裂。同时, 随着机翼俯仰角速度的提高, 机翼抵抗旋涡非对称破裂的能力明显增强, 机翼运动的稳定性也明显提高。     

完全非定常气动力计算的新途径

 本文提出了一种亚音速、完全非定常机翼气动力的数值计算方法（暂称为“势差法”）,该方法对于机翼所作的任意运动,能在时间域内计算出气动力和空间速度场的变化规律。算例结果表明势差法计算结果与其它方法计算结果和实验值是一致的。     

一种改进型松耦合方法在机翼摇滚计算中的应用

提出了一种简单实用的松耦合机翼摇滚计算方法,该方法采用隐式和显式公式相结合来计算刚体动力学方程并升级网格,使计算结果的收敛性和稳定性随时间步长的增大更好,可以在较大的时间步长下得到合理的计算结果。该方法是在隐式和显式的松耦合方法基础上得到的,通过分析机翼摇滚计算中松耦合方法时间精度低的原因,发现如果基于隐式方法计算角速度和角度,由于计算力矩时网格位置滞后,导致气动力相对于运动滞后,时间步长较大时计算的机翼摇滚振幅偏大;而如果基于显式方法计算角速度和角度,计算力矩时网格位置超前,导致计算的机翼摇滚振幅偏小。通过计算80°后掠三角翼的机翼摇滚,证明了改进型松耦合方法的有效性。     

基于伴随算子的翼尖小翼优化设计

采用基于伴随算子的优化设计方法,对跨音速翼尖小翼进行了优化设计。其中,采用基于多块结构网格的流场解算器计算气动力。解算器的空间离散采用Osher格式,湍流模型采用k-ω两方程模型。优化设计结果表明,采用优化设计后的翼尖小翼,能改善机翼的气动特性,在巡航状态下翼尖小翼能增加机翼的升阻比5.8%,从数值计算的角度表明了方法的有效性。     

折叠机翼变体飞机纵向操纵性与稳定性研究

针对一种飞翼布局折叠机翼变体飞机方案,建立了相应的研究模型,使用涡格估算方法计算出机翼折叠角度对全机纵向静稳定性的影响.利用飞机纵向小扰动运动方程,分别得到了机翼展开和折叠状态的长、短周期模态,并对其动稳定性进行了讨论.结合工程估算方法和风洞试验方法,计算出内段机翼折叠过程全机力矩系数、升降舵操纵导数以及升降舵配平偏角...     

飞行器刚体运动与气动弹性相互作用问题研究

通过联立求解空气动力学基本方程、飞行动力学运动方程和弹性结构振动方程,在时间域内模拟和分析了大展弦比飞机纵向动力学稳定性问题。结合动网格技术,气动力计算采用基于欧拉方程的计算流体力学方法,结构变形和飞行姿态位置变化统一为模态表示方法,通过松耦合将飞行器姿态稳定性和结构变形稳定性施行了模拟。以某大展弦比机翼飞机为算例,研究了其刚体运动和机翼的弹性振动的相互影响。结果表明：对于具有大展弦比机翼的飞机,其机翼的低阶弹性模态易与飞机飞行中本身的刚体模态发生耦合,从而导致飞机机翼的气动弹性发散以及飞机本身刚体运动稳定性的改变。对于这类飞机,在其气动弹性和飞行稳定性的分析和设计中必须充分考虑到两种运动的相互影响。     

机翼亚音速非定常气动力计算

分析机翼非定常气动力特性既可用实验方法[1～3],也可用数值计算方法[4]。现采用非定常涡格时间历程法对机翼沉浮、俯仰、偏航和滚转振荡时的气动力进行计算。1　基本原理与计算公式设图1中t=0时机翼作定常运动;t=Δt时,机翼攻角突变1°,则机翼上产生一涡环ΔΓ1,并以V0向下游运动;t=2Δt时,机翼上又产生一涡环ΔΓ2;随着时间的推进,不断产生以V0向下游运动的涡环,从而模拟机翼的非定常绕流。当涡环以V0运动时(图2),其速度势为φP0=V0Γ4πlimy→0y∫b/2-b/2dz∫x1+V0tx11rdx上式对x、y、z分别求导得到非定常涡环的诱导速度Vx=V0Γb4πx…     

应用前缘填角消除气流分离的设计方法

提出利用前缘填角消除机翼一机身接合处气流分离的设计方法，设计方法应用了蒙皮摩擦变化与表面曲率变化相关的设计规则。求出在机翼前缘之前板的当前蒙皮摩擦分布之后，计算可以消除气流分离的目标蒙皮摩擦分布。应用设计规则，把当前的蒙皮摩擦分布转化为目标蒙皮摩擦分布，延伸机翼前缘形成前缘填角，利用流动解算器分析机翼填角，过程迭代直到收敛。     

超声速机翼－尾翼对突风动态响应的数值方法

本文以Ｇｒｅｅｎ函数为基础的边界元法，结合时间历程法，计算机翼－尾翼组合体对突风动态响应的气动性能。本文介绍的方法，特别适于研究、分析各类飞行器非定常动态响应，为飞行器的设计及其非定常气动特性分析提供了一个行之有效的计算方法。编制的计算程序快捷、精度好。本文提供了单个机翼对突风动态响应的几个算例，有的与Ｌｏｍａｘ解进行了比较，吻合得很好，说明本文所使用的方法是可行的，程序是成功的。对于机翼－尾翼组合体的数值结果，则进一步提供了尾翼受机翼洗流影响的过程。     

二维机翼非定常运动的涡流场显示——沉浮运动

介绍了在水洞中进行的二维机翼沉浮运动的流动显示实验。实验结果揭示了机翼正弦运动时尾迹涡街的特性以及与Strouhal数的关系；发现了在不同参数下存在三种前缘涡的形式；结合数值计算结果说明了二维机翼沉浮运动中产生非定常推力的机理。     

亚音速复杂组合体表面压力分布及气动力特性数值计算方法

本文从亚音速小扰动位势方程出发,给出了计算任意复杂组合体表面压力分布和纵向气动力特性的基本解的数值计算方法。把全机表面划分成许多网格,在每个网格上布置基本解以模拟物体,用物面上没有法向气流穿透的边界条件决定基本解的未知强度,从而计算出物面压力系数,积分出总气动力特性。该方法适用于任意形状的机身及任意平面形状、任意厚度、弯度和扭转分布的机翼(尾翼),机翼(尾翼)可以带有两段后缘操纵面。为降低对电子计算机存贮量、速度和稳定时间的要求,本方法采用了机翼、机身、尾翼部件之间的循环迭代求解。     

超声速机翼—尾翼对突风动态响应的数值方法

本文以Ｇｒｅｅｎ函数为基础的边界元法，结合时间历程法，计算机翼－尾翼组合体对突风动态响应的气动性能。本文介绍的方法，特别适于研究、分析各类飞行器非定常动态响应，为飞行器的设计及其非定常气动特性分析提供了一个行之有效的计算方法。编制的计算程序快捷、精度好。本文提供了单个机翼对突风动态响应的几个算例，有的与Ｌｏｍａｘ解进行了比较，吻合得很好，说明本文所使用的方法是可行的，程序是成功的。对于机翼－尾翼组     

大攻角机翼非线性颤振分析方法

 本文提出了一种计算方法分析大攻角、带分离涡机翼的非线性颤振特性。在机翼变形与广义气动力之间建立一种描述函数关系,利用传统v-g法对机翼运动方程进行求解。此外,还发展了一种迭代运算过程,使得输入的振幅幅值比例调整到与求解颤振方程所得到的比例相协调。计算结果表明本文方法与实验是一致的。     

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

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