失速型风力机翼型的设计方法与流程

本发明涉及水平轴风力机叶片大厚度翼型的设计领域，具体涉及一种失速型风力机翼型的设计方法。
背景技术：
：目前针对风力机叶片较大厚度翼型的研究，重点都是研究失速之前的气动性能，很少有学者研究失速之后翼型气动性能的提高。然而，由于水平轴风力机叶片大厚度翼型(例如最大相对厚度30％、40％等)通常运行在失速工况下。因此，提高失速工况下风力机叶片气动性能显得尤为重要。一般风力机翼型都是在一定设计攻角情况下设计出来的，失速之前翼型气动性能较好，但是失速之后翼型气动性能下降得较为厉害。现有技术中并没有采用一种独特的翼型吸力面廓线局部表达方法，进而改变翼型失速之后的流场分布。因此，需要改进现有技术中的翼型吸力面廓线，从而提高风力机叶片大厚度翼型失速之后的气动性能。技术实现要素：本发明的目的在于提供一种失速型风力机翼型的设计方法，通过本方法设计出来的翼型能够显著提高风力机叶片大厚度翼型失速之后的气动性能。本发明解决上述技术问题所采用的方案是：一种失速型风力机翼型的设计方法，其特征在于，包括：s1：在原始翼型廓线的基础上，采用在吸力面局部使用椭圆函数生成内凹形状的廓线进行表征，并确定设计变量以及约束条件；s2：综合考虑失速之前及失速之后翼型的气动性能，给定失速之前和失速之后的设计功角以及设计雷诺数，建立在光滑条件下翼型升阻比最大作为目标函数；s3：以步骤s2中的目标函数为优化目标，并通过约束条件过滤不符合要求的廓线，给定优化算法对符合要求的延迟失速翼型廓线进行优化计算得到目标翼型廓线。进一步地，步骤s1中的椭圆形函数的表达式为：上式中，系数a>b>0，x为翼型在xoy平面横坐标，y为翼型xoy平面纵坐标，c,d为平面坐标偏移参数，θ为变换角度。进一步地，步骤s1中以a、b、c、d及θ五个参数作为设计变量：x＝(a,b,c,d,θ)；约束条件为：将翼型局部椭圆控制函数两个系数进行如下约束：xmin≤x≤xmax设计变量约束范围为：abcdθ/°最大值0.00.00.00.0-30.0最小值1.01.01.01.00.0。进一步地，步骤s2中的目标函数的表达式为：f(x)＝max(f1(x),f2(x))其中，f1(x)＝cl/cd，f2(x)＝cl/c'd式中，cl为失速之前设计攻角在一定角度时的升力系数，cd为失速之前设计攻角在一定角度时的阻力系数；c'l为失速之后设计攻角在一定角度时的升力系数，cd'为失速之后设计攻角在一定角度时的阻力系数。进一步地，步骤s2中，给定失速之前设计攻角为6°，失速之后的设计攻角为14°，设计雷诺数为re＝3.0×106。进一步地，步骤s3还包括：(a)初始化设计变量，生成采用椭圆函数产生内凹形状的翼型廓线，通过约束条件过滤掉非延迟失速翼型廓线；(b)调入符合要求的延迟失速翼型廓线，获取延迟失速翼型廓线上的数据点坐标信息；(c)设置需要计算的翼型气动性能参数，对上述符合要求的一系列延迟失速翼型族，计算收敛后输出每个翼型廓线的数据点信息及其对应的气动性能参数(d)对上述输出的翼型族及其对应的气动性能参数，通过目标函数计算适应度值；(e)将上述四个步骤的计算模块耦合到模拟退火多目标进化算法中进行优化迭代，当满足迭代条件，则输出目标翼型廓线；若不满足迭代条件，则对翼型廓线进行自适应调整，回转执行步骤(a)。进一步地，上述步骤(c)中气动性能参数包括设计攻角在一定角度范围内的升力系数以及阻力系数。进一步地，模拟退火算法相关参数如下：种群数量为30，最大迭代数200，降温速率0.9，初始温度1000，结束温度0.001，链长220。进一步地，将原始翼型吸力面在靠近50％弦长处用半椭圆形状生成内凹形状。本发明的另一个目的是提供一种根据上述的失速型风力机大厚度翼型的设计方法优化设计得出的延迟失速型风力机叶片翼型。与现有技术相比，本发明至少具有以下有益效果：本发明方法针对水平轴风力发电机叶片翼型，提出在翼型吸力面局部采用椭圆内凹廓线的翼型设计方法，以失速之前最大升阻比及失速之后最大升阻比为多目标函数，在du97-w-300翼型基础上，翼型吸力面局部采用椭圆函数表达，并利用多目标模拟退火算法进行翼型局部廓线优化设计，从而使得设计出来的新翼型具有明显的延迟失速特性，其后失速气动性能明显提高，相比传统的水平轴风力发电机叶片翼型而言，该新翼型具有更优越的后失速气动性能，使得叶片能够在失速状态下也能产生较多的风能；本发明方法可以推广应用到水平轴风力发电机叶片翼型上面，采用该新翼型替换传统的水平轴风力发电机叶片翼型，具有良好的社会价值和经济效益。附图说明图1为本发明实施例翼型优化设计的流程图；图2为本发明实施例翼型吸力面局部半椭圆形状的示意图；图3为本发明实施例icem优化延迟失速翼型网格划分的示意图；图4为本发明实施例延迟失速wt-e-300新翼型廓线示意图；图5为本发明实施例wt-e-300翼型升力系数；图6为本发明实施例wt-e-300翼型阻力系数图7为本发明实施例wt-e-300翼型升阻比。具体实施方式为更好的理解本发明，下面的实施例是对本发明的进一步说明，但本发明的内容不仅仅局限于下面的实施例。鉴于现有技术中一般风力机翼型都是在一定设计攻角情况下设计出来的，失速之前翼型气动性能较好，但是失速之后翼型气动性能下降得较为厉害，本发明提供一种失速型风力机大厚度翼型的设计方法，本发明设计出来的翼型能够显著提高风力机叶片大厚度翼型失速之后的气动性能。本发明在通用风力机专用翼型du97-w-300(该翼型最大相对厚度为30％)的基础上，提出以设计工况条件下及失速工况条件下翼型升阻比最大为多目标函数的风力发电机叶片翼型优化设计方法，采用模拟退火多目标进化算法，实现水平轴风力发电机翼型吸力面局部廓线的控制与参数优化，并将翼型流体计算软件fluent耦合到该算法程序中，计算翼型的气动性能。最后，将优化设计出来的wt-e-300新翼型与风力发电机叶片常用的du97-w-300翼型进行了气动性能对比分析。具体地，本实施例的设计流程如图1所示，包括如下步骤：s1：在原始翼型廓线的基础上，采用在吸力面局部使用椭圆函数生成内凹形状的廓线进行表征，并确定设计变量以及约束条件；见图2，本发明在风力机叶片常用翼型du97-w-300廓线1的基础上，将翼型吸力面靠近50％弦长处用半椭圆2内凹形状代替原廓线上该处的形状得到新的翼型廓线形状，其吸力面局部廓线椭圆表达式(1)所示：式(1)中，a>b>0；x为翼型在xoy平面横坐标，y为翼型xoy平面纵坐标；选取适当的系数a，b取值，并通过平面坐标偏移c,d参数与旋转变换θ角度，可得如图2所示得翼型形状，该翼型廓线吸力面具有局部椭圆内凹型特性，能够有效提高翼型后失速气动性能。确定设计变量根据翼型吸力面廓线局部半圆椭圆函数表达的思想，通过控制椭圆函数a、b的系数便可得到叶片翼型吸力面局部廓线内凹形状，本发明选取翼型吸力面局部椭圆廓线控制函数a、b、c、d及θ五个参数作为设计变量：x＝(a,b,c,d,θ)(2)确定约束条件为了使水平轴风力发电机叶片翼型吸力面局部椭圆廓线在可控制的范围内变化，将翼型局部椭圆控制函数两个系数进行如下约束：xmin≤x≤xmax(3)设计变量约束范围如表1所示：表1设计变量范围abcdθ/°最大值0.00.00.00.0-30.0最小值1.01.01.01.00.0。s2：综合考虑失速之前及失速之后翼型的气动性能，给定失速之前和失速之后的设计功角以及设计雷诺数，建立在光滑条件下翼型升阻比最大作为目标函数；综合考虑失速之前及失速后翼型的气动性能，本发明提出在翼型失速之前设计攻角为6°与失速后翼型设计攻角为14°工况下，在设计雷诺数为re＝3.0×106，马赫数ma＝0.15的条件下，以光滑条件下最大升阻比作为多目标函数：f(x)＝max(f1(x),f2(x))(4)f1(x)＝cl/cd(5)f2(x)＝c’l/c'd(6)式(5)中，cl为失速之前设计攻角在6°时的升力系数，cd为失速之前设计攻角在6°时的阻力系数；c'l为失速后设计攻角在14°时升力系数，c'd为失速后设计攻角在14°时的阻力系数。s3：以步骤s2中的目标函数为优化目标，并通过约束条件过滤不符合要求的廓线，给定优化算法对符合要求的延迟失速翼型廓线进行优化计算得到目标翼型廓线；为实现延迟失速翼型的集成耦合优化设计，编制四个模块化程序：延迟失速翼型廓线生成、icem自适应网格技术(如图3所示)、fluent数值计算及优化设计模块。通过编制数据接口程序将这四个模块耦合到模拟退火算法中进行优化迭代。具体的包括如下步骤：(a)初始化变量，生成采用椭圆函数产生内凹形状的翼型廓线，并通过约束条件过滤掉非延迟失速翼型廓线；(b)编制数据传递程序将符合要求的延迟失速翼型廓线数据文件调入到icem自适应网格划分模块中，icem自适应网格划分模块读入数据文件，并对延迟失速型翼型廓线自动生成结构化网格从而读取廓线上的数据点坐标信息，且将该数据点坐标信息输出fluent可读的网格格式；(c)fluent计算模块调用步骤(b)输出的网格信息，设置需要计算的翼型气动性能参数，该气动性能参数包括设计攻角在一定角度范围内的升力系数以及阻力系数，对上述符合要求的一系列延迟失速翼型族，计算收敛后输出每个翼型廓线上的数据点信息及每个翼型对应的气动性能参数；(d)对上述输出的翼型族数据及其对应的气动性能参数，通过目标函数计算适应度值；(e)将上述(a)、(b)、(c)、(d)四个子程序模块耦合到模拟退火算法优化程序中进行优化迭代，从而进行线路规划、升温、恒温及降温等操作，模拟退火算法相关参数如下：种群数量为30，最大迭代数200，降温速率0.9，初始温度1000，结束温度0.001，链长220；当满足迭代条件，则输出目标翼型廓线；若不满足迭代条件，则对翼型廓线进行自适应调整，回转执行步骤(a)通过模拟退火算法优化迭代，最终输出最大相对厚度为30％的延迟失速新翼型，取名wt-e-300，如图4所示。为了验证该新翼型具有较高的气动性能，尤其是失速之后的气动性能。图5、图6和图7给出了新翼型wt-e-300与传统的du97-w-300翼型的气动性能对比图。表2给出了wt-e-300翼型、du97-w-300翼型的最大升力系数、最大升阻比及一定攻角范围内的平均气动数据。由图表可知：当攻角大于14°时，翼型处于失速状态，wt-e-300新翼型的气动性能优于du97-w-300翼型，明显表现出延迟失速特性。新翼型wt-e-300的最大升力系数为1.581，相比du97-w-300翼型增加了约6.0％；wt-e-300翼型的平均升力系数为1.262，而du97-w-300翼型的平均升力系数为1.126，增加了12.1％；虽然wt-e-300新翼型的最大升阻比略有降低，但是其平均升阻比为37.128，而du97-w-300翼型的平均升阻比为35.069，增加了5.9％。表2翼型气动性能参数对比注：括号内表示攻角位置或者范围，cl,max为最大升力系数，cl,aver为一定攻角范围内平均升力系数，l/d,max为最大升阻比，l/d,aver为一定攻角范围内平均升阻比。以上所述是本发明的优选实施方式而已，当然不能以此来限定本发明之权利范围，应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和变动，这些改进和变动也视为本发明的保护范围。当前第1页12