基于N-S方程的飞机静气动弹性计算方法研究

摘要

本文分别基于传统的马蹄涡格网法以及N-S方程数值模拟弹性机翼的气动力，然后采用柔度法发展出了气动力作用下机翼静气动弹性特性的计算方法。 涡格法是升力面理论中一种比较实用的数值计算方法。它所采用的计算模型是：不仅沿展向分布离散的马蹄涡，在弦向也分布离散的马蹄涡，整个机翼用有限多个离散马蹄涡系来代替。控制点的诱导速度为附着涡和左右自由涡对控制点诱导速度之和。 在网格生成方面，采用代数方法及椭圆型方程优化相结合生成三维机翼的贴体网格，一种快捷的弹性变形技术用于生成弹性变形后的网格。 非定常渐近解法求解定常问题。基于Roe格式的有限体积法空间离散N-S方程、五步Runge-Kutta法进行显式时间推进。采用雷诺平均N-S方程模拟紊流，紊流模型选用SA一方程模型。应用了当地时间步长、焓阻尼、隐式残值光顺等加速收敛措施。 分别基于马蹄涡格网法和N-S方程，采用柔度法对机翼模型进行了考虑弹性影响的静气动弹性特性计算。不可压时两种方法计算得到的静气动弹性特性吻合较好，数值计算结果具有一定的实际参考价值。