CN102680226A - 一种水平轴风力机专用翼型性能评估方法 - Google Patents

一种水平轴风力机专用翼型性能评估方法 Download PDF

Info

Publication number
CN102680226A
CN102680226A CN2012101933806A CN201210193380A CN102680226A CN 102680226 A CN102680226 A CN 102680226A CN 2012101933806 A CN2012101933806 A CN 2012101933806A CN 201210193380 A CN201210193380 A CN 201210193380A CN 102680226 A CN102680226 A CN 102680226A
Authority
CN
China
Prior art keywords
aerofoil profile
parameter
performance
design
aeroperformance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012101933806A
Other languages
English (en)
Other versions
CN102680226B (zh
Inventor
杨科
李星星
白井艳
徐建中
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongke Guofeng science and Technology Co Ltd
Original Assignee
Institute of Engineering Thermophysics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Engineering Thermophysics of CAS filed Critical Institute of Engineering Thermophysics of CAS
Priority to CN201210193380.6A priority Critical patent/CN102680226B/zh
Publication of CN102680226A publication Critical patent/CN102680226A/zh
Application granted granted Critical
Publication of CN102680226B publication Critical patent/CN102680226B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种水平轴风力机专用翼型性能评估方法,首先获取翼型及其气动数据;分析翼型的几何结构参数,根据相对厚度的不同划分翼型类型为薄翼型和厚翼型系列;然后获得翼型的设计点性能参数项、非设计点性能参数项、失速特性参数项和稳定性参数项;继而与翼型性能的参考指标相比较,得出相对值,判断翼型各项特性的优良;最后根据翼型类型不同赋予翼型不同的权重系数,结合相对值大小,计算出翼型的综合性能评估值,判断风力机翼型的性能优良水平。本发明可以为风力叶片设计时不同部位选择不同翼型时提供判断的依据和方法,也可用于在新翼型的设计和优化中评估中间翼型的性能。

Description

一种水平轴风力机专用翼型性能评估方法
技术领域
本发明属于水平轴风力机翼型性能评估领域,具体涉及一种水平轴风力机叶片不同部位的翼型的性能评估方法。
背景技术
在全球性的能源危机和环境恶化这一背景下,风电产业迅猛发展。风力发电机本身也朝着运行环境多样化、单机容量变大的方向发展。从近海到沙漠,从平原到丘陵都有已经建成运行的风电场,新建成的风电场中风力机的叶片日趋大型化。这些趋势对风力机捕获风能的能力提出了更高要求。风力机获取风能的核心部件是叶片,叶片性能的优劣直接影响风力机风能转换的效率高低。从而叶片设计技术是风力机发电系统设计的核心技术之一。叶片可以视为一系列翼型沿着叶片展向叠加而成,叶片的性能由构成其气动外形的翼型的气动性能决定。叶片设计时,其不同部位由于结构需求和气动性能需求的差异,需要选择不同类型、不同性能的翼型。准确、完整地评估待选翼型的各项性能的优劣和综合性能的差异可以使得叶片设计时翼型选择的更加准确、合理。另一方面,设计出适合应用于在不同环境下保持高的风能转换率、低气动载荷、宽广的高性能工作范围且性能稳定的长叶片的新一代翼型系列对提高叶片性能具有基础性的作用。而在翼型的优化和设计中,准确、完整地评估中间翼型的气动性能,判断其是否达到设计目标是不可缺少的重要环节。可知,无论是进行叶片设计翼型的选择还是在翼型优化再设计时中间翼型的取舍,对翼型的性能进行完整、准确的评估都具有重要的意义。
目前,国内外风电行业上尚没有对风力机专用翼型性能的评估形成统一标准和完整的参数体系。现有的评估方法一般只能对翼型一部分气动参数进行定量比较:最大升阻比、最大升力系数、设计升力系数、阻力系数和力矩系数等;而对翼型的非设计点性能、失速特性以及在表面粗糙情况变化和流动雷诺数变化条件下性能的稳定性通常只根据翼型基本气动参数的曲线进行简单定性的比较分析,没有给出能够直观表征这些翼型性能的参数,以进行定量的比较分析。也就说,目前的风力机翼型性能评估体系中既没有完整的表征翼型各方面性能特征的参数体系,也没有合理且完整地评估翼型性能的方法;而这些恰是准确而完整地评估翼型性能的关键。
发明内容
鉴于以上现有技术的缺点和不足,本发明提出了一种水平轴风力机专用翼型性能评估方法,可以对风力机专用翼型的气动性能进行有效、合理且完整的评估。
本发明为实现上述技术目的所采取的技术方案为:
一种水平轴风力机专用翼型性能评估方法,其特征在于,包括以下步骤:
(1)获取翼型的几何坐标数据及其基本气动数据;
(2)利用翼型的几何坐标数分析得到翼型的几何参数项,根据所述几何参数项中的最大相对厚度判断翼型类型,即判断其属于薄翼型还是厚翼型;优选地,所述几何参数项包括最大相对厚度、最大相对厚度的弦向位置、最大相对弯度、最大相对弯度的弦向位置、前缘半径、尾缘厚度、尾缘角、上表面最大厚度值以及其弦向位置和下表面最大厚度值以及其弦向位置。
(3)分析翼型的基本气动数据,得到翼型的各气动性能参数项;
(4)根据翼型类型选择所述各气动性能参数项的相应参考值,并计算得出所述各气动性能参数项相对于其参考值的相对值,据此判断翼型各项气动特性的优劣;
(5)根据叶片不同部位对翼型性能要求的不同赋予翼型各气动性能参数项以相应的权重系数,结合步骤(4)中得到的翼型各气动性能参数项的相对值,计算出翼型的综合性能评估值,据此判断翼型的性能优良水平。
优选地,所述的步骤(1)中翼型的基本气动数据是风洞试验结果或者数值模拟计算结果,包括自然转捩和固定转捩两种情况下(研究翼型性能受表面粗糙状况不同的影响)不同流动雷诺数时的各项气动数据,所述各项气动数据主要包括不同流动攻角下的翼型的升力系数、阻力系数、力矩系数、上表面分离点位置、下表面分离点位置等数据。风洞试验中,自然转捩对应试验段叶片表面光滑无污染的情况,固定转捩对应试验段叶片靠近前缘某一处加粗糙带的情况(如在上表面5%弦长处贴zz型粗糙带,粗糙带的具体尺寸取决于叶片试验段的尺寸);数值模拟计算时,自然转捩意味着不限制流动转捩和分离点位置,让其随着流动规律自然发生;固定转捩是通过限制翼型上下表面转捩点的位置(如通常设置翼型固定转捩的位置为:上表面1%弦长处,下表面10%弦长处)来达到实现的,以模拟表面粗糙程度不同所造成影响。
优选地,所述的步骤(2)中,利用翼型的几何坐标数据可以通过多种翼型分析软件得到翼型的几何参数项;根据翼型的最大相对厚度可以确定翼型的类型和可适用于叶片上的位置。所述薄翼型是指最大相对厚度小于25%的翼型,所述厚翼型是指相对厚度大于或等于25%的翼型。
进一步地,所述步骤(3)中,所述翼型的各气动性能参数项包括翼型的设计点性能参数项、非设计点性能参数项、失速特性参数项和稳定性参数项。
所述翼型的设计点性能参数项,其所包括的各参数来自于风洞试验结果或者数值计算结果,可以给出翼型在流动中获得的升力、阻力以及气动效率;翼型的非设计点性能参数项,将翼型最大升阻比所对应攻角和最大升力系数攻角之间翼型的气动性能参数化,来表征翼型从设计点到失速点间的性能特征;翼型的失速特性参数项,根据翼型的升力特性曲线,以最大升力系数为基础,用失速后某一攻角处的升力系数平均变化率和升力系数变化量的乘积来表征翼型的失速后升力曲线的平缓特性;翼型的气动性能稳定性参数项:包括翼型的气动性能随转捩方式不同的稳定性(此为表征翼型的表面粗糙敏感性)和随流动雷诺数Re的不同的稳定性。
进一步地,所述设计点性能参数项包括:最大升阻比、设计攻角、设计升力系数、设计阻力系数、设计力矩系数;所述非设计点性能参数项包括:最大升力系数和最大升力系数攻角、设计攻角到失速攻角的间隔Δα、设计攻角到失速攻角间升力系数的变化率和升阻比的变化率;所述失速特性参数项包括翼型失速特性参数Mstall;所述稳定性参数项包括:翼型性能对表面粗糙影响的稳定性参数psf,对流动雷诺数变化影响的稳定性参数psRe
进一步地,设计攻角到失速攻角的间隔Δα、设计攻角到失速攻角间升力系数的变化率CRcl和升阻比的变化率CRld可分别表示为:
Δα=αstalldesign CRcl = cl max - cl design Δα , CRld = | l / d stall - l / d max | Δα ,
其中αstall(失速攻角)表示翼型升力系数初次达到最大时的攻角,clmax即为升力系数首次达到最大的值,l/dstall即为相应的升阻比;αdesign表示设计攻角,即升阻比达到最大时的攻角,cldesign即为最大升阻比对应攻角下的升力系数(设计升力系数),l/dmax即为最大升阻比。
上述翼型失速特性参数Mstall为翼型失速后十度攻角范围内某一攻角处的升力系数相对于最大升力系数处的平均变化率和升力系数变化量的乘积的最大值:
M stall = max { ( cl - cl max ) 2 &alpha; - &alpha; stall , 0 < &alpha; - &alpha; stall &le; 10 and&alpha; &Element; Z }
其中α,cl分别为失速后某一攻角和该攻角下的升力系数;clmax表示最大升力系数。
翼型性能对表面粗糙影响的稳定性通过稳定性参数psf来量化:
p sf = | p nt - p ft | p nt &times; 100 %
其中p和psf分别表示翼型的某一气动性能参数项(翼型设计点性能参数项、非设计点性能参数项或失速特性参数项)及其对表面粗糙影响的稳定性参数。pnt为翼型所在光滑叶片表面的实验值或自然转捩的数值计算结果,pft为翼型所在粗糙叶片表面的实验值或设置固定转捩的数值计算结果。
翼型性能对流动雷诺数变化时的稳定性通过稳定性参数psRe来量化:
p sRe = { max ( | q Re - q Redesign | q Redesign &times; 100 % ) , Re = 1 2 Re design orRe = 2 Re design }
该式表示当雷诺数分别为设计雷诺数的一半和二倍时,分别计算翼型性能相对变化率,翼型性能对流动雷诺数变化影响的稳定性参数psRe取二者中的较大值;其中qRe为翼型在一定雷诺数下某性能参数,qRe design为对应设计雷诺数下的性能参数。
所述的步骤(4)中,评估翼型所需的所述各气动性能参数项的相应参考值来自:如果对翼型的性能评估适用于叶片设计选择翼型的过程中,则从已有的风力机专用翼型数据库中,选择相同相对厚度或者相同翼型类型的典型翼型的气动性能参数作为参考值;或,如果对翼型的性能评估适用于翼型的优化和设计过程中,则以给定的目标翼型的气动性能参数作为参考值,该目标翼型的气动性能参数的给定取决于待优化设计翼型的气动性能和拟设计的方向。进一步地,得到翼型的所述各气动性能参数项的相应参考值后,将待评估翼型的各气动性能参数项与相应参考值比较得到各气动性能参数项的相对值,形式如下:
Figure BDA00001757622300061
其中p代表翼型的某一气动性能参数项,pr代表相应的参考值,p0代表相对值。p0的大小反应了待评估翼型对参考翼型气动性能的偏离,从而可以反应翼型各项气动性能的优良水平。翼型的一部分性能参数当p0越大该项性能越优良,另一部分性能参数则反之。
所述的步骤(5)中不同类型的翼型所匹配的权重系数是依据叶片不同部位对翼型气动性能要求的差异来决定的,薄翼型和厚翼型有不同的权重系数。将翼型的各参数项的权重系数Ai设定为第一级权重系数,满足
Figure BDA00001757622300062
即各参数项的权重系数Ai之和为1;将各参数项中所包括的各参数的权重系数Aij设定为第二级权重系数,满足
Figure BDA00001757622300063
即各参数项中所包括的各参数的权重系数Aij之和为1。
然后结合已经得到的待评估翼型的各气动性能参数项的相对值,可以计算翼型的综合性能评估参数P:
Figure BDA00001757622300064
其中p0ij为各参数项中某一性能参数的相对值,与权重系数相对应;qij与该某一性能参数的期望方向有关,如果参数越大性能越优良取值为1,反之取值为-1。
进一步地,所述的步骤(5)中,对风力机专用翼型性能优良的评估还包括:
a)将翼型性能优良程度分为“优异”、“良好”、“合格”、“缺陷”、“恶劣”等五个等级;
b)通过各气动性能参数项的相对值p0的大小,判断翼型各项气动性能参数与参考值的偏离,根据翼型各项气动性能的特征设定各自的允许范围,当翼型具有一项或两项气动性能参数项的相对值超出所述允许范围时,判定翼型为“缺陷”,即其性能某方面存在缺陷;当翼型有至少三项气动性能参数项超出所述允许范围时,判定翼型为“恶劣”,即翼型性能存在严重缺陷;如果翼型各气动性能参数项的相对值均在所述允许范围以内时,则初步判定翼型为“合格”;
c)当初步判定翼型为“合格”后,计算综合性能评估参数P,得出P值;根据P值大小划分翼型属于“合格”、“良好”或“优异”。
附图说明
图1是风力机专用翼型性能参数体系图;
图2是风力机专用翼型性能评估流程图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明进一步详细说明。
本发明包含两个相联系的方面:风力机专用翼型性能参数体系和性能评估评估方法。下面结合附图做出进一步的说明。
如图1所示,水平轴风力机专用翼型性能评价参数体系包括:翼型几何参数项,翼型的设计点性能参数项,翼型的非设计点性能参数项,翼型的失速特性参数项,翼型的气动性能稳定性参数项,翼型的比较权重系数项。
其中,翼型几何参数项包括最大相对厚度、最大相对厚度的弦向位置、最大相对弯度、最大相对弯度的弦向位置、前缘半径、尾缘厚度、尾缘角、上表面最大厚度值以及其弦向位置和下表面最大厚度值以及其弦向位置。
具体地,翼型几何参数项中各几何参数的符号和定义如下表一所示。
所述设计点性能参数项包括:最大升阻比、设计攻角、设计升力系数、设计阻力系数、设计力矩系数。所述非设计点性能参数项包括:最大升力系数和最大升力系数攻角、设计攻角到失速攻角的间隔Δα、设计攻角到失速攻角间升力系数的变化率和升阻比的变化率。所述失速特性参数项包括翼型失速特性参数Mstall。所述稳定性参数项包括:翼型性能对表面粗糙影响的稳定性参数psf,对流动雷诺数变化影响的稳定性参数psRe
具体地,翼型的设计点性能参数项、非设计点性能参数项和失速特性参数项所包含的各项参数的符号及其定义如下表二所示。
表一翼型的几何参数项
  翼型的几何参数表   符号   定义
  相对厚度   T   翼型的最大厚度与弦长的比值
  最大相对厚度位置   Xt   翼型最大相对厚度处所在弦向的相对位置
  相对弯度   Cm   翼型的最大弯度与弦长的比值
  最大相对弯度位置   Xcm   翼型最大弯度处所在的弦向的相对位置
  前缘半径   Rle   翼型前缘内切圆的半径大小与弦长比值
  尾缘厚度   Ttr   翼型尾缘处上下表面间距与弦长的比值
  尾缘角   θ   翼型尾缘处上下表面切线的夹角大小
  上表面最大厚度   Yup   翼型上表面点离弦线最远距离与弦长的比值
  上表面最大厚度位置   Xup   翼型上表面最大厚度处所在弦向相对位置
  下表面最大厚度   Ybel   翼型下表面点离弦线最远距离与弦长的比值
  下表面对大厚度位置   Xbel   翼型下表面最大厚度处所在弦向相对位置
表二翼型的设计点性能参数项、非设计点性能参数项和失速特性参数项
Figure BDA00001757622300081
稳定性参数项包括:翼型性能对表面粗糙影响的稳定性参数psf,对流动雷诺数变化影响的稳定性参数psRe。翼型性能对表面粗糙影响的稳定性的量化是通过翼型性能参数受表面粗糙影响的相对变化率来给定参数psf;翼型性能对流动雷诺数变化影响的稳定性的量化是通过当雷诺数为设计雷诺数的一半或二倍时,翼型性能相对变化率的较大值给定参数psRe
如图2所示,本发明的风力机专用翼型族的性能评估方法包括以下步骤:(1)获取翼型的几何坐标数据及其基本气动数据;(2)利用翼型的几何坐标数分析得到翼型的几何参数项,根据所述几何参数项中的最大相对厚度判断翼型类型,即判断其属于薄翼型还是厚翼型;(3)分析翼型的基本气动数据,得到翼型的各气动性能参数项;(4)根据翼型类型选择所述各气动性能参数项的相应参考值,并计算得出所述各气动性能参数项相对于其参考值的相对值,据此判断翼型各项气动特性的优劣;(5)根据叶片不同部位对翼型性能要求的不同赋予翼型各气动性能参数项以相应的权重系数,结合步骤(4)中得到的翼型各气动性能参数项的相对值,计算出翼型的综合性能评估值,据此判断翼型的性能优良水平。
上述步骤(4)中,评估翼型所需的气动性能参数的参考值来自下述其中一个方面:如果对翼型的评估适用于叶片设计选择翼型的过程中,则选择相同相对厚度或者应用于叶片相同位置已经设计成型的风力机专用翼型系列中,性能优良并取得广泛应用的典型翼型(如荷兰的DU翼型系列)的气动参数为参考值;如果对翼型的评估适应于翼型的优化和设计过程中,则以给定的目标翼型的气动性能为参考值,该目标翼型气动数据的给定取决于待优化设计翼型的气动性能和拟设计的方向。得到完整的翼型性能参数的参考值后,将待评估翼型的各性能参数项与参考值比较得到各项性能的相对值p0,形式如下:
Figure BDA00001757622300091
其中p代表翼型的某一项气动性能参数,pr代表相应的参考值,p0代表相对值。p0的大小反应了待评估翼型对参考翼型气动性能的偏离,从而可以反应翼型各项气动性能的优良水平。翼型的一部分性能参数当p0越大该项性能越优良,另一部分性能参数则反之。
翼型各气动性能参数项相对值的符号表示和各项性能参数的期望方向、翼型性能参数相对值的可接受范围见表三中所示。需要说明的是,表格中期望方向一列,如果我们希望翼型某一项的性能参数比参考翼型的性能参数大,即p0>1,则记为“+”,反之则记为“-,“/”表示不作约束,或者说性能参数的变化有其他条件决定;“可接受范围”一列表示翼型性能参数的相对值的可以接受的变化范围。在稳定性参数中,如果出现参考值pr=0,当且仅当待评估参数值p也为0时,相对值定义为1,其他均定义为∞。
上述步骤(5)中,根据叶片不同位置对翼型性能要求的差异匹配不同厚度翼型的权重系数。根据叶片设计时不同部位采用的翼型的相对厚度范围,把翼型相对厚度在25%以下的称为风力机薄翼型系列(典型的相对厚度为:15%、18%、21%);大于该相对厚度的属于风力机厚翼型系列(典型的相对厚度为:25%、30%、35%、40%;45%以上为大厚度翼型)。在长叶片设计中薄翼型适用于叶片外侧,厚翼型适用于叶片中部和内侧。叶片外侧、中部以及内侧的翼型有各自的权重系数。不同类型翼型的第一级权重系数Ai分配如下表四所示。
不同类型的翼型各项性能参数项中,第二级权重系数Aij的分配如表五所示(稳定性参数的子项所属的第二级权重系数按照所含参数的数目平权分配,即此时的Aij=1/11)。
然后结合已经得到的待评估翼型的相对值,可以计算翼型的综合性能评估参数P:
Figure BDA00001757622300101
其中p0ij为某一性能参数的相对值,与权重系数相对应。qij与该项性能参数的期望方向有关,期望方向为“+”取值为1,期望方向为“-”取值为-1.
表三翼型各项性能参数相对值的期望方向和可接受范围
  相对值   期望方向   可接受范围   相对值   期望方向   可接受范围
  p0   +/-   ()   p0   +/-   ()
  l/dmax0   +   αstall,sf0   -   (0,1.25)
  αdesign0   /   /   Δαsf0   -   (0,1.25)
  cldesign0   +   CRclsf0   -   (0,1.25)
  cddesign0   -   (0,1.25)   CRldsf0   -   (0,1.25)
  cmdesign0   -   (0,1.25)   Mstallsf0   -   (0,1.25)
  clmax0   +   l/dmax,sRe0   -   (0,1.25)
  αstall0   /   (0,1.25)   αdesign,sRe0   -   (0,1.25)
  Δα0   +   cldesign,sRe0   -   (0,1.25)
  CRcl0   -   (0,1.25)   cddesign,sRe0   -   (0,1.25)
  CRld0   -   (0,1.25)   cmdesign,sRe0   -   (0,1.25)
  Mstall0   -   (0,1.25)   clmax,sRe0   -   (0,1.25)
  l/dmax,sf0   -   (0,1.25)   αstall,sRe0   -   (0,1.25)
  αdesign,sf0   -   (0,1.25)   ΔαsRe0   -   (0,1.25)
  cldesign,sf0   -   (0,1.25)   CRclsRe0   -   (0,1.25)
  cddesign,sf0   -   (0,1.25)   CRldsRe0   -   (0,1.25)
  cmdesign,sf0   -   (0,1.25)   Mstall,sRe0   -   (0,1.25)
  clmax,sf0   -   (0,1.25)
表四翼型权重系数Ai分配表
Figure BDA00001757622300111
表五第二级权重系数Aij分配表
Figure BDA00001757622300121
对风力机专用翼型性能优良的判断还包含以下具体内容:
a)将翼型性能优良程度分为“优异”、“良好”、“合格”、“缺陷”、“恶劣”等五个等级;
b)通过各气动性能参数项的相对值p0的大小,判断翼型各气动性能参数项与参考值的偏离,根据翼型各项气动性能的特征设定各自的允许范围,当翼型具有一项或两项气动性能参数项的相对值超出所述允许范围时,判定翼型为“缺陷”,即其性能某方面存在缺陷;当翼型有至少三项气动性能参数项超出所述允许范围时,判定翼型为“恶劣”,即翼型性能存在严重缺陷;如果翼型各气动性能参数项的相对值均在所述允许范围以内时,则初步判定翼型为“合格”;
c)当初步判定翼型为“合格”后,计算综合性能评估参数P,得出P值;根据P值大小划分翼型属于“合格”、“良好”或“优异”。
翼型性能的优劣评定等级如表六所示。
表六翼型性能的优劣评定等级
  等级   评价   性能状况说明
  1   优异   各项气动性能参数满足“+”:p0∈(1.2,+∞);“-”:p0∈(0,0.8)
  2   良好   各项气动性能参数满足“+”:p0∈(1,1.2);“-”:p0∈(0.8,1)
  3   合格   各项气动性能参数满足“+”:p0∈(0.8,1);“-”:p0∈(1,1.25)
  4   缺陷   有一到两项翼型气动性能参数超出允许范围
  5   恶劣   有三项以上翼型气动性能参数超出允许范围
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的范围之内。

Claims (10)

1.一种水平轴风力机专用翼型性能评估方法,其特征在于,所述评估方法包括以下步骤:
(1)获取翼型的几何坐标数据及其基本气动数据;
(2)利用翼型的几何坐标数分析得到翼型的几何参数项,根据所述几何参数项中的最大相对厚度判断翼型类型,即判断其属于薄翼型还是厚翼型;
(3)分析翼型的基本气动数据,得到翼型的各气动性能参数项,包括翼型设计点性能参数项、非设计点性能参数项、失速特性参数项和稳定性参数项;
(4)根据翼型类型选择所述各气动性能参数项的相应参考值,并计算得出所述各气动性能参数项相对于其参考值的相对值,据此判断翼型各项气动特性的优劣;
(5)根据叶片不同部位对翼型性能要求的不同赋予翼型各气动性能参数项以相应的权重系数,结合步骤(4)中得到的翼型各气动性能参数项的相对值,计算出翼型的综合性能评估值,据此判断翼型的性能优良水平。
2.根据权利要求1所述的评估方法,其特征在于:上述步骤(1)中,翼型的基本气动数据是风洞试验结果或者数值模拟计算结果,包括自然转捩和固定转捩两种情况下不同流动雷诺数时的各项气动数据,所述各项气动数据包括不同流动攻角下的翼型的升力系数、阻力系数、力矩系数、上表面分离点位置和下表面分离点位置。
3.根据权利要求1所述的评估方法,其特征在于:所述的步骤(2)中,所述薄翼型为最大相对厚度小于25%的翼型,所述厚翼型为最大相对厚度大于或等于25%的翼型。
4.根据权利要求1所述的评估方法,其特征在于:所述的步骤(4)中,评估翼型所需的所述各气动性能参数项的相应参考值来自:如果对翼型的性能评估适用于叶片设计选择翼型的过程中,则从已有的风力机专用翼型数据库中,选择相同相对厚度或者相同翼型类型的典型翼型的气动性能参数作为参考值;或,如果对翼型的性能评估适用于翼型的优化和设计过程中,则以给定的目标翼型的气动性能参数作为参考值,该目标翼型的气动性能参数的给定取决于待优化设计翼型的气动性能和拟设计的方向。
5.根据权利要求4所述的评估方法,其特征在于:所述步骤(4)中,得到翼型的所述各气动性能参数项的相应参考值后,将待评估翼型的各气动性能参数项与相应参考值比较得到各气动性能参数项的相对值p0,形式如下:
Figure FDA00001757622200021
其中,p为翼型的某一气动性能参数项,pr为相应的参考值。
6.根据权利要求5所述的评估方法,其特征在于:所述的步骤(5)中,不同类型的翼型所匹配的权重系数是依据叶片不同部位对翼型气动性能要求的差异来决定的,叶尖、叶中和叶根处的翼型有不同的权重系数;将翼型的各气动性能参数项的权重系数Ai设定为第一级权重系数,满足
Figure FDA00001757622200022
将各气动性能参数项中所包括的各参数的权重系数Aij设定为第二级权重系数,满足
Figure FDA00001757622200023
结合已经得到的待评估翼型的各气动性能参数项的相对值,可以计算翼型的综合性能评估参数P:
Figure FDA00001757622200024
其中p0ij为各参数项中某一性能参数的相对值,与权重系数相对应;qij与该某一性能参数的期望方向有关,如果参数越大性能越优良取值为1,反之取值为-1。
7.根据权利要求6所述的评估方法,其特征在于:所述的步骤(5)中,对风力机专用翼型性能优良的评估还包括:
a)将翼型性能优良程度分为“优异”、“良好”、“合格”、“缺陷”、“恶劣”等五个等级;
b)通过各气动性能参数项的相对值p0的大小,判断翼型各气动性能参数项与参考值的偏离,根据翼型各项气动性能的特征设定各自的允许范围,当翼型具有一项或两项气动性能参数项的相对值超出所述允许范围时,判定翼型为“缺陷”,即其性能某方面存在缺陷;当翼型有至少三项气动性能参数项超出所述允许范围时,判定翼型为“恶劣”,即翼型性能存在严重缺陷;如果翼型各气动性能参数项的相对值均在所述允许范围以内时,则初步判定翼型为“合格”;
c)当初步判定翼型为“合格”后,计算综合性能评估参数P,得出P值;根据P值大小划分翼型属于“合格”、“良好”或“优异”。
8.根据权利要求上述任一项所述的评估方法,其特征在于:所述设计点性能参数项包括:最大升阻比、设计攻角、设计升力系数、设计阻力系数、设计力矩系数;所述非设计点性能参数项包括:最大升力系数、最大升力系数攻角(失速攻角)、设计攻角到失速攻角的间隔Δα、设计攻角到失速攻角间升力系数的变化率和升阻比的变化率;所述失速特性参数项包括翼型失速特性参数Mstall;所述稳定性参数项包括:翼型性能对表面粗糙影响的稳定性参数psf,对流动雷诺数变化影响的稳定性参数psRe
9.根据上述任一项权利要求所述的评估方法,其特征在于:设计攻角到失速攻角的间隔Δα、设计攻角到失速攻角间升力系数的变化率CRcl和升阻比的变化率CRld可分别表示为:
Δα=αstalldesign CRcl = cl max - cl design &Delta;&alpha; , CRld = | l / d stall - l / d max | &Delta;&alpha; ,
其中αstall表示翼型升力系数初次达到最大时的攻角(失速攻角),clmax即为升力系数首次达到最大的值,l/dstall为该攻角下的升阻比;αdesign表示设计攻角,即升阻比达到最大时的攻角,cldesign为最大升阻比对应攻角下的升力系数(设计升力系数),l/dmax为最大升阻比;
上述翼型失速特性参数Mstall为翼型失速后十度攻角范围内某一攻角处的升力系数相对于最大升力系数处的平均变化率和升力系数变化量的乘积的最大值: M stall = max { ( cl - cl max ) 2 &alpha; - &alpha; stall , 0 < &alpha; - &alpha; stall &le; 10 and&alpha; &Element; Z } , 其中:α,cl分别为失速后某一攻角和该攻角下的升力系数;clmax表示最大升力系数;
翼型性能对表面粗糙影响的稳定性通过稳定性参数psf来量化:
p sf = | p nt - p ft | p nt &times; 100 %
其中p和psf分别表示翼型的某一气动性能参数和相应稳定性参数,pnt为翼型所在光滑叶片表面的实验值或自然转捩的数值计算结果,pft为翼型所在粗糙叶片表面的实验值或设置固定转捩的数值计算结果;
翼型性能对流动雷诺数变化的稳定性通过稳定性参数psRe来量化:
p sRe = { max ( | q Re - q Redesign | q Redesign &times; 100 % ) , Re = 1 2 Re design orRe = 2 Re design }
该式表示当雷诺数分别为设计雷诺数的一半和二倍时,分别计算翼型性能相对变化率,翼型性能对流动雷诺数变化影响的稳定性参数psRe取二者中的较大值;其中qRe为翼型在一定雷诺数下某性能参数,qRe design为对应设计雷诺数下的性能参数。
10.根据上述任一项权利要求所述的评估方法,其特征在于:翼型的几何参数项包括最大相对厚度、最大相对厚度的弦向位置、最大相对弯度、最大相对弯度的弦向位置、前缘半径、尾缘厚度、尾缘角、上表面最大厚度值以及其弦向位置和下表面最大厚度值以及其弦向位置。
CN201210193380.6A 2012-06-12 2012-06-12 一种水平轴风力机专用翼型性能参数测试方法 Active CN102680226B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210193380.6A CN102680226B (zh) 2012-06-12 2012-06-12 一种水平轴风力机专用翼型性能参数测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210193380.6A CN102680226B (zh) 2012-06-12 2012-06-12 一种水平轴风力机专用翼型性能参数测试方法

Publications (2)

Publication Number Publication Date
CN102680226A true CN102680226A (zh) 2012-09-19
CN102680226B CN102680226B (zh) 2014-06-04

Family

ID=46812491

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210193380.6A Active CN102680226B (zh) 2012-06-12 2012-06-12 一种水平轴风力机专用翼型性能参数测试方法

Country Status (1)

Country Link
CN (1) CN102680226B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104612892A (zh) * 2014-12-30 2015-05-13 中国科学院工程热物理研究所 一种风力机翼型的多学科优化设计方法
CN104990684A (zh) * 2015-06-26 2015-10-21 西安理工大学 一种基础翼型的水动力特性测试装置及测试方法
CN105115692A (zh) * 2015-08-12 2015-12-02 中国航天空气动力技术研究院 一种cfd数值模拟与风洞试验相结合的气动预测方法
CN105223006A (zh) * 2015-09-17 2016-01-06 中国运载火箭技术研究院 一种动力试车体襟翼抗冲刷试验的操控装置及试验方法
CN108303228A (zh) * 2017-08-08 2018-07-20 北京空天技术研究所 一种风洞试验的边界层流态控制装置及方法
CN109145506A (zh) * 2018-09-20 2019-01-04 中国科学院工程热物理研究所 一种高气动性能低噪声水平风力机外侧翼型的优化设计方法
CN111639404A (zh) * 2020-05-28 2020-09-08 中国科学院工程热物理研究所 带锯齿尾缘的仿生降噪翼型优化方法、装置及电子设备
CN112362291A (zh) * 2020-09-30 2021-02-12 成都飞机工业(集团)有限责任公司 一种飞翼布局飞机纵向气动力系数的雷诺数效应修正方法
WO2022011961A1 (zh) * 2020-11-05 2022-01-20 中国长江三峡集团有限公司 基于ptfe膜对风机叶片气动特性影响的数值模拟方法
CN114611333A (zh) * 2022-05-10 2022-06-10 中国航发上海商用航空发动机制造有限责任公司 压气机效率评估方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101532906A (zh) * 2009-04-27 2009-09-16 东南大学 风力发电机叶片的流体动力学和结构力学分析方法
US20100263448A1 (en) * 2007-12-14 2010-10-21 Alliance For Sustainable Energy, Llc. Dual- axis resonance testing of wind turbine blades
US20110010108A1 (en) * 2008-10-31 2011-01-13 General Electric Company System and method for monitoring health of airfoils
CN102004838A (zh) * 2010-12-17 2011-04-06 中国航天空气动力技术研究院 一种基于有限差分法的风力机叶片结构确定方法
CN102184294A (zh) * 2011-05-11 2011-09-14 北京动力机械研究所 小型涡轮发动机涡轮气动性能虚拟试验系统及试验方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100263448A1 (en) * 2007-12-14 2010-10-21 Alliance For Sustainable Energy, Llc. Dual- axis resonance testing of wind turbine blades
US20110010108A1 (en) * 2008-10-31 2011-01-13 General Electric Company System and method for monitoring health of airfoils
CN101532906A (zh) * 2009-04-27 2009-09-16 东南大学 风力发电机叶片的流体动力学和结构力学分析方法
CN102004838A (zh) * 2010-12-17 2011-04-06 中国航天空气动力技术研究院 一种基于有限差分法的风力机叶片结构确定方法
CN102184294A (zh) * 2011-05-11 2011-09-14 北京动力机械研究所 小型涡轮发动机涡轮气动性能虚拟试验系统及试验方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨科等: "《基于CFD技术的高性能风力机翼型最优化设计方法》", 《工程热物理学报》, vol. 28, no. 4, 31 July 2007 (2007-07-31), pages 586 - 588 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104612892A (zh) * 2014-12-30 2015-05-13 中国科学院工程热物理研究所 一种风力机翼型的多学科优化设计方法
CN104612892B (zh) * 2014-12-30 2017-03-29 中国科学院工程热物理研究所 一种风力机翼型的多学科优化设计方法
CN104990684A (zh) * 2015-06-26 2015-10-21 西安理工大学 一种基础翼型的水动力特性测试装置及测试方法
CN105115692A (zh) * 2015-08-12 2015-12-02 中国航天空气动力技术研究院 一种cfd数值模拟与风洞试验相结合的气动预测方法
CN105115692B (zh) * 2015-08-12 2017-06-09 中国航天空气动力技术研究院 一种cfd数值模拟与风洞试验相结合的气动预测方法
CN105223006A (zh) * 2015-09-17 2016-01-06 中国运载火箭技术研究院 一种动力试车体襟翼抗冲刷试验的操控装置及试验方法
CN105223006B (zh) * 2015-09-17 2017-08-29 中国运载火箭技术研究院 一种动力试车体襟翼抗冲刷试验的操控装置及试验方法
CN108303228A (zh) * 2017-08-08 2018-07-20 北京空天技术研究所 一种风洞试验的边界层流态控制装置及方法
CN109145506A (zh) * 2018-09-20 2019-01-04 中国科学院工程热物理研究所 一种高气动性能低噪声水平风力机外侧翼型的优化设计方法
CN111639404A (zh) * 2020-05-28 2020-09-08 中国科学院工程热物理研究所 带锯齿尾缘的仿生降噪翼型优化方法、装置及电子设备
CN111639404B (zh) * 2020-05-28 2024-01-23 中国科学院工程热物理研究所 带锯齿尾缘的仿生降噪翼型优化方法、装置及电子设备
CN112362291A (zh) * 2020-09-30 2021-02-12 成都飞机工业(集团)有限责任公司 一种飞翼布局飞机纵向气动力系数的雷诺数效应修正方法
CN112362291B (zh) * 2020-09-30 2021-12-10 成都飞机工业(集团)有限责任公司 一种飞翼布局飞机纵向气动力系数的雷诺数效应修正方法
WO2022011961A1 (zh) * 2020-11-05 2022-01-20 中国长江三峡集团有限公司 基于ptfe膜对风机叶片气动特性影响的数值模拟方法
GB2604759A (en) * 2020-11-05 2022-09-14 China Three Gorges Renewables Group Co Ltd Numerical simulation method based on influence of PTFE membrane on aerodynamic characteristics of fan blade
GB2604759B (en) * 2020-11-05 2023-03-01 China Three Gorges Renewables Group Co Ltd Numerical simulation method based on influence of PTFE membrane on aerodynamic characteristics of fan blade
CN114611333A (zh) * 2022-05-10 2022-06-10 中国航发上海商用航空发动机制造有限责任公司 压气机效率评估方法及系统
CN114611333B (zh) * 2022-05-10 2022-08-09 中国航发上海商用航空发动机制造有限责任公司 压气机效率评估方法及系统

Also Published As

Publication number Publication date
CN102680226B (zh) 2014-06-04

Similar Documents

Publication Publication Date Title
CN102680226B (zh) 一种水平轴风力机专用翼型性能参数测试方法
CN109190283A (zh) 一种考虑高湍流自由来流效应的风力机翼型气动优化方法
CN109460566B (zh) 一种风力机叶片内侧厚翼型的气动稳健优化设计方法
CN104778327B (zh) 基于人工神经网络的飞机翼型优化设计方法
CN104318482A (zh) 一套智能配电网综合评估体系和方法
CN110111024A (zh) 基于ahp模糊综合评价模型的科技成果市场价值评估方法
CN103617371A (zh) 一种基于灰色理论的电能质量综合评估方法
CN103839192A (zh) 基于层次分析法及距离综合评价法的风电场综合评价方法
CN104933627A (zh) 机床产品制造系统的能效组合评价方法
CN106126791B (zh) 一种考虑几何不确定性的高超声速机翼气动力/热分析方法
CN104933483A (zh) 一种基于天气过程划分的风电功率预测方法
CN109145506A (zh) 一种高气动性能低噪声水平风力机外侧翼型的优化设计方法
CN113408044B (zh) 一种多旋翼无人机桨叶优化设计方法
CN105005878B (zh) 一种坚强智能电网的综合评估方法
CN203374428U (zh) 一族大厚度钝尾缘风力机翼型
CN106126860A (zh) 一种考虑加工误差的高超声速机翼鲁棒优化设计方法
CN110288196A (zh) 输电线路舞动分级预警方法及系统
CN109409019A (zh) 一种风力机叶片气动外形的预弯优化方法
CN105046041A (zh) 一种压气机可控扩散叶型稳健优化设计方法
Castelli et al. Numerical analysis of the influence of airfoil asymmetry on VAWT performance
CN105678453A (zh) 一种基于电网指标评价体系的多维度评估方法
CN108038292A (zh) 一种基于双代理技术的高效自适应采样方法
Hariharan et al. Helicopter Aerodynamic Modeling of S-76 Rotor with Tip-Shape Variations: Review of AIAA Standardized Hover Evaluations
CN109447026B (zh) 一种基于支持向量机的桩身完整性类别自动识别方法
CN104298706A (zh) 一种基于数据挖掘的输电杆塔塔材实际强度计算方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20170818

Address after: 300392 Tianjin City Huayuan Industrial Zone Branch Road No. 15 Building No. 5 Room 501

Patentee after: Zhongke Guofeng science and Technology Co Ltd

Address before: 100190 Beijing Haidian District North Fourth Ring Road West, No. 11 A202

Patentee before: Institute of Engineering Thermophysics, Chinese Academy of Sciences

TR01 Transfer of patent right