17基于CFDCSD弱耦合旋翼悬停状态气动特性分析.docVIP

第二十八届（2012）全国直升机年会论文 基于CFD/CSD弱耦合的旋翼悬停状态气动特性分析 肖 宇 招启军 王俊毅 徐国华 （南京航空航天大学旋翼动力学国家级重点实验室，南京，210016） 摘要：本文在非惯性系下建立了一套高效高精度的适合于悬停旋翼的CFD/CSD弱耦合分析模型气动载荷通过基于嵌套网格的CFD方法计算求得，结构变形通过基于中等变形梁的有限元方法求解。为提高网格生成效率，采用基于网格点坐标转换的网格变形方法。CFD/CSD耦合分析计算了先进的UH-60A直升机悬停状态旋翼的表面压强及气动载荷，计算结果表明，与刚性旋翼CFD模拟结果比较，本文建立的CFD/CSD 耦合分析模型可以更准确的预估旋翼气动载荷。 关键词：，， 0 引言 在直升机旋翼气弹分析领域，旋翼气动力的准确预测一直是个难题，直升机综合分析软件中为了处理旋翼复杂气动力计算，常常对此进行简化，采用低精度的基于升力线理论的气动模型配合查表来进行气动分析，局限性较大，而CFD计算基于第一性原理，可以较为准确的计算各种飞行状态的气动载荷，因此越来越多的研究人员将之代入到直升机综合分析平台中，替换原有的气动力计算模块，提高了分析精度，得到了巨大的成功[1]。 1986年，Tung等人最早进行了耦合计算，其中CFD部分采用小扰动计算[2]，随后将之改进为基于全位势的CFD方法[3]，在这些计算中，都需要从CSD模块提供入流角，通过修正物面无穿透条件来代替网格变形，同时入流角也考虑了CFD计算域外尾迹的影响，采用该处理方式导致早期的CFD计算精度不足，在进行耦合时存在着收敛性问题[]。随着高速计算机的发展，基于Euler和N-S方程的CFD慢慢进入到耦合计算中，现在一般采用基于嵌套网格的CFD计算与CSD计算相耦合，可以直接模拟背景流场和所有的桨叶流场，避免早期耦合所需的入流角修正，如Potsdam等人[]，另外，Datta等人采用了基于单桨叶的混合CFD方法进行耦合计算[，对UH60A直升机旋翼的气动和结构载荷进行了预测，同样取得了成功，该方法虽然计算效率高，但是其中自由尾迹方法由于经验参数较多，无法对各类型旋翼做统一处理。 鉴于此，本文尝试将基于嵌套网格系统的C 1 旋翼CFD计算模型 1.1 计算网格 计算网格采用嵌套网格系统，如图1所示。整个网格系统分为两大部分：一是包含弹性桨叶的贴体网格，二是包含整个计算流场域的背景网格，桨叶贴体网格采用C-H型网格，用来计及旋翼桨叶的刚性运动及弹性变形，背景网格则采用圆柱形网格。 图1 计算网格系统 由于桨叶网格需要扑捉物面的粘性流动，因此使用N-S方程求解，而背景网格由于一般只需扑捉旋翼尾迹，从计算效率上考虑，采用Euler方程求解。为了与CSD相耦合，桨叶网格需要进行动态变形，弹簧法虽然生成网格质量较高，但是计算需要迭代，耗时较长，因此本文采用基于坐标转换的网格变形[7]。该方法变形后网格可直接通过对变形前网格坐标的矩阵运算得到，无需迭代，因此具有较高的效率，便于CFD/CSD的耦合运算。 1.2 控制方程 针对直升机悬停状态，在旋转坐标系下建立控制方程，此时流场可视为定常状态，具体形式可以表达为 (1) 式中对流通量，Jameson中心差分格式，时间离散格式则采用显式Runge-Kutta积分求解。另外采取了适当加速措施，如当地时间步长及隐式残值光顺等。 2 旋翼CSD计算模型 在直升机综合分析平台中，旋翼CSD计算模型一般采用二维线性截面分析结合一维非线性梁分析构成，如图2所示，这是为了在保证分析精度的同时，避免直接处理三维非线性有限元而造成的时间和内存的浪费。本文同样选用该方法构建旋翼结构计算模型。 = + 图2 旋翼CSD计算模型分解示意图 2.1 截面分析模型 本文截面分析部分采用基于二维有限元的分析方法，与早期解析法相比，有限元法能处理各种截面且无需做任何假设，因此得到广泛的应用。本文有限元单元采用四边形单元和三角形单元两种，分别采用Gauss积分和Hammer积分进行处理。图3显示了本文截面分析的有限元网格示意图，图4是后缘混合网格示意图。 图3 截面示意图 图4 后缘混合单元示意图 本文选用基于中等变形梁理论的改进的YF模型进行一维梁分析[8]。原始的YF模型为23自由度，采用了基于圣维南原理的截面分析模块，具有一定的局限性。本文在此基础上，对YF模型的自由度进行了精简改进，针对欧拉梁和铁木辛柯梁，自由度数分别缩减为14和20，从而更有针对性的处理桨叶模型。图5和图6为本文采用的梁模型自由度示意图。 图5 欧拉梁模型14自由度示意图 图6 铁木辛柯梁20自由度示意图