发明名称 风力机叶片气动外形协同设计方法
摘要 本发明公开了一种风力机叶片气动外形协同设计方法,包括对风力机叶片的气动外形设计,以及利用风力机叶片的数学模型,针对沿叶片展向不同位置风力机专用翼型的设计要求,以二维风能利用系数为优化目标的风力机翼型设计;其中叶片设计方法有别于现有叶片设计方法,不考虑叶片设计过程中轴向和周向诱导因子的求解方法,通过对<img file="DEST_PATH_IMAGE002.GIF" wi="14" he="16" />和<img file="224869DEST_PATH_IMAGE002.GIF" wi="14" he="16" />’直接赋值,进而求解对应的最大风能利用系数,计算方法直接简便,同时在计算过程中将翼型的设计与优化考虑进来,实现了叶片与翼型设计的共用模型,为风力机叶片与翼型的协同设计奠定了基础,整套方法构思巧妙,易于操作,为风力机叶片和翼型的优化设计提供了新的思路。
申请公布号 CN102322407B 申请公布日期 2013.06.12
申请号 CN201110252679.X 申请日期 2011.08.30
申请人 重庆大学 发明人 陈进;程江涛;汪泉;庞晓平
分类号 F03D11/00(2006.01)I;G06F17/50(2006.01)I 主分类号 F03D11/00(2006.01)I
代理机构 北京同恒源知识产权代理有限公司 11275 代理人 赵荣之
主权项 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>&pi;</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>&times;</mo><mfrac><mrow><mi>R</mi><mo>-</mo><mi>r</mi></mrow><mrow><mi>r</mi><mi>sin</mi><mi>&phi;</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>&pi;</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>&phi;</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>&phi;</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>&phi;</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)计算叶片当地其他几何参数,包括叶片扭角,实度,弦长和厚度,得到翼风动机叶片气动外形的相关参数。
地址 400044 重庆市沙坪坝区沙坪坝正街174号