CN109325274B - 基于三因子拟合积分法的风力机风轮气动性能评价方法 - Google Patents

基于三因子拟合积分法的风力机风轮气动性能评价方法 Download PDF

Info

Publication number
CN109325274B
CN109325274B CN201811038027.4A CN201811038027A CN109325274B CN 109325274 B CN109325274 B CN 109325274B CN 201811038027 A CN201811038027 A CN 201811038027A CN 109325274 B CN109325274 B CN 109325274B
Authority
CN
China
Prior art keywords
wind turbine
blade
coefficient
interference factor
wind
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
Application number
CN201811038027.4A
Other languages
English (en)
Other versions
CN109325274A (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.)
Yinchuan College China of CUMT
Original Assignee
Yinchuan College China of CUMT
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 Yinchuan College China of CUMT filed Critical Yinchuan College China of CUMT
Priority to CN201811038027.4A priority Critical patent/CN109325274B/zh
Publication of CN109325274A publication Critical patent/CN109325274A/zh
Application granted granted Critical
Publication of CN109325274B publication Critical patent/CN109325274B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Wind Motors (AREA)

Abstract

本发明涉及一种基于三因子拟合积分法的风力机风轮气动性能评价方法，采用拟合方式对轴向推力、转矩和功率三者无量纲化后对应的系数进行计算，考虑到了传统的计算方法未考虑的被积函数是不连续的，严格意义上讲是不能积分的问题，以及计算出来各叶素面的参数均不同直接用积分的形式来计算性能结果偏差大的问题。本方法相比传统的计算风力机性能的计算模型，更为科学精确，处理得到的结果可靠性更高。

Description

基于三因子拟合积分法的风力机风轮气动性能评价方法
技术领域
本发明涉及风力发电技术领域，具体涉及一种基于三因子拟合积分法的风力机风轮气动性能评价方法。
背景技术
目前，风力发电是当前世界生态经济中最具活力、增长最快的清洁能源产业。随着大型风电机组的迅猛发展，离网型小型风力机的发展也在蒸蒸日上，保证其具有更好的经济性和更高的安全稳定性尤为重要。
近年来，中小型风力机因调速和调向系统不灵(市场中没有调速装置的风力机屡见不鲜)，暴露出一些严重问题，事故或故障频繁发生：如小风速时，风轮经常不对风，机头转动困难；风大时(如风速超12m/s以上)，风轮不能及时偏转限速，使风轮长时间工作在超速旋转状态下，致使风力机工作稳定性变坏，在远短于设计寿命期内频繁发生疲劳损伤事故(如发电机烧毁、风力机叶根损伤和断裂及由其引起的中小型风力机倒塔事故)，给风力发电工作带来巨大损失。为此，有关风力机各参数分析与性能评价工作开始展开，但是，由于对风轮在非设计工况下的性能评价不准确且评价方法不科学，导致调速机构的的结构上设计及性能优化工作停滞不前。
由此可见，如何提供一种准确科学的风力机风轮气动性能评价方法成为本领域技术人员亟待解决的技术问题。
发明内容
本发明所要解决的技术问题是提供一种准确科学的风力机风轮气动性能评价方法。
本发明解决上述技术问题的技术方案如下：一种基于三因子拟合积分法的风力机风轮气动性能评价方法，该方法包括以下步骤：
利用Profili软件和Matlab编程软件结合各翼型的气动数据完成风力机叶片翼型的选取工作，得到修正后叶片的弦长C和扭角θ；
根据修正后风力机的弦长C、扭角θ和Profili翼型数据库中所含数据，计算轴向干涉因子a、周向干涉因子b、叶梢损失系数F及入流角
Figure BDA0001791386960000023
构建风力机的轴向推力T的数学模型，数学模型公式如下：
dT＝4πrV1 2a(1-a)ρFdr   (1)
其中，T为轴向推力，r为叶片各叶素面距轮毂中心的距离；V1为受扰动后的风速；a为轴向干涉因子；ρ为空气密度；F为叶梢损失系数；
将轴向推力T无量纲化，得到其对应的系数CT，计算公式为：
Figure BDA0001791386960000021
其中，CT为轴向推力T无量纲化后对应的系数；λ0为设计叶尖速比；λ为尖速比；T为轴向推力；V1为受扰动后的风速；a为轴向干涉因子；ρ表示空气密度；F表示叶梢损失系数；R为风轮半径；
构建风力机的转矩M的数学模型，数学模型公式如下：
dM＝4πr3V1Ωb(1-a)ρFdr   (3)
其中，M为风力机的转矩；r为叶片各叶素面距轮毂中心的距离；V1为受扰动后的风速；a为轴向干涉因子；b为周向干涉因子；ρ为空气密度；F为叶梢损失系数；
将风力机的转矩M无量纲化，得到其对应的系数CM，计算公式为：
Figure BDA0001791386960000022
其中，CM为风力机的转矩M无量纲化后对应的系数；λ0为设计叶尖速比；λ为尖速比；M为风力机的转矩；V1为受扰动后的风速；a为轴向干涉因子；b为周向干涉因子；ρ表示空气密度；F表示叶梢损失系数；R为风轮半径；
构建风力机的功率P的数学模型，数学模型公式如下：
dP＝ΩdM＝4πr3V1Ω2b(1-a)ρFdr   (5)
其中，P为风力机的功率；M为风力机的转矩；r为叶片各叶素面距轮毂中心的距离；V1为受扰动后的风速；a为轴向干涉因子；b为周向干涉因子；ρ为空气密度；F为叶梢损失系数；
将风力机的功率P无量纲化，得到其对应的系数CP，计算公式为：
Figure BDA0001791386960000031
其中，CP为风力机的功率P无量纲化后对应的系数；CM为风力机的转矩M无量纲化后对应的系数；λ0为设计叶尖速比；λ为尖速比；a为轴向干涉因子；b为周向干涉因子；F为叶梢损失系数；
由于轴向干涉因子a、周向干涉因子b及叶梢损失系数F均是叶片各叶素面距轮毂中心的距离r的单值函数，利用常规的拟合方法，得到如下f(r)的多项表达式：
bF(1-a)＝f1(r)   (7)
其中，a为轴向干涉因子；b为周向干涉因子；F为叶梢损失系数；r为叶片各叶素面距轮毂中心的距离；
将公式(7)代入公式(6)中，得到如下公式：
Figure BDA0001791386960000032
其中，λ0为设计叶尖速比；r为叶片各叶素面距轮毂中心的距离；R为风轮半径；
同理，将公式(7)分别代入公式(2)和公式(4)，得到轴向推力T无量纲化后对应的系数CT以及风力机的转矩M无量纲化后对应的系数CM的数学模型公式，完成风轮性能的计算。
本发明的有益效果是：风力机的三个性能指标轴向推力、转矩和功率均是轴向干涉因子、周向干涉因子及叶梢损失系数的函数，计算性能就是求取这三个参数的过程，本方法采用拟合方式对轴向推力、转矩和功率三者无量纲化后对应的系数进行计算，考虑到了传统的计算方法未考虑的被积函数是不连续的，严格意义上讲是不能积分的问题，以及计算出来各叶素面的参数均不同直接用积分的形式来计算性能结果偏差大的问题。本方法依托“三因子拟合积分法”的算法，找到轴向干涉因子a与叶素面位置r/R的函数关系式，确定风轮推力系数，并非将其按照传统的做法取为常数，相比传统的计算风力机性能的计算模型，更为科学精确，处理得到的结果可靠性更高。
在上述技术方案的基础上，本发明还可以做如下改进。
进一步，利用Profili软件和Matlab编程软件结合各翼型的气动数据完成风力机叶片翼型的选取工作，涉及的具体步骤为：
(1)收集国内外翼型气动数据完整的翼型信息，将其导入Profili软件中，分析其在各个雷诺数Re下，升阻比随攻角的变化关系Alpha＝f(Re)，得到Profili翼型数据库；
(2)将各个雷诺数Re下升阻比随攻角的变化关系Alpha＝f(Re)的数据以记事本文件形式自动导入Matlab软件中，形成m文件的形式；
(3)基于Wilson法编制叶片的设计程序，将步骤(1)中筛选出来的各个翼型分别自动导入设计程序中，计算利用各种翼型设计得到的叶片的弦长和扭角数据；
(4)结合各翼型的实际气动数据，将步骤(2)得到的叶片的弦长和扭角数据，代入风轮气动性能计算程序中，计算出各种翼型对应的整个风轮的功率系数CP，选择功率系数CP最大的风轮对应的翼型作为被选择翼型。
进一步，所述入流角
Figure BDA0001791386960000051
的计算公式为：
Figure BDA0001791386960000052
其中，λ为尖速比；V1为受扰动后的风速；r为叶片各叶素面距轮毂中心的距离；a为轴向干涉因子；b为周向干涉因子。
进一步，轴向干涉因子a和周向干涉因子b的计算公式为：
Figure BDA0001791386960000053
Figure BDA0001791386960000054
其中，a为轴向干涉因子；b为周向干涉因子；
Figure BDA0001791386960000055
为入流角；r为叶片各叶素面距轮毂中心的距离；C为修正后的弦长；Cx为弦长C在x轴方向上的分量；Cy为弦长C在y方向上的分量；B为叶片数；F为叶梢损失系数。
进一步，弦长C在x轴方向上的分量Cx的计算公式和弦长C在y方向上的分量Cy的计算公式分别为：
Figure BDA0001791386960000056
Figure BDA0001791386960000057
其中，CL为升力系数；CD为阻力系数；
Figure BDA0001791386960000058
为入流角。
进一步，升力系数CL和阻力系数CD在Matlab中采用二维线性插值方法计算得到，计算公式分别为：
CL＝interp2(ALPHA，Re，CL，ALPHA，Re)   (14)
CD＝interp2(ALPHA，Re，CD，ALPHA，Re)   (15)
其中，ALPHA表示主程序中的攻角α，Re表示雷诺数。
进一步，攻角α的计算公式为：
Figure BDA0001791386960000059
其中，
Figure BDA00017913869600000510
为入流角，θ为修正后的扭角。
进一步，叶梢损失系数F的计算公式为：
Figure BDA0001791386960000061
其中，f为风轮扫掠面与塔架轴线的距离。
进一步，风轮扫掠面与塔架轴线的距离f的计算公式为：
Figure BDA0001791386960000062
其中，B为叶片数；
Figure BDA0001791386960000063
为入流角；r为叶片各叶素面距轮毂中心的距离；R为风轮半径。
采用上述进一步方案的有益效果是针对某一指定状态叶尖速比λ0，对于每个给定的r和λ，使用While循环迭代法可求取a和b，例如赋初值a＝0.3，b＝0.2，反复迭代直到收敛，最终得到a、b、F，进而拟合得到f(r)多项表达式，将其代入相应公式中，即可完成风轮的性能计算。对风轮的性能计算过程简便快捷，且计算结果精确可靠。
对本发明中涉及到的有关风力机的参数名称的具体含义做如下说明：
(1)风轮：风轮是叶片安装在轮毂上的总称；
(2)风轮直径：叶尖旋转圆的直径D；
(3)叶片：叶片是接受风能的基本部件；
(4)翼弦：前缘点与后缘点的连线，即弦长C。
附图说明
图1为本发明一种基于三因子拟合积分法的风力机风轮气动性能评价方法的方法流程图；
图2为本发明利用Profili软件和Matlab编程软件结合各翼型的气动数据完成风力机叶片翼型的选取工作的方法流程图；
图3为本发明风力机功率系数与叶尖速比的关系图；
图4为本发明风力机叶尖速比与推力系数的关系图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
如图1所示，本实施例提供了一种基于三因子拟合积分法的风力机风轮气动性能评价方法，该方法包括以下步骤：
S1:利用Profili软件和Matlab编程软件结合各翼型的气动数据完成风力机叶片翼型的选取工作，得到修正后叶片的弦长C和扭角θ；
S2:根据修正后风力机的弦长C、扭角θ和Profili翼型数据库中所含数据，计算轴向干涉因子a、周向干涉因子b、叶梢损失系数F及入流角
Figure BDA0001791386960000071
S3:构建风力机的轴向推力T的数学模型，数学模型公式如下：
dT＝4πrV1 2a(1-a)ρFdr   (1)
其中，T为轴向推力，r为叶片各叶素面距轮毂中心的距离；V1为受扰动后的风速；a为轴向干涉因子；ρ为空气密度；F为叶梢损失系数；
构建风力机的转矩M的数学模型，数学模型公式如下：
dM＝4πr3V1Ωb(1-a)ρFdr   (3)
其中，M为风力机的转矩；r为叶片各叶素面距轮毂中心的距离；V1为受扰动后的风速；a为轴向干涉因子；b为周向干涉因子；ρ为空气密度；F为叶梢损失系数；
构建风力机的功率P的数学模型，数学模型公式如下：
dP＝ΩdM＝4πr3V1Ω2b(1-a)ρFdr   (5)
其中，P为风力机的功率；M为风力机的转矩；r为叶片各叶素面距轮毂中心的距离；V1为受扰动后的风速；a为轴向干涉因子；b为周向干涉因子；ρ为空气密度；F为叶梢损失系数；
S4:将轴向推力T无量纲化，得到其对应的系数CT，计算公式为：
Figure BDA0001791386960000081
其中，CT为轴向推力T无量纲化后对应的系数；λ0为设计叶尖速比；λ为尖速比；T为轴向推力；V1为受扰动后的风速；a为轴向干涉因子；ρ表示空气密度；F表示叶梢损失系数；R为风轮半径；
将风力机的转矩M无量纲化，得到其对应的系数CM，计算公式为：
Figure BDA0001791386960000082
其中，CM为风力机的转矩M无量纲化后对应的系数；λ0为设计叶尖速比；λ为尖速比；M为风力机的转矩；V1为受扰动后的风速；a为轴向干涉因子；b为周向干涉因子；ρ表示空气密度；F表示叶梢损失系数；R为风轮半径；
将风力机的功率P无量纲化，得到其对应的系数CP，计算公式为：
Figure BDA0001791386960000083
其中，CP为风力机的功率P无量纲化后对应的系数；CM为风力机的转矩M无量纲化后对应的系数；λ0为设计叶尖速比；λ为尖速比；a为轴向干涉因子；b为周向干涉因子；F为叶梢损失系数；
S5:由于轴向干涉因子a、周向干涉因子b及叶梢损失系数F均是风轮半径r的单值函数，利用常规的拟合方法，得到如下f(r)的多项表达式：
bF(1-a)＝f1(r)   (7)
其中，a为轴向干涉因子；b为周向干涉因子；F为叶梢损失系数；r为叶片各叶素面距轮毂中心的距离；
S6:将公式(7)代入公式(6)中，得到如下公式：
Figure BDA0001791386960000091
其中，λ0为设计叶尖速比；r为叶片各叶素面距轮毂中心的距离；R为风轮半径；
同理，将公式(7)分别代入公式(2)和公式(4)，得到轴向推力T无量纲化后对应的系数CT以及风力机的转矩M无量纲化后对应的系数CM的数学模型公式，完成风轮性能的计算。
参见附图2，在一个具体的实施例中，利用Profili软件和Matlab编程软件结合各翼型的气动数据完成风力机叶片翼型的选取工作，涉及的具体步骤为：
a、收集国内外翼型气动数据完整的翼型信息，将其导入Profili软件中，分析其在各个雷诺数Re下，升阻比随攻角的变化关系Alpha＝f(Re)，得到Profili翼型数据库；
b、将各个雷诺数Re下升阻比随攻角的变化关系Alpha＝f(Re)的数据以记事本文件形式自动导入Matlab软件中，形成m文件的形式；
c、基于Wilson法编制叶片的设计程序，将步骤a中筛选出来的各个翼型分别自动导入设计程序中，计算利用各种翼型设计得到的叶片的弦长和扭角数据；
d、结合各翼型的实际气动数据，将步骤b得到的叶片的弦长和扭角数据，代入风轮气动性能计算程序中，计算出各种翼型对应的整个风轮的功率系数CP，选择功率系数CP最大的风轮对应的翼型作为被选择翼型。
具体地，入流角
Figure BDA0001791386960000092
的计算公式为：
Figure BDA0001791386960000093
其中，λ为尖速比；V1为受扰动后的风速；r为叶片各叶素面距轮毂中心的距离；a为轴向干涉因子；b为周向干涉因子。
具体地，轴向干涉因子a和周向干涉因子b的计算公式为：
Figure BDA0001791386960000101
Figure BDA0001791386960000102
其中，a为轴向干涉因子；b为周向干涉因子；
Figure BDA00017913869600001010
为入流角；r为叶片各叶素面距轮毂中心的距离；C为修正后的弦长；Cx为弦长C在x轴方向上的分量；Cy为弦长C在y方向上的分量；B为叶片数；F为叶梢损失系数。
具体地，弦长C在x轴方向上的分量Cx的计算公式和弦长C在y方向上的分量Cy的计算公式分别为：
Figure BDA0001791386960000104
Figure BDA0001791386960000105
其中，CL为升力系数；CD为阻力系数；
Figure BDA0001791386960000106
为入流角。
具体地，升力系数CL和阻力系数CD在Matlab中采用二维线性插值方法计算得到，计算公式分别为：
CL＝interp2(ALPHA，Re，CL，ALPHA，Re)   (14)
CD＝interp2(ALPHA，Re，CD，ALPHA，Re)   (15)
其中，ALPHA表示主程序中的攻角α，Re表示雷诺数。
上述数学模型中求解升力系数CL和阻力系数CD属于MATLAB中的二维线性插值问题，输入量为主程序中的攻角ALPHA和雷诺数Re，输出量为升力系数CL和阻力系数CD。
具体地，攻角α的计算公式为：
Figure BDA0001791386960000107
其中，
Figure BDA0001791386960000108
为入流角，θ为修正后的扭角。
具体地，叶梢损失系数F的计算公式为：
Figure BDA0001791386960000109
其中，f为风轮扫掠面与塔架轴线的距离。
具体地，风轮扫掠面与塔架轴线的距离f的计算公式为：
Figure BDA0001791386960000111
其中，B为叶片数；
Figure BDA0001791386960000112
为入流角；r为叶片各叶素面距轮毂中心的距离；R为风轮半径。
下面通过具体的实施例并结合相关数据对本方法做分析验证。
赋初值a＝0.3，b＝0.2，反复迭代直到收敛，最终得到轴向干涉因子a、周向干涉因子b以及叶梢损失系数F的值，进而可计算其性能。
对计算得到的数据进行整理，得到风力机功率系数与叶尖速比的关系以及风力机叶尖速比与推力系数的关系，具体可参见附图3和附图4，由图3和图4可以看出，在有效叶尖速比范围内，功率系数都在0.3以上且变化趋势较平缓，在额定尖速比λ为6.0时，功率系数为0.474，推力系数为0.78，达到设计要求。
本实施例提供的方法采用拟合方式对轴向推力、转矩和功率三者无量纲化后对应的系数进行计算，考虑到了传统的计算方法未考虑的被积函数是不连续的，严格意义上讲是不能积分的问题，以及计算出来各叶素面的参数均不同直接用积分的形式来计算性能结果偏差大的问题。本方法相比传统的计算风力机性能的计算模型，更为科学精确，处理得到的结果可靠性更高。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种基于三因子拟合积分法的风力机风轮气动性能评价方法，其特征在于，包括以下步骤：
利用Profili软件和Matlab编程软件结合各翼型的气动数据完成风力机叶片翼型的选取工作，得到修正后叶片的弦长C和扭角θ；
根据修正后风力机的弦长C、扭角θ和Profili翼型数据库中所含数据，计算轴向干涉因子a、周向干涉因子b、叶梢损失系数F及入流角
Figure FDA0004124251960000013
构建风力机的轴向推力T的数学模型，数学模型公式如下：
Figure FDA0004124251960000014
其中，T为轴向推力，r为叶片各叶素面距轮毂中心的距离；V1为受扰动后的风速；a为轴向干涉因子；ρ为空气密度；F为叶梢损失系数；
将轴向推力T无量纲化，得到其对应的系数CT，计算公式为：
Figure FDA0004124251960000011
其中，CT为轴向推力T无量纲化后对应的系数；λ0为设计叶尖速比；λ为尖速比；T为轴向推力，R为风轮半径；V1为受扰动后的风速；a为轴向干涉因子；ρ表示空气密度；F表示叶梢损失系数；
构建风力机的转矩M的数学模型，数学模型公式如下：
dM＝4πr3V1Ωb(1-a)ρFdr     (3)
其中，M为风力机的转矩；r为叶片各叶素面距轮毂中心的距离；V1为受扰动后的风速；a为轴向干涉因子；b为周向干涉因子；ρ为空气密度；F为叶梢损失系数；
将风力机的转矩M无量纲化，得到其对应的系数CM，计算公式为：
Figure FDA0004124251960000012
其中，CM为风力机的转矩M无量纲化后对应的系数；λ0为设计叶尖速比；λ为尖速比；M为风力机的转矩；R为风轮半径；V1为受扰动后的风速；a为轴向干涉因子；b为周向干涉因子；ρ表示空气密度；F表示叶梢损失系数；
构建风力机的功率P的数学模型，数学模型公式如下：
dP＝ΩdM＝4πr3V1Ω2b(1-a)ρFdr     (5)
其中，P为风力机的功率；M为风力机的转矩；r为叶片各叶素面距轮毂中心的距离；V1为受扰动后的风速；a为轴向干涉因子；b为周向干涉因子；ρ为空气密度；F为叶梢损失系数；
将风力机的功率P无量纲化，得到其对应的系数CP，计算公式为：
Figure FDA0004124251960000021
其中，CP为风力机的功率P无量纲化后对应的系数；CM为风力机的转矩M无量纲化后对应的系数；λ0为设计叶尖速比；λ为尖速比；a为轴向干涉因子；b为周向干涉因子；F为叶梢损失系数；
由于轴向干涉因子a、周向干涉因子b及叶梢损失系数F均是r的单值函数，利用常规的拟合方法，得到如下f1(r)的多项表达式：
bF(1-a)＝f1(r)    (7)
其中，a为轴向干涉因子；b为周向干涉因子；F为叶梢损失系数；r为叶片各叶素面距轮毂中心的距离；
将公式(7)代入公式(6)中，得到如下公式：
Figure FDA0004124251960000022
其中，λ0为设计叶尖速比；R为风轮半径；r为叶片各叶素面距轮毂中心的距离；
同理，将公式(7)分别代入公式(2)和公式(4)，得到轴向推力T无量纲化后对应的系数CT以及风力机的转矩M无量纲化后对应的系数CM的数学模型公式，完成风轮性能的计算。
2.根据权利要求1所述一种基于三因子拟合积分法的风力机风轮气动性能评价方法，其特征在于，利用Profili软件和Matlab编程软件结合各翼型的气动数据完成风力机叶片翼型的选取工作，涉及的具体步骤为：
(1)收集国内外翼型气动数据完整的翼型信息，将其导入Profili软件中，分析其在各个雷诺数Re下，升阻比随攻角的变化关系Alpha＝f(Re)，得到Profili翼型数据库；
(2)将各个雷诺数Re下升阻比随攻角的变化关系Alpha＝f(Re)的数据以记事本文件形式自动导入Matlab软件中，形成m文件的形式；
(3)基于Wilson法编制叶片的设计程序，将步骤(1)中筛选出来的各个翼型分别自动导入设计程序中，计算利用各种翼型设计得到的叶片的弦长和扭角数据；
(4)结合各翼型的实际气动数据，将步骤(3)得到的叶片的弦长和扭角数据，代入风轮气动性能计算程序中，计算出各种翼型对应的整个风轮的功率系数CP，选择功率系数CP最大的风轮对应的翼型作为被选择翼型。
3.根据权利要求1所述一种基于三因子拟合积分法的风力机风轮气动性能评价方法，其特征在于，所述入流角
Figure FDA0004124251960000031
的计算公式为：
Figure FDA0004124251960000032
其中，λ为尖速比；V1为受扰动后的风速；r为叶片各叶素面距轮毂中心的距离；a为轴向干涉因子；b为周向干涉因子。
4.根据权利要求3所述一种基于三因子拟合积分法的风力机风轮气动性能评价方法，其特征在于，轴向干涉因子a和周向干涉因子b的计算公式为：
Figure FDA0004124251960000033
Figure FDA0004124251960000034
其中，a为轴向干涉因子；b为周向干涉因子；
Figure FDA0004124251960000035
为入流角；r为叶片各叶素面距轮毂中心的距离；C为修正后的弦长；Cx为弦长C在x轴方向上的分量；Cy为弦长C在y方向上的分量；B为叶片数；F为叶梢损失系数。
5.根据权利要求4所述一种基于三因子拟合积分法的风力机风轮气动性能评价方法，其特征在于，弦长C在x轴方向上的分量Cx的计算公式和弦长C在y方向上的分量Cy的计算公式分别为：
Figure FDA0004124251960000036
Figure FDA0004124251960000037
其中，CL为升力系数；CD为阻力系数；
Figure FDA0004124251960000038
为入流角。
6.根据权利要求5所述一种基于三因子拟合积分法的风力机风轮气动性能评价方法，其特征在于，升力系数CL和阻力系数CD在Matlab中采用二维线性插值方法计算得到，计算公式分别为：
CL＝interp2(ALPHA，Re，CL，ALPHA，Re)   (14)
CD＝interp2(ALPHA，Re，CD，ALPHA，Re)      (15)
其中，ALPHA表示主程序中的攻角α，Re表示雷诺数。
7.根据权利要求6所述一种基于三因子拟合积分法的风力机风轮气动性能评价方法，其特征在于，攻角α的计算公式为：
Figure FDA0004124251960000041
其中，
Figure FDA0004124251960000042
为入流角，θ为修正后的扭角。
8.根据权利要求4所述一种基于三因子拟合积分法的风力机风轮气动性能评价方法，其特征在于，叶梢损失系数F的计算公式为：
Figure FDA0004124251960000043
其中，f为风轮扫掠面与塔架轴线的距离。
9.根据权利要求8所述一种基于三因子拟合积分法的风力机风轮气动性能评价方法，其特征在于，风轮扫掠面与塔架轴线的距离f的计算公式为：
Figure FDA0004124251960000044
其中，B为叶片数；
Figure FDA0004124251960000045
为入流角；r为叶片各叶素面距轮毂中心的距离；R为风轮半径。
CN201811038027.4A 2018-09-06 2018-09-06 基于三因子拟合积分法的风力机风轮气动性能评价方法 Active CN109325274B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811038027.4A CN109325274B (zh) 2018-09-06 2018-09-06 基于三因子拟合积分法的风力机风轮气动性能评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811038027.4A CN109325274B (zh) 2018-09-06 2018-09-06 基于三因子拟合积分法的风力机风轮气动性能评价方法

Publications (2)

Publication Number Publication Date
CN109325274A CN109325274A (zh) 2019-02-12
CN109325274B true CN109325274B (zh) 2023-04-25

Family

ID=65263610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811038027.4A Active CN109325274B (zh) 2018-09-06 2018-09-06 基于三因子拟合积分法的风力机风轮气动性能评价方法

Country Status (1)

Country Link
CN (1) CN109325274B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112084594B (zh) * 2020-09-07 2023-03-24 河海大学 一种改进的潮流能水轮机性能预测方法
CN116628895A (zh) * 2023-07-20 2023-08-22 北京航空航天大学 一种带分流叶片的斜流或径流压气机叶轮性能预估方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MX2016016942A (es) * 2016-12-13 2018-06-12 Centro De Ingenieria Y Desarrollo Ind Aditamento neumático para limitación de fuerzas aerodinámicas en palas de turbinas eólicas de eje horizontal.

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
基于Wilson 法风力机叶片设计及试验研究;贺玲丽等;《可再生能源》;20130228;第10-14页 *
小型水平轴风力发电机叶片设计研究;杨新彦;《中国优秀博硕士学位论文全文数据库 (硕士)工程科技Ⅱ辑》;20061215;第2-6章 *
杨新彦.小型水平轴风力发电机叶片设计研究.《中国优秀博硕士学位论文全文数据库 (硕士)工程科技Ⅱ辑》.2006,第C042-47页第2-6章. *
水平轴风力机不同叶片设计方法的研究;刘雄飞等;《机械设计与制造》;20170708(第07期);第70-73页 *

Also Published As

Publication number Publication date
CN109325274A (zh) 2019-02-12

Similar Documents

Publication Publication Date Title
AU2020100825A4 (en) Method for evaluating aerodynamic performance of wind wheel of wind turbine based on three-factor fitting integral method
Ingram Wind turbine blade analysis using the blade element momentum method. version 1.1
US6503058B1 (en) Air foil configuration for wind turbine
CN102322407B (zh) 风力机叶片气动外形协同设计方法
Wood Dual purpose design of small wind turbine blades
CN110110427B (zh) 一种大功率风力机叶片的气动外形设计方法
Biadgo et al. Aerodynamic Design of Horizontal Axis Wind Turbine Blades.
CN109325274B (zh) 基于三因子拟合积分法的风力机风轮气动性能评价方法
CN109902384B (zh) 一种基于气弹模型的风力机柔性叶片预弯预扭设计方法
Tenguria et al. Investigation of blade performance of horizontal axis wind turbine based on blade element momentum theory (BEMT) using NACA airfoils
Maniaci An investigation of WT_Perf convergence issues
Chaudhary et al. Modeling and optimal design of small HAWT blades for analyzing the starting torque behavior
CN104405596A (zh) 一种风力发电机组低风速翼型族
CN113742861B (zh) 一种适用于风力发电机风洞试验的叶片模型优化设计方法
Wang et al. Optimized chord and twist angle distributions of wind turbine blade considering Reynolds number effects
CN106777525B (zh) 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法
Duquette et al. Solidity and blade number effects on a fixed pitch, 50 W horizontal axis wind turbine
Kuang et al. Systematic investigation of effect of rotor solidity on vertical-axis wind turbines: Power performance and aerodynamics analysis
Battisti et al. A generalized method to extend airfoil polars over the full range of angles of attack
CN103216381B (zh) 一种风力发电机组叶片
CN101749193A (zh) 可设定启动风速的高效风力发电机及其叶片
Rijs et al. Blade element theory for performance analysis of slow running wind turbines
Sutikno et al. Design and blade optimization of contra rotation double rotor wind turbine
Gertz An evaluation testbed for alternative wind turbine blade tip designs
CN102705173B (zh) 一种风力发电机及其叶片

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