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

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

不完全概率信息的压杆稳定可靠性鲁棒设计

将可靠性设计理论和鲁棒设计方法相结合, 讨论了不完全概率信息的压杆稳定可靠性鲁棒设计问题, 提出了稳定可靠性鲁棒设计的计算方法.把可靠性灵敏度溶入可靠性优化设计模型之中, 将稳定可靠性鲁棒设计归结为满足可靠性要求的多目标优化问题.在基本随机参数的前四阶矩已知的情况下, 通过计算机程序可以实现不完全概率信息的压杆稳定可靠性鲁棒设计, 迅速准确地得到不完全概率信息的压杆稳定可靠性鲁棒设计信息.      

螺旋桨滑流对带后缘襟翼机翼气动特性影响的数值分析

采用动量理论结合有限体积法求解N-S方程的方法,进行了螺旋桨滑流对带后缘襟翼机翼气动特性干扰的数值模拟分析.螺旋桨模型简化为一个激励盘,即一个无厚度可穿透的圆盘,以给定的圆盘拉力载荷,即压强差来模拟螺旋桨滑流的效果.采用嵌套网格技术处理桨盘与机翼间复杂的几何关系.算例计算表明,该计算方法能够较好的模拟滑流对带后缘襟翼机翼气动特性的影响.     

TRIP2.0软件非定常模块的开发与确认

在TRIP软件定常模块的基础上,采用双时间法和动网格技术开发了非定常模块.通过NACA0012翼型小攻角俯仰振荡和深失速两个算例的计算,对比分析了不同网格规模、湍流模型、时间步长等的模拟结果,并与试验数据做了比较,结果表明计算结果与试验数据基本一致,验证了TRIP2.0软件的非定常模块的正确性和精度.     

阵风响应及减缓的非定常数值模拟

基于网格速度和动网格方法,通过求解非定常N-S方程计算翼型副翼偏转时的阵风减缓响应,研究开环控制律参数和副翼缝隙对阵风减缓效果的影响。对于单个1-cos阵风,随着偏转角度的增大,减缓效果变大;由于粘性的影响,CFD计算得到的升力峰值减小量明显呈非线性的趋势,减缓效果低于准定常线性理论。对于连续sin阵风,通过调整副翼偏转的频率和最大偏转角,升力能达到很好的减缓效果,此时,副翼偏转频率应接近阵风频率,最大偏转角存在一个最优值。在阵风环境下,副翼间隙对升力系数影响不大,而力矩系数增大一倍,间隙会降低副翼偏转对升力的减缓效果。     

基于CFD数值模拟的拉瓦尔喷管流场分析

结合一维等熵管流理论,采用CFD方法对某型拉瓦尔喷管典型工作的流场进行整场数值模拟,阐述了网格生成,流场初始化,以及RNG湍流模型使用方法,并对拉瓦尔喷管工作特性进行了气动分析。研究表明,CFD计算结合一维管流理论能明显提高计算效率,数值计算结果与理论相符较好。     

临近空间桨梢小翼螺旋桨布局气动增效研究

基于雷诺平均Navier—Stokes控制方程和非结构网格旋彬静止滑移面技术,开展了临近空间螺旋桨桨梢小翼构型气动布局与增效作用机理研究,采用静止与旋转域的多块网格布局策略及转／静滑移面技术,发展了高效的螺旋桨流场与性能分析方法;通过开展上翘小翼、上下翘小翼等多种桨梢小翼螺旋桨布局研究,初步建立了基础参数的性能影响规律,研究表明：桨梢小翼螺旋桨布局是临近空间工况环境提高气动效率的一种合理可行的技术途径。     

用于变体翼梢小翼的伸缩栅格研究

翼梢小翼能抑制翼尖涡流形成、降低机翼的诱导阻力,翼梢小翼的高度是对减阻效果影响较大的参数之一.传统翼梢小翼仅针对巡航状态设计,而在起降、爬升等非设计状态减阻效果不佳.变体翼梢小翼能根据飞行状态主动改变几何外形和尺寸,实时优化减阻效果.为了实现变形,设计了一种用于变体翼梢小翼的伸缩栅格,通过步进电机驱动,可使翼梢小翼的高...     

A new grid deformation technology with high quality and robustness based on quaternion

Quality and robustness of grid deformation is of the most importance in the field of aircraft design, and grid in high quality is essential for improving the precision of numerical simulation. In order to maintain the orthogonality of deformed grid, the displacement of grid points is divided into rotational and translational parts in this paper, and inverse distance weighted interpolation is used to transfer the changing location from boundary grid to the spatial grid. Moreover, the deformation of rotational part is implemented in combination with the exponential space mapping that improves the certainty and stability of quaternion interpolation. Furthermore, the new grid deformation technique named ‘‘layering blend deformation' is built based on the basic quaternion technique, which combines the layering arithmetic with transfinite interpolation(TFI) technique. Then the proposed technique is applied in the movement of airfoil, parametric modeling, and the deformation of complex configuration, in which the robustness of grid quality is tested. The results show that the new method has the capacity to deal with the problems with large deformation, and the ‘‘layering blend deformation' improves the efficiency and quality of the basic quaternion deformation method significantly.     

基于子结构法的大型结构数值敏度计算技术

针对大型结构直接精确敏度分析方法解方程规模大、求解时间长的问题,发展一套基于有限元子结构法的处理技术。大型结构中,数值敏度计算的伪载荷仅与涉及设计变量的结构单元有关,因此存在大量零元的可能性。根据此特点通过节点重排将与设计变量有关的节点位移排到总位移列阵序的后面,按重排后的顺序投放刚度矩阵,然后对其进行区域分块,并聚缩得到规模较小的矩阵。算例计算表明:利用此矩阵求解结构位移场的导数,不但保持了精确法的精度,而且在优化设计的敏度计算中,这些较小规模的数据又被多次使用,从而显著提高计算效率。