基于耦合隐式方法的复杂外形绕流计算研究

飞机复杂外形的气动特性数值模拟在飞机型号研制中有重要的作用,也是数值模拟中比较复杂的问题。笔者采用基于Octree技术的方法生成非结构四面体网格,用耦合隐式离散方法进行有限体积离散,利用此方法对某型飞机带外挂和舵面偏转情况的气动特性进行了数值模拟,将计算结果与实验结果作了比较,收到较好效果。     

机翼静气动弹性数值模拟研究

通过求解基于三维非结构网格的欧拉方程计算机翼气动特性,采用结构柔度影响系数方法进行机翼静弹性变形的计算,利用三维非结构弹性网格技术进行几何外形变形下的网格修正,将以上三部分进行耦合迭代.实现了机翼的静弹性特性的数值模拟.采用开发的数值分析程序,进行了某支线飞机机翼静弹性特性的计算分析,计算结果合理,表明所发展的数值模拟方法可用于机翼静弹性气动特性的分析研究.     

某栅格翼导弹减阻外形设计及其阻力特性数值研究

介绍和分析了国内外栅格翼的研究现状和技术难点。为了详细了解栅格翼气动布局栅格内部的流动特性,以及其与阻力设计满足飞行器性能要求的栅格翼,本文开展了几种不同单栅格形状外形的气动特性数值计算。针对传统数值计算中遇到的栅格翼计算网格生成困难、工作量大的问题,利用非结构空间直角网格对于复杂外形易于生成的特点对单个栅格的外形进行了初步选型研究。对于常规导弹和栅格翼导弹的气动性能进行了对比。     

结构/非结构混合网格数值模拟栅格翼

为了详细了解栅格翼内部的流动特性,设计满足飞行器性能要求的栅格翼,本文开展了多种不同单栅格外形的气动特性数值计算,对单个栅格的外形进行了初步选型研究;然后采用结构网格、结构与非结构混合网格求解NS方程,对带有两片栅格翼的导弹进行了数值模拟,结构网格计算结果与结构非结构混合网格计算结果一致,并与国外文献的计算结果和风洞试验结果进行了对比,三者的结果吻合较好.采用混合网格比结构网格的网格规模小许多,可以节省大量计算机内存和计算时间,是一种比较好的栅格翼数值模拟方法.     

应用混合网格分析复杂积冰翼型气动性能

通过发展一种点-点对接的结构/非结构混合网格生成方法,避免了复杂积冰外形难以划分结构化网格的问题,并对复杂积冰翼型的气动性能进行了分析计算.在靠近积冰边界的内层采用了非结构网格以拟合复杂的边界,在非结构网格外采用结构化网格以节省存储空间和计算时间.为了检验网格对选用湍流模型的影响,整个流动区域分别采用SST和SA湍流模型,求解了雷诺平均N-S方程.数值计算结果表明SST模型更适于模拟复杂分离流动,积冰对翼型的气动性能造成了严重的影响.     

基于非结构网格求解器的CHN-T1标模气动特性计算研究

高可信度计算流体力学(CFD)技术在大型运输机精细化气动设计中发挥着越来越重要的作用。虽然目前基于RANS方程的数值模拟技术已经应用于大型运输机复杂外形的气动性能计算与分析中,但其可信度仍需进一步确认。针对AeCW-1组委会发布的CHN-T1运输机标模,采用基于RANS方程的非结构网格求解器TNS进行了气动性能计算,并与组委会提供的风洞试验结果进行了对比。此外,还研究了网格量和湍流模型等因素对巡航气动特性计算结果的影响,以及支撑干扰和静气弹变形对跨声速抖振特性的影响。结果表明,TNS预测的巡航升阻比及跨声速抖振特性等气动性能与试验值吻合良好,计算结果完全位于试验值误差带内,验证了所采用的非结构网格求解器在大型运输机复杂外形气动力预测中的可靠性。     

带反馈控制的基本带翼外形的耦合气动/运动数值模拟研究

介绍了一种先进的、时间精度的多学科数值模拟平台的开发和应用,该平台可用于预测飞行器机动过程中的非定常气动特性和飞行力学特性.通过耦合求解非定常Navier-Stokes方程、刚体动力学方程及控制律,对基本带翼外形进行了虚拟飞行数值模拟研究.求解器采用基于结构网格的有限体积法,时间方法采用双时间步方法.使用反馈控制器获得舵偏控制指令,重叠网格方法用来模拟物体间的相对运动.模拟结果显示控制器的增益必须正确选择才能获得正确的响应.仿真结果表明本文发展的虚拟飞行数值模拟方法能够预测包含动边界的复杂非定常流动,具有较强的工程适用性.     

带控制律导弹投放数值模拟

基于非结构重叠网格技术,耦合舵面控制律模块,采用非定常雷诺平均NS方程,建立了模拟控制律作用下导弹投放分离动态特性的CFD数值方法.非结构重叠网格生成技术以“物面相交”方法确定洞边界,在多套网格重叠部分采用算术平均进行物理量的计算传递.采用机翼/挂架/带舵外挂物模型,以及不同的舵面控制规律展开了对比计算,分析了结果.并比较了某导弹有无控制律作用时从某战斗机分离的动态特性,研究了增加控制律的必要性.     

复杂积冰翼型气动性能分析

为了克服复杂外形难以划分结构化网格的问题, 采用结构/非结构混合多块网格对复杂积冰翼型的气动性能-进行了分析计算.在靠近积冰边界的内层采用了非结构网格, 在非结构网格外采用结构化网格以节省存储空间和计算时间.不同类型网格之间采用了点点对接的方式以简化通量重构的处理方式.整个流动区域采用SST湍流模型, 求解了雷诺平均N-S方程.数值计算结果表明结构/非结构混合多块网格在复杂积冰翼型气动性能计算上的可行性, 积冰对翼型的气动性能造成了严重的影响.      

外挂物干扰流场特性数值仿真研究

飞行器研制中对于带外挂物的复杂外形气动特性分析研究一直是个难点,对外挂物与挂架不同缝隙的复杂机翼-挂架-外挂物组合外形采用结构重叠网格数值方法,快捷高质地完成了网格生成,通过求解带有k-e湍流模型的Navier-Stokes方程组,得到多个计算工况下的绕流外流场,分析了外挂物在不同工况下的气动特性,结果表明,CFD数值方法可方便快捷地对复杂干扰流场进行数值计算,可为飞行器相关设计提供技术支持。     

偏转弹头导弹动力学建模

偏转弹头控制是一种新型的导弹控制方式,精确动力学模型的建立对于导弹性能、稳定性和控制的分析是必不可少的。本文采用多体动力学建模方法中的凯恩(Kane)方法对偏转弹头导弹建立了动力学模型。首先建立合适的坐标系,选取8个变量的速率为广义速率,之后推导了弹身弹头的运动学方程,接着得到广义主动力和广义惯性力,其中将控制弹头偏转的作动力矩视为主动力矩而非理想约束考虑进广义主动力的计算中,最终得到导弹的8个标量动力学方程。推导过程显示凯恩方法比牛顿-欧拉方法推导过程更为简化。利用得到的精确动力学模型进行动态仿真,通过比较研究了采用弹头质心位于弹头偏转铰链中心假设下的简化模型对导弹动态响应的影响,结果显示,在弹头偏转较慢时简化模型会带来较大的动态响应误差。     

用隐式多重网格法计算实用外形无粘绕流

网格展弦比的好坏严重地影响多重网格的效率。文中细致地处理了多重网格算法的细节,了网格展弦比的影响,发展了一种计算实用外形绕流的隐式多重网格解法,计算效率提高了５倍。     

翼-身-尾组合体的网格生成及 Euler 方程数值模拟

采用代数方法快速生成翼-身-尾组合体三维O-C型贴体网格。该代数网格生成方法使用双边界法,可保证物面附近近似正交,真实地模拟了机身头部及与机翼、尾翼交界处的形状,提高了代数方法生成网格的质量。并用中心格式有限体积法求解跨声速Euler方程,以某翼-身-尾组合体模型为例,计算结果与实验符合良好。     

亚跨声速风洞机弹干扰试验技术

介绍了一种利用小尺寸风洞进行亚跨声速挂弹/载机干扰试验的试验技术,该试验技术是在小尺寸风洞内进行挂弹/载机局部干扰试验,可以测量挂弹投放过程中作用于挂弹及其折叠尾翼上的气动力和舵效率。为验证该试验技术,选取了某空射型导弹进行了挂机干扰试验验证,试验前对局部干扰流场和完整干扰流场进行了数值计算,计算结果表明局部干扰流场与完整干扰流场流态基本一致。试验结果表明,小风洞对局部流场的模拟可以取得合理的试验结果,该试验技术具有可行性。     

基于预估校正和嵌套网格的虚拟飞行数值模拟

 针对导弹虚拟飞行数值模拟问题,发展了空气动力学/飞行力学数值计算方法和软件。控制方程为非定常雷诺时均Navier-Stoker(RANS)方程和刚体六自由度运动方程;流场求解器为有限体积法结构网格求解器,时间推进采用双时间步法,湍流模型为Spalart-Allmaras一方程模型;采用Adams预估校正法实现飞行力学方程与流场控制方程的耦合计算;使用嵌套网格方法模拟多体运动。首先模拟了美国国家航空航天局(NASA)窄条翼导弹模型纵向虚拟飞行,研究耦合方式和时间步长的影响。仿真结果表明,双时间步三阶Adams耦合方法,同等精度下可以显著增大时间步长,缩短仿真时间。最后,采用该方法模拟了导弹自由摇滚特性和纵向虚拟飞行,模拟结果与试验值吻合较好。     

一种高效的组合动导数分量数值模拟算法

基于传统网格速度法,发展了在流场网格不变的情况下,数值模拟沉浮运动、纯俯仰运动和俯仰振荡运动的算法,提出了一套直接求解旋转阻尼导数和洗流时差导数的方法,该方法不需要实时更新网格,减少了计算时间和所需内存,避免了负体积的出现。首先通过NACA 0006的阵风响应和NACA 0012翼型的强迫俯仰振荡,验证了程序和方法的正确性;然后数值模拟了BFM(basic finner missile)标准模型的沉浮运动、纯俯仰运动和俯仰振荡运动,并计算了其旋转阻尼导数、洗流时差导数和组合动导数,所得结果与动网格方法和风洞试验值吻合,最大误差不超过4%,进一步验证了该方法的正确性。     

复杂外形飞行器多块结构贴体网格生成

通过求解椭圆型方程生成三维贴体坐标网格,采用分块技术生成多块结构网格体系。每个子区域转换成计算空间内的立体域,具有自身的曲线坐标系,从而所形成的整个网格系统不会过分扭曲和不规则,这将有利于对于复杂边界的真实飞行器绕流流场的求解。文中讨论了控制函数对网格生成的影响。     

Nonlinear Estimation Via Kalman Filtering

The extended Kalman-Bucy filtering algorithm is applied to the problem of estimating the system parameters of a single-channel missile attitude control system. In particular, noisy observations of body angular rate, normal acceleration, and control surface deflec tion are used in the algorithm to obtain estimates of body angular rate, angle of attack, control surface deflection, air density, and missile velocity. Computational results obtained for this problem are presented. These results are evaluated to determine the applicapility of this estimator as a part of an adaptive control scheme. In addition, a general approach to the formulation of practical nonlinear estimation problems is suggested.     

Changing law of launching pitching angular velocity of rotating missile

In order to provide accurate launching pitching angular velocity（LPAV） for the exterior trajectory optimization design, multi-flexible body dynamics（MFBD） technology is presented to study the changing law of LPAV of the rotating missile based on spiral guideway. An MFBD virtual prototype model of the rotating missile launching system is built using multi-body dynamics modeling technology based on the built flexible body models of key components and the special force model.The built model is verified with the frequency spectrum analysis. With the flexible body contact theory and nonlinear theory of MFBD technology, the research is conducted on the influence of a series of factors on LPAV, such as launching angle change, clearance between launching canister and missile,thrust change, thrust eccentricity and mass eccentricity, etc. Through this research, some useful values of the key design parameters which are difficult to be measured in physical tests are obtained. Finally,a simplified mathematical model of the changing law of LPAV is presented through fitting virtual test results using the linear regression method and verified by physical flight tests. The research results have important significance for the exterior trajectory optimization design.     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。