基于非结构网格数值模拟的气动载荷计算的研究.pdfVIP

粥 巷 第二届全国低跨超声速空气动力学交流会论文集 基于非结构网格数值模拟的 气动载荷计算研究 王世安 周喜军 尹贵鲁 (贵州飞机设计石Ji究所贵州安顺 561000) 摘要 获取飞机气动载荷通常蔓进行全机测力、部件测力、全机测压等实验，通过实验 结果，对飞行参数的仿真模拟，对实验结果的插值、修正、积分等．晟终获得气动载荷，i仃往 飞机的初步设计阶段和经费缺乏的情况下小训能进行大量的实验，但需要计算飞机的载荷平¨寿 命。笔者利用基于非结构嘲格的数值模拟方法来确定飞机的气动载荷，且对于目前不太容易计 算准确的偏转舵面的情况的气动载荷也进行了数值计算，将计算结果与部件测力和全tJl,钡9压结 果进行了比较，收到较好效果。 关键词 非结构恻格 全机测力 全机测压 部件测力 法向载荷系数 压力分布 1 引 言 在飞机型号的研制中，计算气动载荷一般先要进行全机测力实验，利用实验结果进行 飞行参数的仿真模拟，进而进行全机测压实验，根据飞行参数对全机测压实验结果进行插 值和积分，最后选取严重载荷情况。在计算中，由于对全机测压实验结果的插值和积分会 带来一定的误差，这就需要通过部件测力实验来校准载荷的总体量级，要对实验结果进行 合理的修正。因此，气动载荷的获取是需要花费大量时间和资金的，而这在飞机设计初期 和研制经费缺乏的情况下，通过多项实验来获取气动载荷是非常困难的。 近年来，随着计算机技术的进步，CFD也获得了快速的发展，据报道，在F．22飞机 的气动载荷计算中，采用了CFD的TEAM程序，在生产型F．22飞机的风洞试验之前，用 TEAM给出了全机气动载荷进行初步设计，加快了进度。在双座飞机的设计中，利用TEAM 计算的单座和政座飞机载荷的差值修正单座的试验数据，得出双座的气动载萄。通过CFD 数值模拟来获取气动载荷可以节省大量的经费开支，囡为通过一次数值模拟可以获得晰如 全机洲力实验、全机测压实验、部件测力实验的结果，而结果是协调一致的，无须修』L处 理。笔者对具有复杂外形的某型机生成了非结构网格，利用有限体积法求解了全机气动力， 征，W号研制中对利删数值模拟确定气动载荷的方法进行了尝试。 2网格生成 篓三旦釜堡!堡堕箜苎j堡墨塑塑鲎茎塑垒!生墨篓 用于复杂外形网格生成的方法不少，分区网格有较好精度，但对复杂外形生成分区网 格比较凼难，且闷格质晕分布不均的话影响计算精度；非结构网格对复杂外形较容易生成， 不需要分区，有较好的1+程计算精度。笔者利用基于Octree技术的方法对某型机生成了非 结构网格(如图1、图2所示)，其生成方法是采用自顼向下的生成方法，首先给定网格最 火尺度和几何外形网格尺度，在几何外形刷同生成一个规则的几何面把几何外形围住，在 几何面中按网格最大尺度生成四面体网格，接着根据几何外形进一步划分网格元素，如果 几何外形某部分的尺度比给定的网格最大尺度小，则穿过该部分的网格元素应进行划分， 同样对于几何外形上的线按给定尺度划分网格元素，穿过线的网格尺度应比穿过面的网格 圈1未偏舵表面网格 图2偏舵时表面网格 尺度小。在完成网格的划分之后，就要切割网格元素以便拟合真实几何外形，切割的基本 方法是将穿过几何外形的面、线、点的网格的边进行细分或将边上距几何外形最近的结点 移到几何外形的面、线、点上。拟合的顺序是先拟合几何外形的面，然后是几何外形的线， 最后是几何外形的点，因此，为了更好的捕捉几何外形，在几何外形的关键都位必须给出 线和点的描述。 在网格生成之后，为了进一步提高网格质量，可以采用诸如Laplace之类的方法对网 格进行光顺处理。 3计算方法 3 1控制方程 由于在歼击系列飞机的计算中，粘性的作_|=lj仅在离机体很薄的附面层内，对气动载荷 的确定不会引起太人的误差，为了加快计算进度，笔者采用Euler方程作为控制方程。 ．二维直角坐栎系。F守恒7理Euler方程可以表达为 望生+兰生+生￡+旦旦：0 it ax av 52 其中p为守匣变量，E、F、G为对流通量，