一种低风速变速变桨风力机叶片优化设计方法与流程

文档序号:12720869阅读:来源:国知局

技术特征:

1.一种低风速变速变桨风力机叶片优化设计方法,其特征在于,包括以下步骤:

1)设定优化目标,目标一为风力机年发电量最大,目标二为叶片材料成本最小;

2)设定叶片优化设计的变量;

3)设定叶片外形约束条件;

4)同时在保证目标一和目标二的前提下,通过寻优算法计算叶片每个截面的弦长和扭角的最优解。

2.根据权利要求1所述的一种低风速变速变桨风力机叶片优化设计方法,其特征在于,所述方法还包括对最佳桨距角进行最优控制。

3.根据权利要求1所述的一种低风速变速变桨风力机叶片优化设计方法,其特征在于,所述方法还包括对转矩-转速进行最优控制。

4.根据权利要求1所述的一种低风速变速变桨风力机叶片优化设计方法,其特征在于,所述步骤1)将目的一和目的二优化为求解以下目的函数的最小值,

<mrow> <mi>F</mi> <mo>=</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mo>{</mo> <mi>&mu;</mi> <mfrac> <mrow> <msub> <mi>X</mi> <mn>1</mn> </msub> <mo>-</mo> <mi>A</mi> <mi>E</mi> <mi>P</mi> </mrow> <mrow> <msub> <mi>X</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> </mrow> </mfrac> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&mu;</mi> <mo>)</mo> </mrow> <mfrac> <mrow> <mi>A</mi> <mi>R</mi> <mi>E</mi> <mi>A</mi> <mo>-</mo> <msub> <mi>X</mi> <mn>4</mn> </msub> </mrow> <mrow> <msub> <mi>X</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>X</mi> <mn>4</mn> </msub> </mrow> </mfrac> <mo>}</mo> </mrow>

其中,μ为权重系数μ∈[0,1];常值X1、X2、X3及X4,分别为约束范围内出现的最大年发电量、最小年发电量、最大叶片面积和最小叶片面积。

5.根据权利要求2所述的一种低风速变速变桨风力机叶片优化设计方法,其特征在于,所述步骤2)具体为叶片每个截面的弦长分布及扭角分布均以贝塞尔曲线表示;所述叶片优化设计的变量包括弦长控制点CPi(i=4,5...8)和扭角控制点TPi(i=1,2...4)及最佳桨距角θopt共10个变量。

6.根据权利要求5所述的一种低风速变速变桨风力机叶片优化设计方法,其特征在于,所述10个变量通过以下方程进行约束,

CPmax≥CP4≥CP5≥CP6≥CP7≥CP8≥CPmin

TPmax≥TP1≥TP2≥TP3≥TP4≥TPmin

θmax≥θopt≥θmin

CPmin、CPmax、TPmin、TPmax、θmin、θmax均为实常数。

7.根据权利要求2所述的一种低风速变速变桨风力机叶片优化设计方法,其特征在于,通过以下方程对最佳桨距角进行最优控制

A(i)=C1A1(i)+C2A2(i)+C3A3(i)

<mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>C</mi> <mrow> <mi>P</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>-</mo> <mi>min</mi> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mrow> <mi>P</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>max</mi> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mrow> <mi>P</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>min</mi> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mrow> <mi>P</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>

<mrow> <msub> <mi>A</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>max</mi> <mrow> <mo>(</mo> <msub> <mi>B</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>B</mi> <mi>i</mi> </msub> </mrow> <mrow> <mi>max</mi> <mrow> <mo>(</mo> <msub> <mi>B</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>min</mi> <mrow> <mo>(</mo> <msub> <mi>B</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <msub> <mi>B</mi> <mi>i</mi> </msub> <mo>=</mo> <msqrt> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>&lambda;</mi> <mo>=</mo> <msub> <mi>&lambda;</mi> <mn>0</mn> </msub> </mrow> <msub> <mi>&lambda;</mi> <mn>1</mn> </msub> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mrow> <mi>P</mi> <mi>i</mi> </mrow> </msub> <mo>(</mo> <mi>&lambda;</mi> <mo>)</mo> <mo>-</mo> <msub> <mi>Cp</mi> <mrow> <mi>max</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow>

<mrow> <msub> <mi>A</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>max</mi> <mrow> <mo>(</mo> <msub> <mi>&lambda;</mi> <mrow> <mi>o</mi> <mi>p</mi> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&lambda;</mi> <mrow> <mi>o</mi> <mi>p</mi> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> </mrow> <mrow> <mi>max</mi> <mrow> <mo>(</mo> <msub> <mi>&lambda;</mi> <mrow> <mi>o</mi> <mi>p</mi> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>min</mi> <mrow> <mo>(</mo> <msub> <mi>&lambda;</mi> <mrow> <mi>o</mi> <mi>p</mi> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>

其中,C1、C2、C3为权重系数,A1(i)为CPmax评价指标,A1(i)为CP曲线顶端平滑评价指标,A1(i)为叶尖速比λ评价指标,CPmax,i为第i个桨距角下CP曲线中的最大CP值;Bi表示第i个桨距角下CPmax,i附近不同叶尖速比λ下CPi(λ)与CPmax,i的标准差;λopt,i为第i个桨距角下CPmax对应的叶尖速比。

8.根据权利要求3所述的一种低风速变速变桨风力机叶片优化设计方法,其特征在于,通过以下方程对转矩-转速进行最优控制,

<mrow> <msub> <mi>&omega;</mi> <mi>g</mi> </msub> <mo>=</mo> <msqrt> <mfrac> <mi>M</mi> <msub> <mi>K</mi> <mrow> <mi>o</mi> <mi>p</mi> <mi>t</mi> </mrow> </msub> </mfrac> </msqrt> </mrow>

<mrow> <msub> <mi>K</mi> <mrow> <mi>o</mi> <mi>p</mi> <mi>t</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>&pi;&rho;R</mi> <mn>5</mn> </msup> <msub> <mi>C</mi> <mrow> <mi>P</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mfrac> <mn>1</mn> <mrow> <msup> <mi>G</mi> <mn>3</mn> </msup> <msup> <msub> <mi>&lambda;</mi> <mrow> <mi>o</mi> <mi>p</mi> <mi>t</mi> </mrow> </msub> <mn>3</mn> </msup> </mrow> </mfrac> </mrow>

其中,M为转矩观测值,ρ为空气密度;R为风轮半径;CPmax为选定桨距角下最大功率系数;ωg为发电机转速;G为齿轮箱传动比;λopt为CPmax对应的叶尖速比;Kopt为最优转矩-转速控制系数。

9.根据权利要求1所述的一种低风速变速变桨风力机叶片优化设计方法,其特征在于,所述步骤4)采用MATLAB编制多岛遗传算法来进行优化,叶片气动性能的计算采用GHBLADED,其计算模型采用叶素动量理论来对气动性能进行计算。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1