风力机在复杂大气环境中使用会经受超出设计值的动态载荷,产生动态失速问题,导致风力机性能及使用寿命下降。为了研究风力机翼型动态失速特性,选取不同风速条件下的6种工况(风速范围为15m/s-40m/s),选用CFD方法,利用S-A湍流模型,对绕翼型1/4弦长点作正弦振荡中的翼型进行模拟计算。获取不同风速下的升力系数迟滞闭环规律,给出了一个完整振荡周期中不同时刻的流场分布,描述脱落涡从翼型表面发生、扩散、破碎的运动过程。通过研究动态失速升力系数迟滞闭环规律,发现其升力系数与风速呈线性增长关系,迟滞闭环的面积在风速小于35m/s下也随之增大。所研究的结果对风力机设计制造和维护具有指导意义。风力机动态特性-电动数控滚圆机滚弧机张家港钢管滚圆机滚弧机折弯机 荡运动的互相作用。2.2叶片网格划分本文计算S809翼型的振荡运动选用了动网格方法[12]，图1分别为翼型计算域网格和放大后翼型附近网格。风力机动态特性-电动数控滚圆机滚弧机张家港钢管滚圆机滚弧机折弯机整个计算域均采用三角形网格，共计约26400个网格单元，计算域分为内部动区域和外部静止区域。从外部静止区域到内部交界处网格相近以此减小插值误差。翼型处的网格密度较大，可在图1(b)放大图呈现。翼型计算域在进口边界设置来流速度，出口边界设置为静压力，在翼型表面边界满足无滑移壁面条件。(a)翼型计算域(b)翼型附近网格图1S809翼型网格模拟计算设置时间步长为0.001秒，且需要检查时间步长的无关性，为了得到周期性稳定重复的流场，每个工况经过6个以上连续周期的计算，并提取最稳定的数据进行后处理以及与OSU实验数据对比。3计算结果与分析3.1算例验证图2是S809翼型在风速32m/s条件下静态失速和动态失速升力系数迟滞图算例验证，相关计算参数如下：雷诺数Re=106，折合频率k=0.050，振荡频率f=1.2Hz，初始攻角α0=8°，振荡幅角α1=5.5°。如图2(a)所示，静态失速中计算值在小攻角下与实验值从图中可以看出升力系数计算值与OSU实验值在8°攻角前十分接近，在攻角超过8°以后计算值略大于实验值，实验值趋于平稳，说明此时出现流动分离。由图2(b)可知，利用S-A湍流模型模拟计算的气动力迟滞闭环曲线与OSU实验值具有良好的一致性，两个迟滞闭环曲线均为“8”字形，CFD计算值在刚开始的上仰过程中升力系数值与OSU实验值相持平，当攻角达到10°以后，计算值略大于实验值并达到峰值，因为在翼型表面出现漩涡，并且漩涡会附着在翼型表面上。在下俯阶段，大攻角下的计算值大于实验值，此时存在流动分离现象，湍流模型对于风力机动态特性-电动数控滚圆机滚弧机张家港钢管滚圆机滚弧机折弯机风力机动态特性-电动数控滚圆机滚弧机张家港钢管滚圆机滚弧机折弯机