主权项 |
1.风力机叶片气动外形协同设计方法,其特征在于:包括风力机叶片的气动外形设计,所述风力机叶片的气动外形设计包括以下步骤:1)设定参数化翼型型线表达式,设定叶尖速比λ、叶片设计半径R,叶素当地半径r,叶片数B,收敛容许偏差tol,计算精度m;2)设定α(0)和α’(0)的初值α(0)=1/n,α’(0)=0;n为自然数;3)采用下式计算入流角φ:<img file="FDA00002899562300011.GIF" wi="481" he="150" />其中ω为叶片旋转角速度,ρ为空气密度,V<sub>0</sub>为风速;4)采用下式分别计算叶尖损失因子F和F1:<maths num="0001"><![CDATA[<math><mrow><mi>F</mi><mo>=</mo><mfrac><mn>2</mn><mi>π</mi></mfrac><mi>arccos</mi><mo>[</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo><mo>,</mo></mrow></math>]]></maths>其中,<maths num="0002"><![CDATA[<math><mrow><mi>f</mi><mo>=</mo><mfrac><mi>B</mi><mn>2</mn></mfrac><mo>×</mo><mfrac><mrow><mi>R</mi><mo>-</mo><mi>r</mi></mrow><mrow><mi>r</mi><mi>sin</mi><mi>φ</mi></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>R为风轮半径;<maths num="0003"><![CDATA[<math><mrow><msub><mi>F</mi><mn>1</mn></msub><mo>=</mo><mfrac><mn>2</mn><mi>π</mi></mfrac><msup><mi>cos</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>g</mi><mfrac><mrow><mi>B</mi><mrow><mo>(</mo><mi>R</mi><mo>-</mo><mi>r</mi><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mi>r</mi><mi>sin</mi><mi>φ</mi></mrow></mfrac><mo>)</mo></mrow><mo>]</mo><mo>,</mo></mrow></math>]]></maths>其中g=exp[-0.125(BX-21)]+0.1,式中B表示为叶片数;5)采用下式分别计算切向力系数C<sub>x</sub>和法向力系数C<sub>y</sub>:<maths num="0004"><![CDATA[<math><mrow><msub><mi>C</mi><mi>x</mi></msub><mo>=</mo><msub><mi>F</mi><mn>1</mn></msub><msub><mi>C</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>sin</mi><mi>φ</mi><mo>-</mo><mfrac><mn>1</mn><mrow><msub><mi>C</mi><mi>l</mi></msub><mo>/</mo><msub><mi>C</mi><mi>d</mi></msub></mrow></mfrac><mi>cos</mi><mi>φ</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><img file="FDA00002899562300016.GIF" wi="684" he="152" />式中C<sub>l</sub>为翼型升力系数,C<sub>d</sub>为翼型阻力系数,C<sub>x</sub>和C<sub>y</sub>分别为切向力系数和法向力系数;6)采用下式计算风能利用系数C<sub>p</sub>(0):C<sub>p</sub>=[(1-a)<sup>2</sup>+x<sup>2</sup>(1+a′)<sup>2</sup>]xC<sub>x</sub>σ;式中α为轴向诱导因子,α’为周向诱导因子,σ表示叶片实度,计算C<sub>p</sub>(0)时,α和α’分别赋值为α(0)和α’(0);x为叶素所在展长处当地速比,x=X*r/R,其中X为叶尖速比,R为风轮半径;7)进入循环,即在一定的轴向和周向诱导因子时,利用协同设计模型迭代求解,直至C<sub>p</sub>(1)与C<sub>p</sub>(0)的变化小于设定容许偏差tol,完成一次计算;8)重复上述步骤2)~7),在步骤2)中,设定α(n)和α’(n)的初值a(n)=(i/n),a’(n)=(0),n为大于1的整数,通过上述步骤,依次分别求得C<sub>p</sub>(2),C<sub>p</sub>(3)…C<sub>p</sub>(n);i即循环变量,为1,2,3,…n;9)比较C<sub>p</sub>(1),C<sub>p</sub>(2)…C<sub>p</sub>(n)大小,求得最大C<sub>p</sub>以及对应的α,α’;10)计算叶片当地其他几何参数,包括叶片扭角,实度,弦长和厚度,得到翼风动机叶片气动外形的相关参数。 |