风机叶片大厚度翼型气动弹性中的分叉现象.pdfVIP

第七届全国风能应用技术年会暨青海省等高原地区风能开发利用研讨会论文集 风机叶片大厚度翼型气动弹性中的分又现象+ 叶正寅蒋跃文武洁 (两北工业大学，翼型、叶栅窄气动力学国防科技霞点实验室，西安，邮编710072) 合翼型弹性振动方程，在时间域内模拟了不同厚度对称翼型的气动弹性过程，重点研究J，大迎角下的人厚度翼型流 场和气动弹性的分叉现象。研究结果发现：在酉万雷诺数条件下，厚度大于20％的对称翼型在一定大迎角范围内， 会f}{现气动弹性振动过程的分叉现象，通过跟踪流场的变化，发现引起翼型振动分叉的原凶在于翼型分离涡平衡态 的转化造成的。 关键词：气动弹性，分离流，颤振，分叉 引 言 当前，风能作为最环保的能源之一，风力发电技术受到了国内外的极人关注[1-3]o为了提高风 力发电的效率，目前大尺寸风机越来越受欢迎(如1给出了风机尺寸的发展变化情况，而A380人型 客机的翼展也只有79．5米)，这对风机叶片强度和刚度提出了新的挑战，采用大厚度翼型是提高毪 机叶片强度和刚度的根本途径，目前，风机翼犁厚度可以达到30％以上。从空气动力学角度讲，大 厚度翼型流动分离的可能性人幅度增加，另外，对于风机而言，为了提高风能效率，风机翼型一般 在较大迎角的T况下运动，因此，风机翼型一般会面临人迎角流动分离的问题。 0AOMWO．751蝴N 1．5MW 2．5MW 3．5MW 7．0MW 1985 1995 2000 2003 2006 2009 图l风力机发展的尺寸变化情况 作为一种转动机械，风机叶片的振动严重影响着叶片和整个风机系统的使用寿命及成本。叶片 的振动问题是空气动力学与结构的耦合力学问题即：气动弹性问题。传统气动弹性分析方法中使用 的nt‘素理论模型f4-51不能精确模拟流动分离问题，而大迎角分离流条件下的气动弹性特性对于风机尤 其重要【6刮。近年来，运用CFD计算气动弹性的丁作成为重要研究途径争川。本文在线性结构假设下， 性和流场特性进行研究，旨在探索大厚度翼型在大迎角条件下非线性特性。 ’本文得到国家自然科学基金资助 第七届全国风能应用技术年会暨青海省等高原地区风能开发利用研讨会论文集 1气动力计算方法 出于计算量和叶片大展弦比特征的考虑，仍然以二维翼型问题为研究对象，为了保证流动的粘 性特征得到很好模拟，同时考虑到非定常运动过程中重新生成网格的效率问题，采用改进的无限插 值理论n21生成绕翼型的C型网格，该方法可以很灵活地控制翼型表面附近的正交性和网格垂直方向 网格尺度变化。 流动基本方程为雷诺平均的可压缩、非定常N--S方程，其积分形式为： 昙妒协+∥矛如o (1) 其中： Q=[P，∥，pv，P】1 P，“，’，，e分别为空气密度、x，Y方向的速度分量和单位体积的总内能，历为线积分的法向单位向量， V为积分域，aV为包含体积分域的边界，F为通量项，它包括无粘项吒和粘性项E两部分： F=FE+Fv (2) 为了在非定常情况下，采用全隐式双时间方法求解非定常N—S方程，伪时间采用LU—SGS隐式 方法以提高计算效率，离散后的N—s方程在时间和空间上均具有二阶精度。此外，为了更好地计算 低速流动，对空气动力学基本方程运用了预处理方法…1。 图2翼型升力系数随迎角变化规律 图3翼型阻力系数随迎角变化 由于目前湍流转捩位置没有较好的确定方法，以流场为全湍流情况进行计算，湍流模型采用目 前认为较适合分离流动模拟的基丁k—CO的SST湍流模型。为检验数值方法和程序的适应性，以 态进行对比，来流马赫数为0．29，雷诺数为1950000，其迎角变化的数学表达式为： 口=口，+gosin(co·，)，减缩频率k=功·∥2·圪·口。表示平均迎角，％表示俯仰运动的迎角幅值。 缈俯仰运动的圆频率。算例中，平均迎角为17度，俯仰运动迎角的幅值为4．2度，减缩频率为O．1。 图2和图3分别为升力系数和阻力系数与迎角的关系，