CN113931807B - 一种风电叶片运行攻角测量方法 - Google Patents
一种风电叶片运行攻角测量方法 Download PDFInfo
- Publication number
- CN113931807B CN113931807B CN202110979566.3A CN202110979566A CN113931807B CN 113931807 B CN113931807 B CN 113931807B CN 202110979566 A CN202110979566 A CN 202110979566A CN 113931807 B CN113931807 B CN 113931807B
- Authority
- CN
- China
- Prior art keywords
- blade
- angle
- attack
- probe
- airflow probe
- 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.)
- Active
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 11
- 239000000523 sample Substances 0.000 claims abstract description 59
- 238000005259 measurement Methods 0.000 claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 238000000034 method Methods 0.000 claims description 14
- 238000009434 installation Methods 0.000 claims description 10
- 239000002131 composite material Substances 0.000 claims description 5
- 238000009530 blood pressure measurement Methods 0.000 claims description 4
- 230000006698 induction Effects 0.000 abstract description 23
- 239000000306 component Substances 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D1/00—Wind motors with rotation axis substantially parallel to the air flow entering the rotor
- F03D1/06—Rotors
- F03D1/065—Rotors characterised by their construction elements
- F03D1/0675—Rotors characterised by their construction elements of the blades
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种风电叶片运行攻角测量方法,包括叶片前缘相对入流速度矢量测量,叶片翼型气动力测量,基于升力环量修正的攻角计算三个步骤,其中叶片前缘相对入流速度矢量测量通过将气流探针安装于叶片前缘,随叶片共同旋转测得,叶片翼型气动力测量指翼型弦线坐标系下的法向力与切向力测量,基于升力环量修正的攻角计算用来修正气流探针直接测量结果中包含的叶片升力环量诱导速度。
Description
技术领域
本发明属于风电叶片空气动力学技术领域,涉及一种空气动力学参数测量方法,更具体涉及一种风电叶片运行攻角测量方法。
背景技术
风电叶片是将风能转换为机械能的核心部件,其气动性能决定了风电机组的风能利用效率(功率系数)。风电叶片气动外形是由一些列不同厚度及弦长的翼型按照一定安装扭角分布积叠而成,评价翼型气动性能的关键参数是攻角和气动力。
尽管风电产业发展迅速,装机体量庞大,但风电叶片在实际运行过程中的气动性能一直没有得到很好的评估与验证,主要原因之一是目前还缺乏有效的测量方法获得风电叶片运行攻角,其难点在于,根据风电机组空气动力学理论可知,风电叶片翼型攻角与风轮远前方来流风速、风轮转速、安装扭角、轴向及周向诱导因子有关,而其中轴向及周向诱导因子是无法直接测量的未知量,因此仅测量来流风速、风轮转速并不能计算得到风电叶片攻角。
发明内容
针对现有技术的以上问题,本发明提出一种风电叶片运行攻角测量方法,提出将气流探针安装于叶片前缘,直接测量叶片所感受的相对入流速度矢量,该速度矢量是远前方来流速度、风轮转速、轴向和周向诱导速度、叶片升力环量诱导速度的合成矢量,其中与叶片运行攻角相关的速度矢量包括远前方来流速度、风轮转速、轴向和周向诱导速度,因此叶片升力环量诱导速度需要计算去除,通过提出相对入流速度矢量和气动力同时测量的方案,并采用基于升力环量修正的计算方法,可以准确去除叶片升力环量诱导速度的影响,从而得到了叶片运行攻角。
本发明为解决其技术问题所采取的技术方案为:
一种风电叶片运行攻角测量方法,其特征在于,所述方法至少依次包括叶片前缘相对入流速度矢量测量、叶片翼型气动力测量、基于升力环量修正的攻角计算等三个步骤,其中,
SS1.叶片前缘相对入流速度矢量测量
在风电叶片的前缘位置安装一沿弦长方向延伸的气流探针,所述气流探针随风电叶片共同旋转,气流探针直接测量的速度为叶片前缘相对入流速度矢量,包括合成速度值Vp、方向角度αp;
SS2.叶片翼型气动力测量
在所述气流探针的流动干扰区域之外的风电叶片表面上选择至少一个待测叶片翼型,并在各所述待测叶片翼型的压力面和吸力面上分别从前缘至尾缘布置多个表面压力测量点以测量所述待测叶片翼型的表面压力分布,并根据各所述待测叶片翼型的弦线坐标系分别积分得到沿翼型弦线的切向力和法向力,并将各所述待测叶片翼型的切向力和法向力的平均值作为所述气流探针安装位置的叶片翼型的切向力A和法向力N;
SS3.基于升力环量修正的攻角计算
所述升力环量修正指由于气流探针安装在叶片附近,其直接测量的速度中包含了叶片升力环量诱导产生的速度,因此气流探针直接测量的速度矢量需要进行升力环量修正,具体修正方法如下子步骤:
SS3.1.将步骤SS1气流探针所测得的叶片前缘相对入流速度矢量的方向角度αp作为风电叶片运行攻角的初始值α0;
SS3.2.根据步骤SS3.1所确定的风电叶片运行攻角的初始值α0和步骤SS2所测的所述气流探针安装位置的叶片翼型的法向力N和切向力A,计算所述气流探针安装位置的叶片翼型的升力L,其计算式为L=N·cosα0-Asinα0;
SS3.6.如果步骤SS3.5计算得到的更新攻角αe与步骤SS3.1初始攻角α0的差值大于设定阈值,则将步骤SS3.5计算得到的更新攻角αe作为步骤SS3.1初始攻角α0,再次重复步骤SS3.1~SS3.5,直到步骤SS3.5计算得到的更新攻角αe与步骤SS3.1初始攻角α0的差值小于设定阈值,则步骤SS3.5更新攻角αe为最终所测攻角,运动攻角测量与计算完成。
优选地,上述步骤SS1中的气流探针可以采用多孔气流探针或多分量热线风速探针等体积紧凑、对下游流动干扰较小的空气流速矢量传感器。
优选地,上述步骤SS2中叶片翼型气动力测量可以采用表面压力分布测量方法或者气动天平测力方法。
优选地,上述步骤SS2中待测的叶片翼型可以选取一个或者多个。
优选地,上述步骤SS2中待测的叶片翼型要处于气流探针的流动干扰区域之外,但应尽可能接近气流探针安装位置,以得到更高的测量精度。
优选地,上述步骤SS2中,在所述气流探针的安装位置的展向上的两侧,分别选定所述气流探针长度l的一半距离作为要进行气动力测量的两个待测叶片翼型,在各所述待测叶片翼型的压力面和吸力面上分别从前缘至尾缘布置多个表面压力测量点以测量所述待测叶片翼型的表面压力分布,并根据各所述待测叶片翼型的弦线坐标系分别积分得到沿翼型弦线的切向力和法向力,并将各所述待测叶片翼型的切向力和法向力的平均值作为所述气流探针安装位置的叶片翼型的切向力A和法向力N。
通过以上技术方案可知,本发明提出将气流探针安装于叶片前缘,直接测量叶片所感受的相对入流速度矢量,该速度矢量是远前方来流速度、风轮转速、轴向和周向诱导速度、叶片升力环量诱导速度的合成矢量,其中与叶片运行攻角相关的速度矢量包括远前方来流速度、风轮转速、轴向和周向诱导速度,因此叶片升力环量诱导速度需要计算去除。为此本发明提出相对入流速度矢量和气动力同时测量的方案,并采用基于升力环量修正的计算方法,可以准确去除叶片升力环量诱导速度的影响,从而得到了叶片运行攻角。可见,本发明所述方法有效规避了风轮轴向和周向诱导因子测量的难题,通过间接方法测得了叶片运行攻角,进一步通过所测的叶片运行攻角还可以推算风轮轴向和轴向诱导因子。
附图说明
图1为本发明的风电叶片运行攻角测量原理示意图;
附图标记说明:
气流探针1,风电叶片2,待测叶片翼型3,表面压力测量点4
具体实施方式
为使本发明目的、技术方案更加清楚明白,以具体实施案例,并参照附图,对本发明进一步详细说明。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例,而不是全部的实施例,旨在用于解释本发明,而不能理解为对本发明的限制。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图一所示,本发明的风电叶片运行攻角测量方法在实施时,主要包括叶片前缘相对入流速度矢量测量,叶片翼型气动力测量,基于升力环量修正的攻角计算等三个步骤,各步骤的具体实施方式如下:
1.叶片前缘相对入流速度矢量测量
所述速度矢量包括速度大小和方向角度,所述叶片前缘相对入流速度指将气流探针安装于叶片前缘,随叶片共同旋转,气流探针直接测量的速度为叶片前缘相对入流速度。具体而言,如图1所示,将七孔气流探针1安装于叶片2前缘,随叶片2共同旋转,七孔气流探针1直接测量的速度为叶片前缘相对入流速度,具体包括合成速度大小Vp、方向角度αp。
2.叶片翼型气动力测量
所述叶片翼型指气流探针的流动干扰区域之外的翼型,由于攻角未知,所述气动力测量指翼型弦线坐标系下的法向力与切向力测量。具体而言,如图1所示,在七孔气流探针1安装位置两侧,选定七孔气流探针1长度l的一半距离作为要进行气动力测量的两个叶片翼型3,在各自叶片翼型3上布置表面压力测量点4,测量表面压力分布,进而按各自叶片翼型3坐标系分别积分得到沿各自翼型弦线的切向力和法向力,并将两个叶片翼型3的切向力和法向力平均值作为七孔气流探针1安装位置翼型的切向力A和法向力N。
3.基于升力环量修正的攻角计算
所述升力环量修正指由于气流探针安装在叶片附近,其直接测量的速度中包含了叶片升力环量诱导产生的速度,因此气流探针直接测量的速度矢量需要进行升力环量修正。具体修正方法如下:
3.1.将步骤1七孔气流探针1所测的方向角度αp作为攻角的初始值α0;
3.2.根据步骤3.1攻角初始值α0和步骤2所测的法向力N和切向力A,计算升力L=N·cosα0-Asinα0;
3.6.如果步骤3.5更新攻角αe与步骤3.1初始攻角α0差值大于0.1°,则将步骤3.5更新攻角αe作为步骤3.1新的初始攻角α0,再次重复步骤3.1~3.5,直到步骤3.5更新攻角αe与步骤3.1初始攻角α0差值小于0.1°,则步骤3.5更新攻角αe为最终所测攻角,运动攻角测量与计算完成。
通过以上技术方案可知,本发明提出将气流探针安装于叶片前缘,直接测量叶片所感受的相对入流速度矢量,该速度矢量是远前方来流速度、风轮转速、轴向和周向诱导速度、叶片升力环量诱导速度的合成矢量,其中与叶片运行攻角相关的速度矢量包括远前方来流速度、风轮转速、轴向和周向诱导速度,因此叶片升力环量诱导速度需要计算去除。为此本发明提出相对入流速度矢量和气动力同时测量的方案,并采用基于升力环量修正的计算方法,可以准确去除叶片升力环量诱导速度的影响,从而得到了叶片运行攻角。可见,本发明所述方法有效规避了风轮轴向和周向诱导因子测量的难题,通过间接方法测得了叶片运行攻角,进一步通过所测的叶片运行攻角还可以推算风轮轴向和轴向诱导因子。
以上所述仅为本发明的一个实施例,并不用以限制本发明,凡在本发明的思路和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的范围之内。
Claims (3)
1.一种风电叶片运行攻角测量方法,其特征在于,所述方法至少依次包括叶片前缘相对入流速度矢量测量、叶片翼型气动力测量、基于升力环量修正的攻角计算这三个步骤,其中,
SS1.叶片前缘相对入流速度矢量测量
在风电叶片的前缘位置安装一沿弦长方向延伸的气流探针,所述气流探针随风电叶片共同旋转,气流探针直接测量的速度为叶片前缘相对入流速度矢量,包括合成速度值Vp、方向角度αp;
SS2.叶片翼型气动力测量
在所述气流探针的流动干扰区域之外且在所述气流探针的安装位置的展向上的两侧的风电叶片表面上,分别选定所述气流探针长度l的一半距离作为要进行气动力测量的两个待测叶片翼型,并在各所述待测叶片翼型的压力面和吸力面上分别从前缘至尾缘布置多个表面压力测量点以测量所述待测叶片翼型的表面压力分布,并根据各所述待测叶片翼型的弦线坐标系分别积分得到沿翼型弦线的切向力和法向力,并将各所述待测叶片翼型的切向力和法向力的平均值作为所述气流探针安装位置的叶片翼型的切向力A和法向力N;
SS3.基于升力环量修正的攻角计算
所述升力环量修正指由于气流探针安装在叶片附近,其直接测量的速度中包含了叶片升力环量诱导产生的速度,因此气流探针直接测量的速度矢量需要进行升力环量修正,具体修正方法如下子步骤:
SS3.1.将步骤SS1气流探针所测得的叶片前缘相对入流速度矢量的方向角度αp作为风电叶片运行攻角的初始值α0;
SS3.2.根据步骤SS3.1所确定的风电叶片运行攻角的初始值α0和步骤SS2所测的所述气流探针安装位置的叶片翼型的法向力N和切向力A,计算所述气流探针安装位置的叶片翼型的升力L,其计算式为L=N·cosα0-Asinα0;
SS3.6.如果步骤SS3.5计算得到的更新攻角αe与步骤SS3.1初始攻角α0的差值大于设定阈值,则将步骤SS3.5计算得到的更新攻角αe作为步骤SS3.1初始攻角α0,再次重复步骤SS3.1~SS3.5,直到步骤SS3.5计算得到的更新攻角αe与步骤SS3.1初始攻角α0的差值小于设定阈值,则步骤SS3.5更新攻角αe为最终所测攻角,运动攻角测量与计算完成。
2.根据权利要求1所述的风电叶片运行攻角测量方法,其特征在于,上述步骤SS1中的气流探针采用多孔气流探针或多分量热线风速探针。
3.根据权利要求1所述的风电叶片运行攻角测量方法,其特征在于,上述步骤SS2中叶片翼型气动力测量采用表面压力分布测量方法或者气动天平测力方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110979566.3A CN113931807B (zh) | 2021-08-25 | 2021-08-25 | 一种风电叶片运行攻角测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110979566.3A CN113931807B (zh) | 2021-08-25 | 2021-08-25 | 一种风电叶片运行攻角测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113931807A CN113931807A (zh) | 2022-01-14 |
CN113931807B true CN113931807B (zh) | 2023-04-21 |
Family
ID=79274437
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110979566.3A Active CN113931807B (zh) | 2021-08-25 | 2021-08-25 | 一种风电叶片运行攻角测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113931807B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1217767A (zh) * | 1995-07-06 | 1999-05-26 | 奥拉夫·沙茨 | 垂直轴涡流风力涡轮机 |
EP2536945A2 (en) * | 2010-02-16 | 2012-12-26 | Technion Research and Development Foundation, Ltd. | Flow control on a vertical axis wind turbine (vawt) |
DE102013103150A1 (de) * | 2013-03-27 | 2014-10-02 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Verfahren zum Ermitteln eines Anstellwinkels |
CN110532604A (zh) * | 2019-07-23 | 2019-12-03 | 南京航空航天大学 | 一种动态失速状态下带后缘小翼桨叶气动载荷的计算方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6068446A (en) * | 1997-11-20 | 2000-05-30 | Midwest Research Institute | Airfoils for wind turbine |
US20050197811A1 (en) * | 2004-02-02 | 2005-09-08 | Ariyur Kartik B. | Prediction of dynamic ground effect forces for fixed wing aircraft |
EP2253836A1 (en) * | 2009-05-18 | 2010-11-24 | Lm Glasfiber A/S | Wind turbine blade |
DK2569535T3 (en) * | 2010-05-10 | 2016-12-05 | Univ Darmstadt Tech | Invention relating to the blade, especially for wind projects |
JP5881491B2 (ja) * | 2011-05-02 | 2016-03-09 | 株式会社東芝 | 風力発電システムおよびその制御方法 |
CN102322407B (zh) * | 2011-08-30 | 2013-06-12 | 重庆大学 | 风力机叶片气动外形协同设计方法 |
US9372201B2 (en) * | 2014-03-31 | 2016-06-21 | Alstom Renewable Technologies | Yaw and pitch angles |
CN107169163A (zh) * | 2017-04-13 | 2017-09-15 | 南京航空航天大学 | 一种适用于机翼气动参数分布实时计算的解耦算法 |
CN108069020A (zh) * | 2017-11-21 | 2018-05-25 | 中国人民解放军空军工程大学 | 飞行器升力优化系统 |
CN109751201B (zh) * | 2019-02-18 | 2020-12-15 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种风力机涡尾迹修正方法 |
-
2021
- 2021-08-25 CN CN202110979566.3A patent/CN113931807B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1217767A (zh) * | 1995-07-06 | 1999-05-26 | 奥拉夫·沙茨 | 垂直轴涡流风力涡轮机 |
EP2536945A2 (en) * | 2010-02-16 | 2012-12-26 | Technion Research and Development Foundation, Ltd. | Flow control on a vertical axis wind turbine (vawt) |
DE102013103150A1 (de) * | 2013-03-27 | 2014-10-02 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Verfahren zum Ermitteln eines Anstellwinkels |
CN110532604A (zh) * | 2019-07-23 | 2019-12-03 | 南京航空航天大学 | 一种动态失速状态下带后缘小翼桨叶气动载荷的计算方法 |
Non-Patent Citations (2)
Title |
---|
芮晓明 ; 马志勇 ; 康传明 ; .并网型垂直轴风力机的基本构成与气动特性.太阳能.2007,(02),全文. * |
陈洪胜 ; .后掠改型风力机叶片叶素的入流特性.科学技术创新.2019,(第13期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN113931807A (zh) | 2022-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Inoue et al. | Behavior of tip leakage flow behind an axial compressor rotor | |
Hathaway et al. | Experimental and computational investigation of the NASA low-speed centrifugal compressor flow field | |
US10202965B2 (en) | Method and system for estimating the surface condition of a rotating blade | |
Sareen et al. | Effects of leading-edge protection tape on wind turbine blade performance | |
CN101725466A (zh) | 使用压力传感器的失速检测 | |
Wagner et al. | Inlet boundary layer effects in an axial compressor rotor: part I—blade-to-blade effects | |
Meyer et al. | A parameter study on the influence of fillets on the compressor cascade performance | |
CN108241767A (zh) | 一种基于叶素理论的海上风机叶片气动载荷分析方法 | |
Walker et al. | Periodic transition on an axial compressor stator: Incidence and clocking effects: Part I—Experimental data | |
CN216342583U (zh) | 一种风电叶片运行攻角测量装置 | |
CN113931807B (zh) | 一种风电叶片运行攻角测量方法 | |
Hamlaoui et al. | Improved bem method for hawt performance predictions | |
Sant et al. | Estimating the angle of attack from blade pressure measurements on the National Renewable Energy Laboratory phase VI rotor using a free wake vortex model: yawed conditions | |
CN114266202A (zh) | 一种模拟风力机尾流的修正致动线模型方法 | |
CN106246459B (zh) | 一种风力机导风筒 | |
Maeda et al. | Surface pressure measurement on a rotating blade of field horizontal axis wind turbine in yawed condition | |
Zeschky et al. | Effects of stator wakes and spanwise nonuniform inlet conditions on the rotor flow of an axial turbine stage | |
CN206071794U (zh) | 一种风力机导风筒 | |
DK2940292T3 (en) | Device for a rotor blade of a wind turbine | |
Soto-Valle et al. | On the Influence of trip strips on Rotor Blade Measurements | |
Pullan et al. | The Influence of Sweep on Axial Flow Turbine Aerodynamics at Mid-Span | |
Soto-Valle et al. | Development and assessment of a method to determine the angle of attack on an operating wind turbine by matching onboard pressure measurements with panel method simulations | |
Cyrus | Effect of the inlet velocity profile in the three-dimensional flow in a rear axial compressor stage | |
Davino et al. | Turbulence characteristics in the annulus-wall boundary layer and wake mixing region of a compressor rotor exit | |
Nishi et al. | Rotor Design Method for a Propeller-Type Wind Turbine that Combines Optimization Method with the Blade Element Momentum Theory |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |