基于流固耦合的机翼颤振数值分析.docxVIP

基于流固耦合的机翼颤振数值分析 1 物理模型的建立 气流流经飞机，导致不确定的气动力，导致飞机的振动。颤振是造成机翼失效的重要原因之一,深入研究机翼颤振问题对设计高性能的机翼是非常重要的。 随着计算机技术的发展,人们着手用数值方法研究机翼的颤振问题。数值模拟机翼颤振问题的难点一方面在于要在计算流体力学(CFD)程序和计算固体力学(CSD)程序之间传递数据,而这两者的程序在算法上有很大区别。另一个难点在于机翼的颤振经常伴随着复杂的分离流动,过去流固耦合一般都采用简化的气动模型,比如在流体区域采用欧拉方程等,采用这些简化的物理模型是很难准确地模拟失速颤振时机翼周围复杂的流场的。 本文采用更准确的流固耦合数学物理模型研究机翼的颤振问题:在流体区域求解非定常的N avierStokes方程和低雷诺数双方程q-ω湍流模型,在固体区域采用两个自由度的PAPA模型(a pitch and plunge airofoil),求解两个自由度的振动方程;每计算一个时间步长后,非定常流场和固体弹性系统都相互传递一次边界条件;计算流场时采用高精度、高收敛率的LU-SGS-GE隐式格式和高阶MUSCL TVD格式)。应用上述方法本文成功地实现了包括失速颤振在内的机翼颤振问题的流固耦合数值模拟,并得出机翼颤振问题的一系列规律。 2 机翼的控制方程 气流流经机翼时产生的非定常气动力引起的机翼振动形式有扭转振动和弯曲振动两种。对于二维问题,机翼的颤振问题可以用PAPA模型简化,即认为在振动中心处有两个弹簧支撑机翼,一个控制机翼的扭转运动(α),另一个弹簧控制机翼的上下运动(φ),(如图1所示)。控制方程包括流体的基本方程和固体的基本方程两部分。 流体部分的控制方程为Farve平均Navierstokes方程和低雷诺数双方程q-ω湍流模型,任意曲线坐标系ξi下控制方程的无量纲形式可写成: 振动机翼的绕流流场为典型的湍流流场,选择合适的湍流模型对于准确模拟振动机翼的流场至关重要。本文采用的低雷诺数双方程q-ω模型能够较准确地模拟大尺度分离流,更有着计算量少,边界条件容易处理,又适应粗糙的初始湍流流场的优点。固体部分的方程为: 3 时间步长的确定 本文采用流固耦合方法研究机翼的颤振问题:在流体区域用计算流体力学的方法(CFD)求解非定常Navier-Stokes方程,在固体区域用计算固体力学的方法(CSD)求解非定常的振动方程。在每一个时间步长结束后,流体和固体都相互交换一次边界条件(流体把计算出来的气动力传递给固体,固体把计算出来的固体边界位置和速度传递给流体),然后进行下。个时间步的计算。 Navier-Stokes方程的数值计算采用了LU-SGS-GE隐式格式和改良型四阶精度MUSCL TVD格式处理。因为机翼的位置随时间不断变化,因此每一个时间步长后都需要生成一次新的网格。经过大量数值计算比较,采用的网格点数为361×76。最靠近壁面的一排网格点到壁面的垂直距离为0.00002c,其中c为弦长。振动方程采用4阶Rounge-Kutta法来求解,具体解法参见数学手册。 4 结果和分析 4.1 其他湍流模型计算结果比较 准确计算失速颤振状态下非定常流场是研究失速颤振的一个难点。为了验证程序正确性,首先计算了NACA0015振动机翼深度动态失速时的非定常二维流场,并和美国NASA Ames研究中心Piziali的实验结果以及其他四种湍流模型的计算结果进行了比较,计算中用到的流动参数完全按照实验给定:来流马赫数Ma=0.29,雷诺数为1.95×106(特征长度为机翼的弦长c,特征速度为入口来流速度)。攻角与无量纲时间t的关系为α=α0+α1sin(2kt),平均攻角α0=15°,α1=4.2°,衰减频率系数k=0.1,t和有量纲t*的关系为,。非定常计算中采用的无量纲时间步长Δt=0.0108,计算一个周期需要10000个时间步长。 从图2可看出,采用本文的算法得到的数值模拟结果和试验吻合很好,尤其在振动的下冲程,本文的结果比以前文献中的结果精度更高。详细的计算结果参见文献。 4.2 失速颤振的频率场 以NACA 0015机翼为例用流固耦合的方法研究了大攻角(α0=0°～50°)范围的颤振,气动参数和结构参数如下:M∞=0.29,Re=1.95×106,Sα=0.2,Iα=0.3,ωφ=0.06,ωα=0.06,阻尼系数ζα=ζφ=0.0,U*=1.83。 当来流攻角α0=0°时,机翼周围没有分离流产生,振动达到平衡状态后,振动的频率谱存在一个明显的主频(ω1=0.012),而且这个主频与机翼的固有频率(ωα=ωα=0.06)相差很远(图3(a),图3中|A|为α的幅值)。0°攻角下振动的极限环是一个封闭的周期性曲线(图4(a)),这时的振动属于稳定的强迫振动。