自由弹性机翼跨音速阵风响应分析.docxVIP

自由弹性机翼跨音速阵风响应分析 基于cfd的阵风响应分析 传统的矩阵响应计算在频域进行，主要采用线性平面气动力，并结合动态测量变量的弹性方程对矩阵效应进行分析。 随着代理模型的发展, 出现了基于CFD的降阶模型方法 1 气动弹性方法和公式 1.1 气动力f的传统描述 机翼的结构运动方程可用矩阵形式描述 式中M为广义质量矩阵;G为广义阻尼矩阵;K为广义刚度矩阵;q为广义坐标向量, 右端项的广义气动力包括阵风引起的广义气动力F 基于线性假设, 采用传统振动模态的方法描述机翼表面的振动变形, 引入状态变量E 将式 (2) 改写为线性方程组 式中0为n阶零方阵, I为n阶单位阵。 采用经典的4步Runge-Kutta法求解线性方程组 (3) : 假设机翼的弹性力在线性范围内, 采用传统振动模态的概念来描述机翼表面点的纵向振动变形, 将结构运动方程从二阶常微分方程改写为线性方程组的形式, 采用经典的4步Runge-Kutta法求解线性方程组, 采用面样条插值方法实现模态插值。 1.2 时间步长及对流项格式 对于带运动边界的非定常流体运动, 本文采用Euler方程的任意拉格朗日欧拉公式 (ALE) , 它的微分形式为 式中W为守恒变量矢量;f, g和q为对流通量矢量。本文采用有限体积法格心格式, 方程采用双时间步长推进方法, 在每个物理时间步内进行虚拟定常迭代时, 采用5步Runge-Kutta显式时间推进。对流项采用带人工耗散项的Jameson中心格式。初始条件取均匀流场, 壁面处采用无穿透边界条件, 远场边界处采用无反射边界条件。 1.3 冲突上的流场模拟 为了计算使翼型攻角发生突变的阵风以及后面的陡边型阵风的响应, 引入“网格速度”概念 式中u, v和w为直角坐标系下的速度分量;x 即 上式就相当于给网格一个向下的网格速度, 由此来模拟向上的流场, 即迎角的突变。 如图1所示, 对于不同的阵风, 在某一时刻t, 赋予位于x处的网格和阵风速度相反的速度w 式中f (x, t) 对应于不同阵风的数学表达式。 1.4 流固耦合法 求解时采用传统的松耦合迭代过程 2 自由弹性机翼在资本特性和阵风作用下的响应 传统的卷积模型首先计算刚性机翼在锐边阵风下的响应, 任意阵风引起的广义气动力响应 (不考虑结构) 可用以下的卷积公式得到 式中 本文发展的方法通过求解方程 (1) 直接计算自由弹性机翼在锐边阵风下的响应, 任意阵风的响应 (包括结构) 可用以下的卷积公式得到: 式中 3 计算和分析 3.1 cfd方法的应用 本算例计算喷气运输机机翼 (BAH Wing) 在方波阵风下的响应过程, 这是MSC.Nastran计算机翼离散阵风响应的标准算例, 文献[13]中有具体的结构建模方法。结构提取6阶模态, 其中第1阶是刚体沉浮模态, 其余5阶是弹性模态, 参考弦长c=4.128m, 半展长b/2=12.7m 采用CFD方法计算机翼在方波阵风下的广义坐标响应, 并和Nastran的结果比较见图3和4所示。 从图3和4可以看出, 本文CFD计算的结果和Nastran计算的结果一致, 验证了本文所采用的CFD方法的正确性。 3.2 cfd方法的计量结果 计算相同机翼 (BAH Wing) 在锐边阵风下的响应过程。计算马赫数Ma=0.85, 阵风幅值取0.05倍的来流速度。用CFD方法计算得到的升力系数, 广义坐标响应曲线如图5～7所示。 从图5～7可以看出, 自由 (沉浮) 弹性机翼在锐边阵风的作用下会一直向上运动, 最后保持匀速直线运动, 运动速度等于阵风的速度。 3.3 卷积模型和cfd模拟的阵风响应分析 本算例考虑相同机翼 (BAH Wing) 在1-cos阵风下的响应过程。计算马赫数Ma=0.85, 阵风形状如图8所示。阵风幅值取0.05倍的来流速度, 阵风宽度2 H=12。使用3.2中的结果和公式 (11) 及 (12) 建立卷积模型, 计算得到的升力系数, 广义坐标响应曲线和翼尖加速度响应曲线和直接CFD模拟的结果比较见图9～11所示。阵风幅值取0.1倍的来流速度, 其他条件不变, 采用卷积模型和直接CFD模拟的广义坐标响应曲线和翼尖加速度响应曲线比较见图12～14所示。 从图9和12可以看出, 卷积模型和CFD计算的升力系数峰值吻合, 在机翼穿过阵风区域后, 两者存在很小的差别。这是因为自由弹性机翼在锐边阵风下作用下, 最终会向一个方向匀速运动, 如图6, 这会导致网格变形较大, 0.5s以后计算得到的锐边阵风响应值会有一定的误差。和升力系数一样, 图10, 11, 13和14中的广义坐标也会有同样的问题, 刚体模态广义坐标响应曲线的差别比弹性模态广义坐标响应曲线明显, 但是这种差别不大。 从图9～11和图12～14可以看出, 对于不同的