CFD结合降阶模型预测阵风响应.pdf

维普资讯 第 40卷 第 2期 力 学 学 报 Vo1．40，No．2 2008年 3月 ChineseJournalofTheoreticalandApplied Mechanics Mar．，2008 CFD 结合降阶模型预测阵风响应 杨国伟 。)王济康 (中国科学院力学研究所高温气体动力学重点实验室，北京 100080) 摘要 传统的阵风响应主要在频域内进行分析，气动载荷基于线性方法计算，不能考虑黏性和跨声速流动影 响．飞机设计需考虑不同频率和不同形状阵风的响应，基于 CFD的阵风响应预测由于计算工况太多，工作量 巨大．本文发展了一种 CFD结合非定常气动力 ARMA (autoregressive—moving—averagemode1)降阶模型的 阵风响应分析方法，CFD只要针对给定频率和形状的一种阵风响应进行计算，对获得的气动力时间历程运用 线性最小二乘法参数辨识ARMA降阶模型的系数，则对任意频率和形状的阵风，代‘入降阶模型即可确定该阵 风的响应，大大提高了计算效率．为效验发展的方法，先计算NACA0012翼型在低马赫数 0．11的阵风响应， 通过对比CFD、ARMA 降阶模型及早期发展的不可压阵风响应预测方法的结果，验证了方法的有效性．再 对比CFD、ARMA在跨声速马赫数 0．8的阵风响应预测结果，证实所发展的方法对跨声速阵风响应预测亦 是有效的． 关键词 CFD，阵风响应，ARMA 降阶模型，参数辨识，最小二乘法 中图分类号： 0355，V211 文献标识码：A 文章编号： 0459—1879(2008)02—0145—09 引 言 应计算的缺点是计算量巨大，不适合于飞行器研制 中需对不同频率和不同形状阵风响应分析的需要， 阵风响应研究飞行器通过大气阵风引起的响应 也难于与气动伺服弹性结合研究飞行器的阵风减缓 问题．阵风响应分为时间和频域两种分析方法．传 设计． 统的频域阵风响应分析方法是将气动弹性运动方程 近 10年来，为了开展气动弹性定性分析和参 在广义坐标系下写为频域形式，输入不同频率的谐 数设计，基于 CFD 的非定常气动力降阶模型 (re— 频阵风，其引起气动力通过修正频域广义空气动力 ducedordermodel，ROM)得到了长足发展，其方 影响系数加以考虑．时间域分析方法在早期针对不 法主要分为本征正交分解法 (properorthogonalde— 可压流动的二维翼型发展了Wagner和 Kussner解 composition，POD)技术 基于系统辨识技术 L5_ 析函数方法 Wagner方法用于预测翼型攻角阶跃 和基于谐波平衡 (harmonicbalance，HB)技术 L6】发 变化的升力时间变化历程， Kussner方法预测阶跃 展的各种非定常气动力降阶模型．本文选择基于系 阵风的升力时间历程．任意形状的阵风响应可以通 统辨识的ARMA模型进行阵风响应研究，主要研 过Duhamel积分获得的卷积公式计算．后来通过对 究 CFD结合该模型模拟阵风响应的可行性，因此， Kussner函数的系数进行经验修改 2【J，使该方法也可 选择 NACA0012翼型作为算例，通过 CFD直接模 用于可压缩流动的阵风响应分析．其缺点是该方法 拟、ARMA模型及其不可压阵风响应解析解的相互 只对二维平板不可压流动有效． 比较，验证本文所发展的方法． 作者在 2004年率先开展了基于CFD的飞行器