基于CFD技术的高效叶栅耦合颤振分析方法.docxVIP

航 空 学 报 Ｄｅｃ．２５ ２０１４ Ｖｏｌ．３５ Ｎｏ．１２ ３２３２－３２４３ ＡｃｔａＡｅｒｏｎａｕｔｉｃａｅｔＡｓｔｒｏｎａｕｔｉｃａＳｉｎｉｃａ ＩＳＳＮ １０００－６８９３ ＣＮ１１－１９２９／Ｖ ｈｔｔｐ：／／ｈｋｘｂ．ｂｕａａ．ｅｄｕ．ｃｎ ｈｋｘｂ＠ｂｕａａ．ｅｄｕ．ｃｎ ｄｏｉ：１０．７５２７／Ｓ１０００－６８９３．２０１４．０７７ 基于 ＣＦＤ 技术的高效叶栅耦合颤振分析方法 苏丹１，张伟伟１，＊ ，全金楼１，马明生２，叶正寅１ １．西北工业大学 翼型叶栅空气动力学国家重点实验室，陕西 西安 ７１００７２ ２．中国空气动力研究与发展中心 计算空气动力研究所，四川 绵阳 ６２１０００ 摘 要：为了考虑叶轮机叶片结构与流体之间的耦合效应 ，同时提高叶轮机颤振数值 研究的效率 ，发展了一种基于非定 常气动力降阶模型（ＲＯＭ）的叶栅耦合颤振分 析方 法。 该方法运用时域计算流体力学 （ＣＦＤ）技术计算少数几个叶片的 非定常气动力 ，通过系统辨识及一些假设构建整个叶栅振动的非定常气动力降阶模型 ，并在状态空间耦合 叶栅结构动力 学方程建立叶栅气动弹性方程 ，采 用特征值和时域仿真分析该系统稳定性。 运用该降阶耦合方法对 ＳＴＣＦ４（Ｓｔａｎｄａｒｄ ＴｅｓｔＣｏｎｆｉｇｕｒａｔｉｏｎ４）以 及 ＮＡＳＡ Ｒｏｔｏｒ６７叶栅系统的稳定性进行了 计算。 通过与直接计算流体力 学／计算结构动力学 （ＣＦＤ／ＣＳＤ）耦合方法和非耦合方法计算结果的比较验证了该方法的准确性 ，且该降价耦合方法的计算效率相对于直接 ＣＦＤ／ＣＳＤ 耦合方法提高了 １～２个 量级 ，为叶轮机气动弹性参数研究 、失谐研究以及多模态耦合计算等提供了便利 。 关键词：叶栅；气动弹性；颤振；降阶模型；计算流体力学；耦合方法 中图分类号 ：Ｖ２１１．４７ 文献标识码 ：Ａ 文章编号 ：１０００－６８９３（２０１４）１２－３２３２－１２ 现代叶轮机械的设计性 能 不断提 高 ，叶 轮 机叶片表 面 的 工 作 载 荷 也 随 之 增 加 ，而 叶 片 却 越来越轻 、越 薄 ，发 生 颤 振 的 可 能 性 日 渐 增 大。 叶轮机叶 片 颤 振 是 一 种 典 型 的 气 动 弹 性 现 象， 它能使叶 片 在 短 时 间 内 发 生 破 坏 ，危 害 性 非 常 大 ，避 免叶 片 颤 振已 成为现代 高性能叶轮机设 计中的一个重要 环 节。 研究 叶片颤振特性的方 法通常有经验法［１］和 数值方法［２－３］。 经验法在 大 量实验数据及工 程 经 验 的基 础上依靠经验公式 及简化模 型 来 进 行 叶 栅 颤 振 的 预 测 ，在 现 代 高 性能叶轮机设计环 节中 存 在一定的局限性 。 随 着计算机 的 发 展 和 数 值 模 拟 技 术 的 提 高 ，数 值 方法［４］在 叶轮机 气 动弹 性分 析领域发挥着越来 越重要的作用。 叶轮机颤振分析的数值方法可分为非耦合方 法和耦合方法［５］。非耦合方法不考虑结构和 流体 之间的互相作用，认为叶片的固有频率和模 态振 型不受气动力的影响 ，叶 片质量比较大时能较为 准确地预测系统的颤振特性。非耦合方法计算量 相对较小，近年 来仍被广泛用于叶轮机气动弹性 分析。能量法［６－７］是最为常见的非耦合方 法 ，它 通 常结合行波模型［８］用 于叶轮机单自由度 （或 单模 态）颤振分析。 由于未来叶轮机械设计指标的进一 步 提 升， 轻质材料的使用以及压比的进一步提高 ，使 得质 量比减小，流体和结构的耦合效应增强 ，传统非耦 合方法的合理性有待进一步验证。Ｂｅｎｄｉｋｓｅｎ 和 Ｆｒｉｅｄｍａｎｎ［９］认 为 结 构／流体之间的耦合影响叶 收稿日期 ：２０１４－０１－０８；退 修日期 ：２０１４－０２－２８；录 用日期 ：２０１４－０４－２４；网 络出版时间 ：２０１４－０６－０３１０：５７ 网络出版地址 ：ｗｗｗ．ｃｎｋｉ．ｎｅｔ／ｋｃｍｓ／ｄｅｔａｉｌ／１１．１９２９．Ｖ１０５７．００１．ｈｔｍｌ 基金项目：国家自然科学基金 ＊通讯作者 ．Ｔｅｌ．：０２９ Ｅ－ｍａｉｌ：ａｅｒｏｅｌａｓｔｉｃ＠ｎｗｐｕ．ｅｄｕ．ｃｎ 引用格式 ：ＳｕＤ，Ｚｈａｎｇ Ｗ Ｗ，ＱｕａｎＪＬ，ｅｔａｌ．ＡｎｅｆｆｉｃｉｅｎｔｃｏｕｐｌｅｄｍｅｔｈｏｄｏｆｃａｓｃａｄｅｆｌｕｔｔｅｒａｎａｌｙｓｉｓｂｙｕｓｉｎｇＣＦＤｔｅｃｈｎｉｑｕｅ［Ｊ］．Ａｃｔａ ＡｅｒｏｎａｕｔｉｃａｅｔＡｓｔｒｏｎａｕｔｉｃａＳｉｎ