• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 12
    • 下载费用:30 金币  

    重庆时时彩开奖结果表: 一种风力机钝尾缘翼型优化设计方法.pdf

    关 键 词:
    一种 风力机 钝尾缘翼型 优化 设计 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201611075093.X

    申请日:

    2016.11.25

    公开号:

    CN106777567A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06F 17/50申请日:20161125|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 天津工业大学
    发明人: 张旭; 王格格; 李伟; 刘海龙; 张孟洁
    地址: 300387 天津市西青区宾水西道399号天津工业大学
    优先权:
    专利代理机构: 代理人:
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201611075093.X

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明公开了一种风力机钝尾缘翼型优化设计方法,包括以下步骤:采用风力机翼型型线集成理论和B样条曲线,形成钝尾缘翼型型线参数化控制方程组;以翼型的形状函数系数、B样条控制参数以及钝尾缘厚度和其在中弧线上侧分配比为设计变量,建立非对称钝尾缘翼型优化设计模型;利用粒子群算法耦合XFOIL软件进行翼型优化设计,提出设计变量计及钝尾缘厚度及其在中弧线上侧分配比的钝尾缘翼型优化设计方法;针对S812翼型优化得到尾缘厚度2.61%弦长、厚度分配比0∶1的钝尾缘改型,采用CFD方法研究原翼型及其改型的升、阻力系数和升阻比。本发明的风力机钝尾缘翼型的气动性能显著提高,更好地提升了风力机的风能利用率,可降低发电成本。

    权利要求书

    1.一种风力机钝尾缘翼型优化设计方法,其特征在于:包括以下步骤:
    步骤(1),钝尾缘翼型型线表达方法:鉴于B样条曲线易于实现局部调控、更逼近特征多
    边形以及能使曲线通过指定点的优势,结合翼型型线集成理论建立钝尾缘翼型型线参数化
    控制方程组;翼型的上翼面距前缘0.4c(c为弦长)之前和下翼面距前缘0.5c之前的型线,均
    采用翼型型线集成理论进行表达,因此该段翼型型线坐标可以表示为:
    <mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>x</mi> <mo>=</mo> <mi>a</mi> <mrow> <mo>(</mo> <mi>&rho;</mi> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> <mo>+</mo> <mn>1</mn> <mo>/</mo> <mi>&rho;</mi> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> <mo>)</mo> </mrow> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&theta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>y</mi> <mo>=</mo> <mi>a</mi> <mrow> <mo>(</mo> <mi>&rho;</mi> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> <mo>-</mo> <mn>1</mn> <mo>/</mo> <mi>&rho;</mi> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> <mo>)</mo> </mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&theta;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
    式中,x为翼型横坐标,y为翼型纵坐标,a为1/4翼型弦长,θ为幅角,ρ(θ)为翼型的形状
    函数,ρ(θ)采用多项式可表达为:
    ρ(θ)=C0+C1θ+C2θ2+…+Ckθk,k=1,2,3,…,n (2)
    式中,C0,C1,C2,…,Ck为形状函数系数,C0=1;通过调整C1,C2,…,Ck,对翼型的形状进行
    控制;
    上翼面距前缘0.4c之后和下翼面距前缘0.5c之后的型线坐标,均采用B样条曲线进行
    表示;为使翼型型线集成理论与B样条曲线在结合点处具有连续、光滑的特性,采用三次B样
    条曲线矩阵形式表示翼型后缘型线坐标,翼型上翼面后缘的型线坐标为:
    <mrow> <msub> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mn>3</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>6</mn> </mfrac> <mo>&lsqb;</mo> <mn>1</mn> <mo>,</mo> <mi>t</mi> <mo>,</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> <mo>,</mo> <msup> <mi>t</mi> <mn>3</mn> </msup> <mo>&rsqb;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>4</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>3</mn> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>3</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>3</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>6</mn> </mrow> </mtd> <mtd> <mn>3</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>3</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>3</mn> </mrow> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mn>6</mn> <mi>P</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <mn>4</mn> <msub> <mi>P</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>6</mn> <mi>P</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>-</mo> <mn>4</mn> <msub> <mi>P</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mi>t</mi> <mo>&Element;</mo> <mo>&lsqb;</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>&rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
    同理,下翼面后缘的型线可以表示为:
    <mrow> <msubsup> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mn>3</mn> </mrow> <mo>&prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>6</mn> </mfrac> <mo>&lsqb;</mo> <mn>1</mn> <mo>,</mo> <mi>t</mi> <mo>,</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> <mo>,</mo> <msup> <mi>t</mi> <mn>3</mn> </msup> <mo>&rsqb;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>4</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>3</mn> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>3</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>3</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>6</mn> </mrow> </mtd> <mtd> <mn>3</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>3</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>3</mn> </mrow> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mn>6</mn> <mi>P</mi> <msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>&prime;</mo> </msup> <mo>-</mo> <mn>4</mn> <msubsup> <mi>P</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <msubsup> <mi>P</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>P</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>P</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>6</mn> <mi>P</mi> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&prime;</mo> </msup> <mo>-</mo> <mn>4</mn> <msubsup> <mi>P</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <msubsup> <mi>P</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mi>t</mi> <mo>&Element;</mo> <mo>&lsqb;</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>&rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
    式中,P0、P1、P2、P3为上翼面后缘型线控制点,P′0、P′1、P′2、P′3为下翼面后缘型线控制
    点;
    步骤(2),钝尾缘翼型优化数学模型:将最大升阻比作为目标函数,选取翼型的形状函
    数ρ(θ)中第2到第12项系数、钝尾缘厚度和其在中弧线上侧分配比以及B样条曲线控制点为
    设计变量,利用粒子群优化算法与XFOIL软件耦合,进行钝尾缘翼型型线的优化设计,钝尾
    缘翼型优化设计流程图为图1;
    步骤(3),钝尾缘翼型优化和气动性能分析:采用上述步骤(2)的钝尾缘翼型优化数学
    模型对翼型进行优化设计,尖尾缘翼型与优化得到钝尾缘翼型的型线如图2所示,利用
    FLUENT软件数值计算其气动性能如图3、图4和图5所示;
    步骤(4),通过步骤(1)至步骤(3)实现风力机钝尾缘翼型优化设计。
    2.根据权利要求1所述的一种风力机钝尾缘翼型优化设计方法,其特征在于:在步骤
    (1)风力机钝尾缘翼型型线表达式建立后,以翼型的形状函数系数、B样条控制参数以及钝
    尾缘厚度和其在中弧线上侧分配比为设计变量,进行非对称钝尾缘翼型优化设计模型的构
    建。
    3.根据权利要求1或2所述的一种风力机钝尾缘翼型优化设计方法,其特征在于:将所
    述翼型采用粒子群算法与XFOIL软件相结合的方法进行钝尾缘翼型优化设计。
    4.根据权利要求3所述的一种风力机钝尾缘翼型优化设计方法,其特征在于:所述
    S812BT翼型的钝尾缘厚度为2.61%c,以及上、下翼面尾缘厚度分配比为0∶1。

    说明书

    一种风力机钝尾缘翼型优化设计方法

    技术领域

    本发明属于翼型改型及优化设计技术领域,尤其是涉及一种利用计算流体力学与
    最优化算法的风力机钝尾缘翼型优化设计方法。

    背景技术

    叶片是风力机捕获风能的重要部件,其翼型气动性能直接影响风力机的风能利用
    率。设计出性能良好的风力机翼型,对提高叶片的风能转换率至关重要。目前,风力机翼型
    的设计方法主要包括反设计和直接优化。直接优化方法能很好解决难以给定恰当目标压力
    和速率分布的问题,并使计算流体动力学(CFD)和最优化方法的多学科交叉设计成为可能。

    国内外学者在风力机翼型直接优化设计研究方面取得了一系列的进展。Ribeiro
    等、Diavareshkian等分别利用Bezier曲线和改进的Hicks-Henne函数参数化表示翼型,采
    用遗传算法耦合CFD以及人工神经网络模型进行翼型优化设计。Liu等修改Bezier曲线基函
    数,建立翼型参数化表达式,采用遗传算法进行优化设计。Sanaye等利用三次Bezier曲线构
    造翼型参数化表示方法,采用遗传算法结合XFOIL与NAFNoise软件对S822翼型进行多目标
    优化设计。Karthikeyan等使用形状分类函数转换法形成翼型参数化表达式,采用Matlab遗
    传算法工具箱调用XFOIL软件进行翼型型线优化设计。白井艳等利用NUMECA软件的
    AUTOBLADE??榻幸硇图负紊杓?,建立风力机专用翼型族。琚亚平等运用Bezier函数建立
    翼型参数化表征方法,采用人工神经网络和遗传算法来优化FX63-167翼型。陈进等基于儒
    可夫斯基保角变换和西奥道生法,提出风力机翼型型线集成理论,并利用遗传算法进行翼
    型优化。该翼型型线集成理论与其他翼型型线表征形式相比,更易于优化和扩展形成新的
    翼型型线,将翼型形状优化问题转化为几何参数优化问题。

    此外,进行钝尾缘改型也是提高叶片的气动性能以及强度、刚度的有效措施。很多
    学者对钝尾缘翼型进行了深入的研究,结果表明钝尾缘改型可以增加最大升力和失速攻
    角,降低最大升力对前缘粗糙的敏感性;并在结构上增大了横截面的面积和弯曲惯性矩,使
    叶片能承担更大的压力载荷。

    鉴于优化设计和钝尾缘改型均能提高翼型的气动性能,且钝尾缘改型同时提高了
    叶片的结构特性,因此进行钝尾缘翼型的优化设计具有重要的意义。陈进等利用翼型型线
    集成理论构造尖尾缘翼型,对尖尾缘翼型的尾缘对称加厚得到钝尾缘翼型,并对其进行优
    化设计,结果显示新翼型具有良好的结构特性和几何兼容性,且气动性能得到提高。上述研
    究在单独考虑翼型型线参数化表达和钝尾缘改型情况下进行,而在钝尾缘翼型型线参数化
    表达中包含钝尾缘厚度及其在中弧线上侧分配比的影响,并进行翼型的优化设计的研究并
    未涉及。

    发明内容

    本发明要解决的问题是提供一种风力机钝尾缘翼型优化设计方法,该方法可采用
    风力机翼型型线集成理论和B样条曲线,形成钝尾缘翼型型线参数化控制方程组;以翼型的
    形状函数系数、B样条控制参数以及钝尾缘厚度和其在中弧线上侧分配比为设计变量,建立
    非对称钝尾缘翼型优化设计模型;利用粒子群算法耦合XFOIL软件进行翼型优化设计,提出
    设计变量计及钝尾缘厚度及其在中弧线上侧分配比的非对称钝尾缘翼型优化设计方法,实
    现风力机钝尾缘翼型优化设计精确性的提高。

    为解决上述技术问题,本发明采用的技术方案是:一种风力机钝尾缘翼型优化设
    计方法,其特征在于:包括以下步骤:

    步骤(1),钝尾缘翼型型线表达方法:鉴于B样条曲线易于实现局部调控、更逼近特
    征多边形以及能使曲线通过指定点的优势,结合翼型型线集成理论建立钝尾缘翼型型线参
    数化控制方程组;翼型的上翼面距前缘0.4c(c为弦长)之前和下翼面距前缘0.5c之前的型
    线,均采用翼型型线集成理论进行表达,因此该段翼型型线坐标可以表示为:


    式中,x为翼型横坐标,y为翼型纵坐标,a为1/4翼型弦长,θ为幅角,ρ(θ)为翼型的
    形状函数,ρ(θ)采用多项式可表达为:

    ρ(θ)=C0+C1θ+C2θ2+...+Ckθk,k=1,2,3,...,n (2)

    式中,C0,C1,C2,...,Ck为形状函数系数,C0=1;通过调整C1,C2,...,Ck,对翼型的形
    状进行控制;

    上翼面距前缘0.4c之后和下翼面距前缘0.5c之后的型线坐标,均采用B样条曲线
    进行表示;为使翼型型线集成理论与B样条曲线在结合点处具有连续、光滑的特性,采用三
    次B样条曲线矩阵形式表示翼型后缘型线坐标,翼型上翼面后缘的型线坐标为:


    同理,下翼面后缘的型线可以表示为:


    式中,P0、P1、P2、P3为上翼面后缘型线控制点,P′0、P′1、P′2、P′3为下翼面后缘型线控
    制点;

    步骤(2),钝尾缘翼型优化数学模型:将最大升阻比作为目标函数,选取翼型的形
    状函数ρ(θ)中第2到第12项系数、钝尾缘厚度和其在中弧线上侧分配比以及B样条曲线控制
    点为设计变量,利用粒子群优化算法与XFOIL软件耦合,进行钝尾缘翼型型线的优化设计,
    钝尾缘翼型优化设计流程图为图1;

    步骤(3),钝尾缘翼型优化和气动性能分析:采用上述步骤(2)的钝尾缘翼型优化
    数学模型对翼型进行优化设计,尖尾缘翼型与优化得到钝尾缘翼型的型线如图2所示,利用
    FLUENT软件数值计算其气动性能如图3、图4和图5所示;

    步骤(4),通过步骤(1)至步骤(3)实现风力机钝尾缘翼型优化设计。

    由于采用上述技术方案,与现有方法相比,本发明的风力机钝尾缘翼型优化设计
    方法,采用风力机翼型型线集成理论和B样条曲线,可形成钝尾缘翼型型线参数化控制方程
    组;以翼型的形状函数系数、B样条控制参数以及钝尾缘厚度和其在中弧线上侧分配比为设
    计变量,可建立非对称钝尾缘翼型优化设计模型;利用粒子群算法耦合XFOIL软件,可实现
    设计变量计及钝尾缘厚度及其在中弧线上侧分配比的非对称钝尾缘翼型优化设计,该优化
    得到的钝尾缘翼型相比于原尖尾缘翼型,其气动性能得到明显改善,从而提高了风力机的
    风能利用率。本发明方法解决了所述的问题,提高了风力机钝尾缘翼型优化设计的精确性,
    为翼型设计提供了技术支持和重要参考。

    附图说明

    下面通过参考附图并结合实例具体地描述本发明,本发明的优点和实现方式将会
    更加明显,其中附图所示内容仅用于对本发明的解释说明,而不构成对本发明的任何意义
    上的限制,在附图中:

    图1是本发明风力机钝尾缘翼型优化设计流程图;

    图2是本发明S812翼型钝尾缘优化前后型线图;

    图3是本发明S812与S812BT翼型的升力系数图;

    图4是本发明S812与S812BT翼型的阻力系数图;

    图5是本发明S812与S812BT翼型的升阻比图。

    具体实施方式

    下面结合实施例及其附图进一步叙述本发明:

    本发明一种风力机钝尾缘翼型优化设计方法基于以下设计思想:

    1、采用风力机翼型型线集成理论和B样条曲线,形成钝尾缘翼型型线参数化控制
    方程组,以防止翼型型线集成理论无法对尾缘形状进行很好的控制;

    2、以翼型的形状函数系数、B样条控制参数以及钝尾缘厚度和其在中弧线上侧分
    配比为设计变量,建立非对称钝尾缘翼型优化设计模型,以防止未计及钝尾缘参数造成翼
    型优化模型的准确度降低;

    3、利用粒子群算法耦合XFOIL软件,实现设计变量计及钝尾缘厚度及其在中弧线
    上侧分配比的非对称钝尾缘翼型优化设计,以防止风力机钝尾缘翼型优化设计的精确性较
    低。

    解决所述技术问题,本发明从钝尾缘翼型型线参数化表达构建和优化数学模型建
    立等方面进行了创新设计:

    1、钝尾缘翼型型线表达方法

    在不改变翼型的最大相对厚度及其位置、弯度和弦长的情况下,利用翼型型线集
    成理论和B样条曲线建立非对称钝尾缘翼型型线参数化控制方程组。

    2、钝尾缘翼型优化数学模型

    风力机钝尾缘翼型优化数学模型的建立必须考虑以下几方面问题:

    (1).提高升力系数的同时,阻力系数变化应较??;

    (2).风力机翼型型线集成理论中形状函数的阶数应能使构成的翼型的几何形状
    与原翼型可以较好地吻合,以利于选取翼型的形状函数ρ(θ)中前几项系数、钝尾缘厚度和
    其在中弧线上侧分配比以及B样条曲线控制点为设计变量;

    (3).控制变量的范围需要进行一定限定,以避免形成的型线不具有翼型的形状特
    征,并可减少多余的迭代次数;

    (4).风力机翼型的相对厚度不宜过大或过小,应使叶片主要功率产生区内翼型的
    最大厚度为0.12~0.25,以避免叶片的结构特性受到影响。

    3、钝尾缘翼型优化和气动性能分析

    利用本发明优化方法,通过Matlab编写优化程序和XFOIL软件计算气动性能,进行
    钝尾缘翼型优化设计,并采用FLUENT软件计算翼型钝尾缘优化前后的升、阻力系数和升阻
    比。

    本发明一种风力机钝尾缘翼型优化设计方法,包括以下步骤:

    步骤(1),钝尾缘翼型型线表达方法:鉴于B样条曲线易于实现局部调控、更逼近特
    征多边形以及能使曲线通过指定点的优势,结合翼型型线集成理论建立钝尾缘翼型型线参
    数化控制方程组;翼型的上翼面距前缘0.4c(c为弦长)之前和下翼面距前缘0.5c之前的型
    线,均采用翼型型线集成理论进行表达,因此该段翼型型线坐标可以表示为:


    式中,x为翼型横坐标,y为翼型纵坐标,a为1/4翼型弦长,θ为幅角,ρ(θ)为翼型的
    形状函数,ρ(θ)采用多项式可表达为:

    ρ(θ)=C0+C1θ+C2θ2+...+Ckθk,k=1,2,3,...,n (2)

    式中,C0,C1,C2,...,Ck为形状函数系数,C0=1;通过调整C1,C2,...,Ck,对翼型的形
    状进行控制;

    上翼面距前缘0.4c之后和下翼面距前缘0.5c之后的型线坐标,均采用B样条曲线
    进行表示,即:


    式中,p(u)为B样条曲线坐标,di为控制点,Ni,k(u)为k次基函数,u为结点矢量;

    为使翼型型线集成理论与B样条曲线在结合点处具有连续、光滑的特性,采用三次
    B样条曲线矩阵形式表示翼型后缘型线坐标,翼型上翼面后缘的型线坐标为:


    翼型下翼面后缘的型线坐标为:


    式中,P0、P1、P2、P3为上翼面后缘型线控制点,P′0、P′1、P′2、P′3为下翼面后缘型线控
    制点;

    上翼面后缘型线P0,3(t)的点P0,3(0)经过翼型型线集成理论形成的上翼面型线中
    最后一个坐标点,P0,3(1)经过钝尾缘翼型上翼面型线终点(1,h×k),其中h为钝尾缘厚度,k
    为上翼面尾缘厚度与钝尾缘厚度的比值;通过公式(4)即可反算出p0和p3,则上翼面后缘型
    线的控制变量只有p1和p2,因此翼型上翼面后缘的型线具体表示为:


    同理,下翼面后缘的型线可以表示为:


    公式(1)、(6)和(7)即为本发明利用翼型型线集成理论和B样条曲线建立的非对称
    钝尾缘翼型型线参数化控制方程组;

    步骤(2),钝尾缘翼型优化数学模型:以翼型的形状函数系数、B样条控制参数以及
    钝尾缘厚度和其在中弧线上侧分配比为设计变量,建立非对称钝尾缘翼型优化设计模型;

    升、阻力系数是研究风力机翼型气动性能的两个重要指标,提高翼型的气动性能
    旨在使升力系数提高的同时,阻力系数变化较小,因此将最大升阻比作为目标函数:

    f(x)=max(CL/CD) (8)

    式中,CL为翼型的升力系数,CD为翼型的阻力系数;

    利用风力机翼型型线集成理论对翼型进行参数化拟合,当形状函数阶数为11阶
    时,构成的翼型的几何形状与原翼型可以较好地吻合,因此选取翼型的形状函数ρ(θ)中第2
    到第12项系数、钝尾缘厚度和其在中弧线上侧分配比以及B样条曲线控制点为设计变量,
    即:

    X=(C1,C2,C3,C4,C5,C6,C7,C8,C9,C10,C11,h,k,P1,P2,P1′,P2′) (9)

    在翼型优化过程中,需要对控制变量的范围进行一定限定,以避免形成的型线不
    具有翼型的形状特征,同时可减少多余的迭代次数,则优化变量边界约束条件为:

    Xmin≤X≤Xmax (10)

    风力机叶片的主要功率产生区为沿翼展方向约75%处,且翼型的相对厚度对叶片
    结构特性至关重要;大部分叶片在该区域采用翼型的最大厚度一般为0.12~0.25,因此翼
    型的相对厚度约束条件为:

    0.12≤t≤0.25 (11)

    粒子群优化算法具有实现容易、精度高、收敛快等优点,解决实际问题的效率很
    高,将该算法与XFOIL软件耦合,进行钝尾缘翼型型线的优化设计,钝尾缘翼型优化设计流
    程图为图1。

    步骤(3),钝尾缘翼型优化和气动性能分析:基于上述钝尾缘翼型优化数学模型,
    利用粒子群算法耦合XFOIL软件,实现设计变量计及钝尾缘厚度及其在中弧线上侧分配比
    的非对称钝尾缘翼型优化设计,采用FLUENT软件计算翼型钝尾缘优化前后的升、阻力系数
    和升阻比如图3、图4和图5所示;

    步骤(4),通过步骤(1)至步骤(3)实现风力机钝尾缘翼型优化设计。

    本发明未述及之处适用于现有技术。

    实施例:

    1、S812翼型广泛应用于风力机叶片的主要功率产生区,具有位于39%c处的最大
    相对厚度21%,以及位于77.1%c处的最大相对弯度1.57%。

    2、优化初始条件为:雷诺数Re=1×106,马赫数Ma=0.04。粒子群算法参数设定:
    学习因子S1、S2均取0.5,变量维数为21,种群大小取20;惯性权重w决定了粒子之前的速度对
    现在速度的影响程度,采用Y.shi和R.Eberhart提出的w线性减少自适应调整公式,使算法
    程序可以在初期具有较强的搜索能力,后期还具有较好的收敛性,该公式为:


    式中,wmax和wmin分别表示惯性权重的最大值和最小值,根据设计经验选取1.2和
    0.4;t和tmax分别表示当前迭代次数和最大迭代次数,最大迭代次数为300。

    3、钝尾缘翼型优化和气动性能分析

    采用本发明优化方法,通过Matlab编写优化程序和XFOIL软件计算气动性能,进行
    美国国家可再生能源实验室(NREL)设计的中等厚度风力机专用翼型S812的钝尾缘优化设
    计,并利用FLUENT软件进行气动性能计算。

    本发明所述风力机钝尾缘翼型优化设计的S812翼型钝尾缘改型S812BT的形状如
    图2所示。对其升、阻力系数和升阻比进行研究分析表明:(1)S812BT翼型的尾缘厚度为
    2.61%c,以及上、下翼面尾缘厚度分配比为0∶1。(2)S812BT翼型的升力系数相较于S812翼
    型明显升高;阻力系数在12.1°攻角之前与S812翼型非常接近,超过12.1°攻角之后高于
    S812翼型。S812翼型的失速攻角为16.2°,S812BT翼型在20.2°攻角时还未失速,说明优化后
    的钝尾缘翼型的失速推迟了。(3)S812BT翼型的升阻比在攻角小于15.3°时明显高于S812翼
    型;攻角大于15.3°之后,与S812翼型十分接近。S812BT翼型的最大升阻比较S812翼型有大
    幅提高。

    关于本文
    本文标题:一种风力机钝尾缘翼型优化设计方法.pdf
    链接地址://www.4mum.com.cn/p-6021130.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 重庆时时彩稳赚杀号 澳门威尼人网站赌博 时时彩组三稳赚方法技巧 一晚上搞10万的门路 pk10技巧规律 幸运飞艇助力计划 大乐透走势图表近50期 双色球稳赚不赔的方法 开奖特马白小姐 黑龙江时时玩法 pk10有什么办法稳赢 幸运飞艇免费计划安卓软件 大乐透走势图表图 时时彩大小历史数据 时时彩计划软件手机版 北京pk10免费全天计划