基于beddoes-leshman模型的动态失速特性分析.docxVIP

基于beddoes-leshman模型的动态失速特性分析 0 动态失速模型 准确的动态动力负荷分析是风力机设计和修订的基础。对动态失速等非定常气动特性的精确建模,是风力机动态气动载荷分析的前提。动态失速是指当翼型经历周期性的或非定常的运动时,发生失速延迟的复杂现象。发生动态失速时的攻角偏离静态失速攻角,气动力系数可以极大地偏离静态值。工作过程中的风力机,叶片翼型的攻角是时变的。进行风力机运行工况仿真时,如果采用静态气动数据进行计算,将使得载荷被极大地低估。因此有必要研究动态失速模型以建立更精确的动态气动载荷分析程序。 国际上对动态失速模型的研究始于20世纪80年代。1981年,Tran和Petot建立了ONERA模型,使用微分方程描述气动力系数。1986年开始,Beddoes和Leishman对动态失速现象进行研究并建立了模型,该模型结合对流动分离的发生和影响的近似描述表现附着流的流动延迟效应。ye于1991年提出的动态失速模型忽略了附着流的瞬态效应,通过对用简单插值关系得到的流动的静态等效附着程度引入一个一阶滤波器来描述动态失速。2004年,丹麦Risφ国家实验室发展了Beddoes-Leishman模型的简化版本,忽略了流动的可压缩效应和开始于前缘的流动分离,引入了类似于ye模型的插值关系,使得模型适用于全范围攻角。 近年来,计算流体力学的发展使得确定运动翼型的气动响应已经不再困难,但求解过程仍然非常复杂。因此,在风力机气动设计领域中,利用半经验模型来求解非定常气动载荷仍将占主要地位。本文基于Beddoes-Leishman模型实现了水平轴风力机的翼型动态失速特性的数值计算方法,并将之与已有的风力机稳态性能分析程序集成,实现了水平轴风力机的动态气动载荷分析软件。 1 气动模型的建立 Beddoes-Leishman模型采用单位阶跃响应模型描述附着流,然后针对流动分离和涡流升力的影响进行了相应的扩展。模型采用固定在翼型上的坐标系,气动力系数分量是法向力系数CN(垂直于弦长方向)和切向力系数CT(弦长方向)。模型基于定常的气动数据CL和CD,因此首先要将其变换到翼型坐标系中。 下面分别从附着流和分离流两方面来描述Beddoes-Leishman模型。 1.1 非绕流升力模型 对于附着流,主要考虑尾流分离效应和气流加速效应这两方面因素对升力的影响,并分别建模。尾流分离效应对升力的影响部分称为绕流升力(circulatory lift),用CN,c表示;气流加速效应对升力的影响部分称为非绕流升力(noncirculatory lift),用CN,I表示。 绕流升力随攻角变化的方程为: 式中,CNα——法向力系数CN曲线的斜率;M——马赫数:φC——反映尾流分离效应的指数函数,近似为: 式中,s——无量纲时间,s=2Vrel·t/c,其中Vrel为风速,c为弦长;A1、A2、b1和b2——经验常数; 非绕流升力可以由活塞理论进行计算: 式中传递函数φ1可以被近似为: 其中,Kα——马赫数的函数,可表示为: 于是总附着流升力可以叠加绕流升力和非绕流升力得到: 1.2 动态动态升力模型 分离流部分主要考虑两方面的影响:①翼型后缘分离点的运动的动态延迟所产生的影响;②涡流的产生和分离所产生的影响。 首先考虑分离点运动的动态延迟。翼型表面的流动分离使绕流产生损失,从而使得由附着流计算得到的气动力系数减小。分离点在翼型表面的位置可以用f表示(f=x/c,x是分离点到后缘的距离,c是弦长)。分离点位置和气动力系数的关系可由基尔霍夫(Kirchhoff)流动理论得到: 式中,α0——CN=0时的攻角。 因此,静态情况下不同攻角的流动分离点可以利用定常气动力系数(式(8))求解。为获得动态流动分离点,采用一个一阶延迟方程来反映分离点相对于静态值的滞后: 式中,Tf——s空间的半经验时间常数。 为得到动态分离点的位置f,需要得到分离时的有效攻角,因此需要考虑由于后缘压力响应滞后导致的攻角偏移。采用一阶延迟方程描述压力响应滞后对附着流升力的影响: 式中,Tp——经验时间常数。于是有效攻角可以表示为: 得到有效攻角αf之后,可以求得静态分离点位置fstat,进而得到动态分离点位置f,由式(8)得到动态升力系数CN,f。 模型使用“涡流升力”来描述涡流的产生和分离对升力的影响。涡流升力的增量与非定常绕流升力CN,C和动态升力CN,f的差成比例,即: 同时,累积的涡流升力随时间按指数规律递减,但也有可能被新的涡流补充,因此涡流升力的变化规律可以表示为: 式中,Tv——经验时间常数。 于是,总升力系数可以表示为: 得到翼型坐标系下的非定常气动力系数CN和CT后,将其投影到垂直和平行于风速的方向,加上最小定常阻力系数CD0,就可以得到非定常升阻力系数CL和CD