用于超长柔性下风向风力机叶片有效扫风半径的计算方法与流程

1.本发明属于风力发电技术领域，特别涉及用于超长柔性下风向风力机叶片有效扫风半径的计算方法。


背景技术：

2.为了使风电与传统化石能源相比具备价格优势，进一步降低其单位发电成本成为风力发电技术研究的目标。为实现这一目标，增大单机容量成为最便捷和最有效的方式，叶片大型化成为风力机发展的必由之路。当前风力发电已经进入10mw级时代，叶片长度达到100m级，超长叶片设计的难点在于叶片气动设计与结构设计之间的矛盾，以及叶片轻质要求与结构性能要求之间的矛盾。传统的上风向叶片必须维持足够的刚度以避免与塔架发生干涉，叶片大型化后质量急剧上升使得在达到临界尺寸之后将不再经济。下风向叶片通过锥角和向后弯曲的设计，离心力距可以抵消一部分推力距，具有被动降载特性，且风载荷将叶片推离塔更远，不存在塔架净空问题，叶片可以更轻更柔，这些优势在超长叶片设计中逐渐突显出来。
3.风载荷作用下，超长柔性下风向叶片会发生大的弯曲变形，且不同风速下其弯曲变形会有很大不同。叶片柔性变形虽然不会打塔，但会影响叶片的扫风面积，进而影响风轮的输出功率。所以在超长柔性下风向叶片设计时，需要考虑柔性变形对叶片有效扫风半径的影响，才能准确预估叶片的设计目标，如年发电量、度电成本等，对于超长柔性下风向风力机叶片的设计与优化至关重要。迄今为止，没有对下风向叶片由于柔性变形计算有效扫风半径的文献及报道。


技术实现要素：

4.针对于上述现有技术的不足，本发明提供了用于超长柔性下风向风力机叶片有效扫风半径的计算方法，该方法用于下风向风力机运行状态下风轮有效扫风半径的计算。所述计算方法包括：采集风力发电机组的运行参数值和必要的超长柔性叶片参数值；根据风力发电机组的气动模型和所采集参数值确定风力作用于风力发电机组的外载参数值；根据风力发电机组的外载参数值、所采集参数值和超长柔性叶片的变形计算模型确定超长柔性叶片的挠度值；根据超长柔性叶片的挠度值、所采集参数值和超长柔性叶片的数学等效模型确定超长柔性下风向风力发电机组风轮有效扫风半径。其中涉及指示风力发电机组外载参数值与所采集参数值关系的风力发电机组气动模型、指示风力发电机组叶片的挠度值与外载参数值关系的超长柔性叶片变形计算模型、指示下风向风力发电机组风轮有效扫风半径与超长柔性叶片的挠度值关系的叶片数学等效模型。使用该方法能较为便捷地计算超长柔性叶片在运行工况下的实际风轮扫风半径，进而计算出风力发电机组的输出功率，在超长柔性下风向风力机叶片的设计及优化方面有很好的应用前景。
5.用于超长柔性下风向风力机叶片有效扫风半径的计算方法，其特征在于，所述计算方法包括：
6.采集风力发电机组的运行参数值和必要的超长柔性叶片参数值；根据风力发电机组的气动模型和所采集参数值确定风力作用于风力发电机组的外载参数值，其中，所述气动模型指示风力发电机组外载参数值与所采集参数值关系；根据风力发电机组外载参数值、所采集参数值和超长柔性叶片的变形计算模型确定超长柔性叶片的挠度值，并将叶片的挠曲线用二阶单项式进行拟合，其中，所述超长柔性叶片变形计算模型指示风力发电机组叶片的挠度值与外载参数值关系；根据超长柔性叶片的挠度值、所采集参数值和超长柔性叶片的数学等效模型确定超长柔性下风向风力发电机组风轮有效扫风半径，其中，所述超长柔性叶片数学等效模型指示超长柔性下风向风力发电机组风轮有效扫风半径与叶片的挠度值关系。
7.进一步的，所述运行参数值包括来流风速、风轮转速以及叶片桨距角，所述必要的超长柔性叶片参数值包括叶片长度、叶片各计算截面距离叶根长度、预弯、刚度等，所述风力发电机组外载参数值为叶片各叶素段气动推力值；其中，所述气动模型指示来流风速、叶尖速比、叶片桨距角以及叶片长度与风力发电机组叶片各叶素段气动推力值之间的对应关系，所述超长柔性叶片变形计算模型指示风力发电机组叶片各叶素段气动推力值、叶片各计算截面距离叶根长度、刚度等与叶片的挠度值之间的对应关系，所述超长柔性叶片数学等效模型指示超长柔性下风向风力发电机组叶片长度、迭代后预弯值与风轮有效扫风半径之间的对应关系；其中，所述叶尖速比根据所述来流风速和所述风轮转速确定。
8.进一步的，所述超长柔性叶片预弯参数为如下分布函数模型：
[0009][0010]
式中，无量纲数为预弯后叶片上任一点的实际展向位置，无量纲数dx为展向处的预弯量，a为叶片预弯段服从分布函数的系数，无量纲数为预弯后叶片叶根线性段的实际展向位置，无量纲数l
f
为预弯后叶片叶尖的实际展向位置。
[0011]
进一步的，所述超长柔性叶片数学等效模型的计算公式为：
[0012][0013][0014]
r'＝r
×
l
f
r0[0015]
式中，无量纲数l
f
为预弯后叶片叶尖的实际展向位置，无量纲数为预弯后叶片上任一点的实际展向位置，无量纲数为预弯后叶片叶根线性段的实际展向位置，a为叶片预弯段服从分布函数的系数，r0为无预弯状态下叶片叶根线性段无量纲位置，r'为风轮有效扫风半径，r为无预弯状态下叶片的长度，r0为风轮轮毂半径。
[0016]
进一步的，所述计算超长柔性下风向风轮有效扫风半径的步骤包括：根据所述运行参数值及叶片参数值、风力发电机组的气动模型、超长柔性叶片的变形计算模型和超长柔性叶片的数学等效模型，通过迭代求解得到所述风轮有效扫风半径。
[0017]
进一步的，所述通过迭代求解得到超长柔性下风向风轮有效扫风半径的步骤包括：
[0018]
步骤1：根据所采集参数值和超长柔性下风向叶片数学等效模型计算风轮有效扫风半径，并更新预弯后的部分叶片展向参数值；
[0019]
步骤2：根据所述运行参数值及更新后的叶片参数值、风力发电机组的气动模型确定风力发电机组的外载参数值；
[0020]
步骤3：根据风力发电机组外载参数值、所采集参数值和叶片的变形计算模型确定超长柔性下风向风力发电机组叶片的挠度值，并将叶片的挠曲线用二阶单项式进行拟合，二阶单项式系数记为a
′
；
[0021]
步骤4：当a
′
与a差值不小于预定阈值时，更新叶片的预弯值a＝a
′
，返回步骤1，重新计算叶片的挠度值，当a
′
与a差值小于预定阈值时，输出当前的超长柔性下风向风轮有效扫风半径，并结束所述迭代求解。
[0022]
进一步的，所述气动模型为基于叶素动量理论、涡方法等针对风力发电机组建立的气动模型；所述超长柔性叶片变形计算模型为基于奇异函数法等针对风力发电机组建立的叶片变形模型。
[0023]
本发明的有益效果为：通过采集较少的叶片参数和简单的叶片所服从分布函数，能有效且准确地构建超长柔性下风向风力发电机组叶片的预弯参数模型；采用该计算方法，能便捷地计算不同风况下的超长柔性下风向风轮有效扫风半径，进而计算实际的风轮输出功率，在超长柔性下风向风力机叶片的设计及优化方面有很好的应用前景。
附图说明
[0024]
图1为本发明用于超长柔性下风向风力机叶片有效扫风半径的计算方法流程图。
[0025]
图2为本发明实施例的超长柔性无预弯叶片与不同预弯叶片的比较图。
[0026]
图3为本发明实施例的超长柔性不同刚度叶片输出功率曲线的比较图。
[0027]
图4为本发明实施例的不同风况下所采集运行参数的示意表格。
[0028]
图5为本发明实施例的不同风况下所计算的超长柔性下风向风轮有效扫风半径结果图。
具体实施方式
[0029]
为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚、明白，以下结合附图和实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
[0030]
以美国可再生能源实验室(nrel)发布的5mw叶片为参考叶片建立超长柔性下风向顺载叶片模型，轮毂半径r0为2m、无预弯状态下的叶片长度r＝66.5m、叶片叶根线性段无量纲位置r0＝0.1，叶片初始预弯服从的二阶单项式系数a＝0.1，如图2所示。不同风况下所采集的运行参数如图3所示。
[0031]
具体计算步骤如下：
[0032]
步骤1：采集以下参数值：来流风速v、风轮转速n以及叶片桨距角β
i
、叶片数b、轮毂半径r0、无预弯叶片长度r、各计算截面距离叶根长度r
i
、弦长c
i
刚度k
i
以及初始预弯模型系数a和叶根线性段无量纲位置r0等。
[0033]
步骤2：根据所采集参数值和超长柔性叶片数学等效模型计算风轮有效扫风半径，
计算公式为：
[0034][0035][0036]
r'＝r
×
l
f
r0[0037]
式中，无量纲数为预弯后叶根线性段的实际展向位置，无量纲数l
f
为预弯后叶片叶尖的实际展向位置，无量纲数为预弯后各计算截面实际展向位置，r'为风轮有效扫风半径，r0为风轮轮毂半径。
[0038]
步骤3：根据无量纲数l
f
更新各计算截面距离叶根长度r
i
：
[0039]
r
i
′
＝r
i
×
l
f
[0040]
步骤4：根据所采集参数值和风力发电机组的气动模型计算风载荷作用于风力发电机组叶片第i个计算截面叶素段dr的气动推力值f
i
：
[0041]
f
i
＝ρbv2c
ti
(λ
i
,β
i
)c
i
dr/2
[0042]
式中，ρ为空气密度，c
ti
为推力系数，叶尖速比λ
i
＝πnr
i
′
/30v，dr＝r
′
i 1
‑
r
i
′
，i＝1,2,3,
…
,。
[0043]
步骤5：根据风力发电机组叶片第i个计算截面叶素段dr的气动推力值f
i
、所采集其他参数值和超长柔性叶片的变形计算模型确定风力发电机组叶片第i个计算截面的挠度值ω
i
：
[0044][0045]
式中，θ
i
为第i个计算截面的转角，m
i
为第i个叶素段固定端面弯矩，f
i
′
为第i个叶素段固定端面反力，i＝1,2,3,
…
,，
[0046]
步骤6：将超长柔性下风向风力发电机组叶片的挠度曲线用二阶单项式函数进行拟合，二阶单项式系数a
′
计算公式为：
[0047][0048]
步骤7：当二阶单项式系数a
′
与初始预弯模型系数a差值不小于预定阈值时，更新叶片的预弯值，a＝a
′
，返回步骤2，重新计算；当二阶单项式系数a
′
与初始预弯模型系数a差值小于预定阈值时，输出超长柔性下风向风轮有效扫风半径r'，并结束计算。
[0049]
上述说明示出并描述了本发明的优选实施例，如前所述，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述发明构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附
权利要求的保护范围内。