基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。     

基于嵌套网格的CHN-T1标模数值模拟

为响应第一届航空CFD可信度研讨会(AeCW-1),对CHN-T1标模翼身-尾翼组合体构型采用两种RANS求解器进行了基于结构化嵌套网格的网格收敛性研究和抖振特性计算,并与对接网格结果进行了对比分析。参考AeCW-1网格生成指南自主研发了一族以约3.3倍规模增长的稀-中等-密-特密嵌套网格系列。相比于分区点对接网格,融合了贴体网格和笛卡尔网格优势的嵌套网格从根本上降低了网格生成难度,单块网格拓扑形式更为合理且正交性较好,同时能够有效地平衡近场和远场的网格量。对于跨音速运输机构型,同等规模下嵌套网格的物面网格密度较对接网格更大,能够对激波和分离等复杂流动进行更理想的预测。对比计算过程中,自研求解器OFS3D表现出了较高的嵌套网格计算效率和可靠性,而NASA CFL3D求解器的嵌套网格计算能力有待进一步挖掘和验证。计算还发现,不同湍流模型得到的表面压力分布、分离区大小以及最大升力系数存在差别,而风洞模型支撑和机翼静气弹变形对飞行器力矩特性预测影响较大。     

飞行仿真气动力数据机器学习建模方法

基于机器学习思想,提出了一种大空域、宽速域的气动力建模方法。该方法利用飞行仿真弹道数据辨识的气动力数据,采用人工神经网络技术,实现了对高度、速度、姿态和舵偏角等多维度强非线性特性的全弹道气动力数据的高精度逼近。首先,分析了神经网络层数、隐含层神经元个数等对建模误差的影响,通过对典型弹道气动数据的神经网络建模计算,确定了较合适的神经网络层数和较优的隐层神经元个数。进而,利用飞行仿真的弹道数据辨识出沿弹道的气动力,采用神经网络建立了包含多个弹道融合的气动力模型,输出量分别为三轴气动力系数和力矩系数。最后通过气动模型输出量与原样本数据的对比,以及4条未参与训练弹道气动数据的预测,验证了该气动力建模方法具有较高的精度。建模结果表明:采用神经网络方法建立的飞行器气动力模型,对拟合多源耦合输入全弹道非线性气动力是可行的和有效的,在样本覆盖的高度、速度、姿态和控制舵偏角范围内,气动力拟合能力较强,并具有一定的外推性。该项研究可以为基于飞行试验数据的气动建模提供新的方法,并且能为飞行器气动力数据挖掘、飞行仿真和总体性能分析提供参考。     

基于一种新的均匀化实施方法的FRP刚度预测

基于一种新的渐近均匀化(AH)实施方法,预测并讨论了单向纤维增强复合材料(FRP)的宏观等效弹性性能及纤维排列方式对等效力学性能的影响。该方法方便地将有限元分析(FEA)软件作为一个工具箱使用,只需在单胞上施加简单位移周期边界条件开展静力学分析,即可经过简单计算得到等效弹性常数,相比传统均匀化实施方法显著降低了实施难度并简化了计算过程。通过对比不同数值方法的结果验证了该方法的有效性和精确性。数值结果表明:六边形排列下单向纤维增强复合材料呈现横观各向同性,而正方形排列下则呈现宏观正交各向异性,经过刚度平均化过程可得到横观各向异性材料性质,纤维体积含量对两种排列方式下材料等效弹性模量影响显著但有所差别。     

固体火箭发动机零维两相燃烧室压强计算方法研究

为了更准确地预估含金属燃料固体火箭发动机的燃烧室压强,在压强计算中考虑两相流的影响,从一维两相喷管流动的求解出发,通过两相平衡流模型、两相常滞后模型、两相等温流模型、颗粒定温模型等模型的简化,分别推导不同模型下喷管中两相混合物的流量计算公式,再把流量公式应用到发动机零维内弹道理论中,推导并简化得到零维燃烧室平衡压强的计算公式。把压强公式用于HTPB推进剂固体火箭发动机和铝冰固体火箭发动机的燃烧室压强计算,结果表明,当固体推进剂中金属含量较高时(如铝含量为21%的HTPB推进剂发动机),用传统零维燃烧室压强公式预估的压强与实验误差较大,而使用合适的两相流模型和对应的零维燃烧室压强计算方法,在HTPB发动机中,能把压强预估结果与实验的误差降低到6%以内。如果使用多维内流场计算的方法,燃烧室压强预测结果的误差将下降到2.5%以内。结论发现在含金属固体火箭发动机的燃烧室压强计算中,考虑两相流的影响是必要的,而使用两相流修正后的零维燃烧室压强计算公式能够快速、较准确地预估这些发动机的燃烧室压强。     

颗粒无序分布的堆积床内部流动与传热分析

为了详细分析单组元发动机催化床内部的流动与传热,针对颗粒无序分布的堆积床,提出了一种网格生成途径,适合于对大量颗粒自然堆积的管道流动进行数值模拟。根据该方法,对不同颗粒直径、不同进口速度的管道堆积床进行了液体介质的流动与非定常传热模拟,并与实验结果和Ergun方程计算进行了对比。结果表明,本文提出的方法可以有效模拟颗粒无序分布的堆积床内部流动与传热,适合于大量颗粒的堆积床仿真;网格尺度小于1/20颗粒直径时,模拟结果与实验结果符合很好;在相同雷诺数下,摩擦系数随颗粒直径增大而减小,Ergun方程在高雷诺数下计算的流阻偏大;非定常传热时,壁面效应会使壁面附近温度先接近流体温度。     

基于CFD和CSM耦合的通用静气弹分析方法

提出了一种适用于有限元精细化建模的流固耦合插值点选择方法,通过RBF(径向基函数)方法实现流固耦合面的数据交换,实现了基于CFD/CSM(computational fluid dynamics/computational structural mechanics)耦合的通用非线性静气弹分析方法。以HIRENASD(high Reynolds number aero-structural dynamics)风洞试验模型为验证对象,数值结果很好地与风洞试验结构变形、气动压力分布吻合,验证了所发展非线性CFD/CSM耦合静气弹求解器的精度。详细研究了HIRENASD模型在大迎角(AOA)流动下的静气动弹性特性,以及该模型弹性变形对机翼气动特性影响规律。研究表明:HIRENASD弹性模型变形后其升力小于刚性模型;在小迎角范围内刚性、弹性模型升力差随迎角增大呈线性增长;当迎角大于4°后,升力差先减小后基本保持不变,呈非线性关系。      

一种部分森林覆盖环境电波传播特性预测方法

针对部分森林覆盖环境中传播媒介的非均匀性,提出一种基于动态网格的抛物方程(PE)预测方法,该方法采用宽角PE法进行大区域电波传播特性预测。为平衡PE求解的准确性与时效性,引入动态网格技术。通过仿真算例,探讨了基于动态网格的PE法在部分森林覆盖环境电波传播特性预测中的有效性。结果表明:在计算精度相似的情况下,相对于均匀细网格,动态网格会大大减少PE所需的计算时间;而相比均匀粗网格,动态网格具有更高的计算精度。     

共轴刚性双旋翼/机身干扰流场数值模拟

采用滑移网格技术求解Navier-Stokes (RANS)方程的方法,研究了共轴刚性双旋翼/机身的干扰问题。通过Caradonna-Tung旋翼、Robin直升机、Maryland直升机旋翼/机身干扰和Harrington 2共轴双旋翼等算例,验证了所提出的旋翼流场数值模拟方法的正确性。在此基础上,以Maryland机身为原型,分析了不同桨距的共轴刚性双旋翼与机身之间的干扰特性。结果表明:所提出的数值模拟方法能够很好地模拟共轴刚性双旋翼/机身的气动干扰特性;由于机身对于共轴刚性双旋翼下洗流场的阻滞作用,旋翼的悬停效率增加5%左右,并且随着拉力系数的增大使得悬停效率的增量更加明显;旋翼的悬停效率增加主要来源于下旋翼0°方位角附近的桨叶升力系数的增大,并且拉力系数的增量由桨根向桨尖方向逐渐减小。      

气压高度误差对适航标准的符合性方法研究

运输类飞机适航标准规定了气压高度误差的可接受范围,按照条款要求,表明符合性时需要将试验高度处的气压高度误差向海平面进行换算。相比ERJ、ARJ在适航审定试飞时采用的等静压误差处理方法,本文从是否考虑压缩性的角度出发,建立了适用性更高的气压高度误差转换方法,并通过某型飞机试飞结果分析进行了验证与说明。