基于cfd技术的swim尾涡气动特性分析.docxVIP

基于cfd技术的swim尾涡气动特性分析 0 气动响应特性仿真模型 飞机通常受到强烈的干扰，在导航和机动性飞机上行驶时，它们通常受到严重的风吹影响。阵风产生确定的或随机的附加气动载荷,引起飞行器的刚体运动,过大幅值的阵风可能诱发飞行器出现严重偏离航线、大迎角失速、剧烈翻滚等问题,从而影响飞行器的飞行性能与安全。为了研究阵风附加的气动载荷与飞行器之间的相互作用,通过CFD技术求解飞行器遭遇阵风时的气动响应,同时求解飞行器飞行的六自由度方程, 实现CFD-6DOF(Degree of Freedom)的耦合来研究飞行器阵风中的响应特性。 早期对飞行器阵风响应的研究主要采用准定常假设法,即飞行器在某一时刻的阵风响应载荷等于由阵风引起的攻角变化乘以当地定常升力系数, 该方法使用方便快捷但计算精度相对较低 尾涡可近似认为是垂直阵风与侧向阵风的叠加,建立尾涡的阵风场模型即可研究尾涡对飞机飞行所带来的影响。当小型飞机进入大型飞机的尾涡时,尤其是在飞机的起降阶段,小型飞机可能出现剧烈地翻滚而导致失事,而飞机在某些特定状况下不可避免地进入其他飞机的尾涡场,因此有必要对飞机在尾涡中的响应特性进行研究。 本文首先采用改进的Lamb-Oseen涡模型建立尾涡形式的阵风场;然后将尾涡引起的扰动速度添加到N-S方程的求 解中,对SWIM (Subsonic Wall Interference Model)尾涡场中的定常气动特性进行验证;最后通过CFD-6DOF的耦合数值模拟SWIM在尾涡场中的飞行轨迹,对飞行过程中出现的抖动、下沉、滚转等动态特性进行研究。 1 cfd-6dof矩阵响应 1.1 双时间步长法q 本文使用的流场求解器是基于非结构混合网格雷诺时均N-S方程求解程序。流动控制方程 式中:守恒变量Q = [ρ,ρu,ρv,ρw,ρE] 对于方程(1),采用有限体积法对其进行数值求解,非定常时间导数项的离散采用双时间步长法,在虚拟时间步中的推进计算采用隐式LUSGS近似法,空间对流项采用Roe的通量差分来离散, 黏性项则采用中心格式离散 1.2 翼型迎角突变作用 假设翼型以一定速度水平飞行,突然受到一上升气流作用 (如图1所示 由于直接赋予翼型迎角突变的边界条件的缺陷,目前阵风响应的数值模拟一般采用“网格速度”思想:根据相对运动关系,若网格的速度为-W 1.3 惯性系下质心平动方程空间位置 飞行器作为刚体在空中运动,一般有六个自由度。相应的有六个动力学 方程,其中三个描述平动,三个描述绕质心转动。另外,还有六个运动学方程,确定飞行器空间位置和姿态。飞行器坐标示意图如图2所示。x轴指向地平面某选定方向,z轴铅垂向下,y轴垂直xz平面,按右手定 则确定方向。 惯性系下质心平动方程 体轴系下绕质心转动方程 空间位置 空间姿态 式中:m为质量;v为速度;F为力;ω为角速度;M为外力矩;I为惯性矩;通常ω 在求解刚体运动方程(2)~方程(5)时,采用的线性四步Adams预估校正法的步骤为:首先求解n时刻的流场,得到气动力并传给预估步刚体运动方程,求解n+1时刻飞行器空间位置和姿态,再进行流场计算,得到气动力并传给校正步刚体运动方程得到修正后的n+1时刻飞行状态量。如此反复,直到仿真结束。因阵风对飞行器飞行速度和角速度带来较大影响,进而引起飞行迎角和侧滑角变化,本文通过n+1时刻质心位移增量和角度增量的变化来考 虑速度和 角速度变 化对气动 特性的影响。 2 cfd计算结果 为了研究 飞行器尾 涡中的气 动特性,R. Struijs等 采用改进Lamb-Oseen涡表达式,使之更符合尾涡特性,改进后的Lamb-Oseen涡表达式为 式中:-e SWIM在涡旋场中被 放置于多 个位置,选取其中五个位置(如图3所示)来捕捉模型在涡旋场中升力和滚转的极限状况,分别对五个不同位置的气动参数进行CFD计算。F29模型产生的翼尖涡参数及位置如表1所示。CFD计算结果与风洞实验结果、文献 [10]CFD计算结果 的对比,如图4所示。 从图4可以看出:SWIM机翼采用NACA4412翼型,在0°迎角下的升力系数C 对实际飞行而言,飞机进入尾涡时,滚转响应特性尤为重要。当SWIM处于y=-0.42m时, 由于位于左舷涡旋中心,因此会承受最剧烈的滚转效应。本文和文献中的CFD计算所得滚转力矩系数与实验值整体趋势一致,大部分数值点符合较好 (如图5所示),表明本文所使用的尾涡模型能够较准确地预测飞机在尾涡中的滚转响应特性。 3 swim的速度特性 为了研究尾涡对飞机飞行的危害性,数值模拟SWIM从右舷涡旋正上方10m处、以-5°的俯仰角俯冲进入尾涡场,迎角、侧滑角均为0,如图6所示。取非定常计算步长为0.005s,总模拟时间为4s,其中0~1.5s的内迭代 步数为150,1.5~