国家数值风洞工程结构网格变形程序SGDP V1.0开发与应用

对结构网格变形程序SGDP V1.0的理论方法、组成架构、主要功能、运行流程、运行性能、变形算例和工程应用进行了简要的介绍。SGDP V1.0是国家数值风洞NNW项目共性基础库网格变形模块下实现的第1款模块化程序,在吸收前期研究代码的基础上,通过代码重构、接口规范、功能扩展等技术实现开发。目前,SGDP V1.0程序依托RBF-TFI方法实现结构网格变形,采用模块化的代码结构和扁平化的参数设计,可以实现4种不同类型运动下的变形网格生成。SGDP V1.0程序每百万节点网格的读取时间约为0.01s、写出时间约为0.008s、网格更新时间约为0.1s,可在数秒内实现千万量级网格的变形更新,具有非常高的运行效率。网格变形算例和工程应用结果表明,SGDP V1.0程序能够处理复杂拓扑网格的变形问题,且有很好的网格质量维持能力,从而能够获得更加精确的数值模拟数据。     

含缺陷铝蜂窝道面板弯曲和疲劳特性

为了研究面板厚度及不同类别、不同程度的先天缺陷对铝蜂窝夹芯结构道面板疲劳性能的影响,对7组铝蜂窝夹芯道面板试样进行三点弯曲静载试验和三点弯曲疲劳试验。结果表明:面芯脱焊缺陷对铝蜂窝道面板的疲劳性能影响最大,含10%和30%面芯脱焊缺陷的试件疲劳性能分别降低了76%和92%,说明防止面芯脱焊缺陷至关重要。经分析,铝蜂窝道面板弯曲疲劳特性服从威布尔分布,相关系数超过0.95以上。对完好试样和含不同类别、不同程度先天缺陷的试样分别建立了反映应力水平和疲劳性能关系的威布尔方程。     

风洞试验中模型迎角视觉测量技术研究

 发展了一种基于双相机的风洞试验模型迎角视觉测量技术。即在风洞模型表面绘制或喷涂一些高对比度的标记点,采用两台400万像素的工业相机采集模型的运动图像,然后求解图像的共线方程获得模型表面标记点的三维空间坐标,利用表面标记点坐标最小二乘拟合出一条空间直线,通过不同时刻空间直线单位方向矢量的变化计算出模型的迎角变化。利用精密电子倾斜仪对迎角视觉测量系统进行静态标定,结果表明:迎角视觉测量系统的准确度在0.01°以内。迎角视觉测量系统在2 m超声速风洞中对一尾支撑无主测力天平铰链力矩测量模型的实时迎角进行了测量,结果表明:系统具有很好的响应特性,可以作为风洞试验中迎角测量的一种有效方法,尤其是用于无主测力天平时模型真实迎角和支杆弹性角的测量。     

有无涵道动力下的类BWB低速布局气动特性研究

随着航空飞行器的不断发展,新型非常规气动布局研究早已成为世界航空大国的关注焦点。基于TRIP 3.0软件平台通过等效盘模型应用到涵道风扇内外流一体化模拟中,完成对类BWB低速布局有无涵道风扇动力下的气动特性和流场影响分析。首先,对某单独螺旋桨验证算例进行分析;然后,对单独涵道风扇进行计算,得到设计推力下的涵道动力叶片数、安装角及转速等参数;最后,将设计推力下涵道动力参数应用到类BWB低速布局全机有无动力模拟中。结果表明:计算拉力和扭矩值与试验值吻合较好;尾部布置的涵道动力对机身后段及尾翼压力分布影响显著;涵道对前方气流抽吸作用,加速机身后段上表面气流流速,减小内侧平尾有效迎角,对机身和平尾升力以及俯仰力矩系数影响较大;在小迎角状态下,涵道动力对机身增升作用明显,会产生明显低头力矩,但对平尾作用正好相反,且两者增量大小相当,使得全机增升效果不显著,且俯仰力矩变化较小。     

DLR-F6/FX2B翼身组合体构型高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6和DLR-F6_FX2B 2种翼身组合体构型的高阶精度数值模拟,计算外形来自AIAA第三届阻力预测研讨会。主要目的是确认WCNS模拟跨声速典型运输机构型和预测局部构型变化引起的气动特性变化量的能力。在固定升力系数条件下,采用粗、中、细3套网格开展了网格收敛性研究,从气动力系数、压力系数分布、表面流态等方面研究了网格规模对DLR-F6和DLR-F6_FX2B翼身组合体数值模拟结果的影响;采用中等网格开展了来流迎角对2种翼身组合体气动特性的影响研究。通过与National Transonic Facility(NTF)的试验结果和CFL3D的计算结果对比,表明采用高阶精度计算方法得到了网格收敛的数值模拟结果,较好地模拟了DLR-F6翼身组合体局部修型引起的微小气动特性变化和翼身结合部流动特性的差异。     

CRM翼/身/平尾组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼/身/平尾组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测方面的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细3套网格,对应的y⁺分别为1、2/3和4/9。通过对CRM模型的计算和分析,研究了网格规模对气动特性、压力分布和翼根后缘局部分离区的影响。通过与DPW4统计结果和部分实验结果的对比,高阶精度数值模拟结果表明,本文的气动特性计算结果与统计平均结果吻合较好;网格密度对机翼内侧的展向站位压力分布影响较小,对机翼外侧展向压力分布影响较大;网格密度对翼根后缘局部分离区略有影响。     

中断起飞极限速度的计算分析

提出了计算中断起飞极限速度的基本公式和程序设计方法,根据飞行部门中断起飞的实际问题,对两种情况的中断起飞过程进行了计算分析。计算方法和过程分析对于飞行中处置中断起飞特情具有指导意义。     

CRM-WBPN风洞试验模型数值模拟

基于雷诺平均Navier-Stokes（RANS）方程和多块对接结构网格技术，采用亚跨超CFD软件平台（TRIP 3.0），研究了共同研究模型（CRM）风洞模型支撑装置对CRM翼/身/架/舱（CRM-WBPN）模型压力分布及气动特性的影响。构造了一个包含小、粗、中、细网格的4套网格，并进行了网格收敛性研究，在升力系数为0.50条件下，开展了固定升力系数下模型支撑对CRM-WBPN模型和CRM-WB模型压力分布及挂架短舱阻力增量的影响，以及不同来流迎角下模型支撑对CRM-WBPN模型气动特性的影响。通过与NASA国家跨声速风洞设备（NTF）试验结果和CRM-WB模型数值模拟结果的对比分析，表明模型支撑对固定升力系数下挂架短舱阻力增量影响较小；在计算构型中考虑模型支撑的影响，引起升力系数和阻力系数下降、俯仰力矩系数增加，数值模拟结果更加接近相应的试验结果。     

风洞试验中模型迎角的视频测量及精度研究

模型迎角的准确测量是获得高精度风洞试验数据的基础（如失速点与最小阻力点的准确迎角值是飞行器研制与改型的关键数据）。为此,提出风洞试验中模型迎角的视频测量方法,分析其测量精度。2m超声速风洞中的多个迎角视频测量实例表明：该方法的迎角测量精度≤0．01°,可清晰分辨幅度为0．01。的迎角振动过程,而且该方法既不破坏模型的外形,又不改变模型的刚度与强度,测量精度不受模型振动影响,具有实用价值。     

湍流模型离散精度对数值模拟影响的计算分析

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS),通过改变物面法向第一层网格间距,开展了剪切应力输运(SST)两方程模型不同离散精度的数值分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR 7301两段翼型和高速RAE2822翼型,研究内容主要包括湍流模型的二阶精度离散和五阶精度离散两种方式对收敛历程、边界层湍流黏性系数分布、边界层速度分布、压力系数分布以及气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于不同的第一层物面法向网格间距,湍流模型离散精度对低速绕流计算结果有比较明显的影响,对于高速小迎角附着流动计算结果影响不明显;相对于湍流模型二阶精度离散,湍流模型高阶精度离散网格敏感性较弱,具有更高的数值模拟精度,但收敛性略差。