一种变速风机叶片的多工况气动优化的改进方法与流程

本发明属于风力机叶片气动优化设计领域，特别是一种变速风机叶片的多工况气动优化的改进方法。
背景技术：
：近年来，变速风机(vswt)已成为大型风电机组的主流机型。在低于额定风速时，风机普遍采用最大功率点跟踪(maximumpowerpointtracking，mppt)控制策略。它会根据风速的变化调整风轮转速，使风轮运行在最佳叶尖速比λopt，并以最大风能利用系数cp,max捕获风能。因此，传统的针对变速风机风轮的气动设计一般仅优化风轮在单一工况点(即λopt)的气动性能。而实际上，由于风轮转动惯量较大，风力机无法快速响应风速的突然变化、从而常运行在跟踪λopt的过程中而非λopt处。考虑到变速风机无法保持运行在单一叶尖速比，为提高其在多工况点的平均气动性能，现有技术已经提出以多个叶尖速处的风能利用系数cp的加权和作为目标函数的多工况气动优化设计方法。该方法通过风力机闭环系统(含风轮、传动链和mppt控制)的数值仿真获得对应于不同运行叶尖速比区间的来流风能、并将其作为量化描述各工况重要程度的指标，然后基于重点提升对应来流风能较大的叶尖速比处的cp的思路，确定目标函数中的多个设计叶尖速比及其对应权值。仿真结果表明，该方法能够提高变速风机在多个运行叶尖速比λope处的平均气动性能。但是，由于该多工况方法是根据初始风轮的λope的分布特性来确定目标函数进而获得优化风轮，为保证该方法有效，必须保证优化风轮与初始风轮的λope分布相同，且优化风轮能够在这些固定的叶尖速比处(即对应初始风轮来流风能较大的叶尖速比)提升cp。而实际上，风力机在不同叶尖速比所能达到的功率系数极限是不同的，也即风力机在不同叶尖速比提升cp的难易程度是不同的。因此，主动调整运行工况使其集中分布于更容易提升cp的叶尖速比处，可以进一步提升风力机的平均气动性能。由上可知，现有技术固定了优化风轮的运行叶尖速比λope分布，使得运行工况单一且固定，不能主动调整至提升风能利用系数cp更容易的工况下，因而风力机的平均气动性能提升比例有待提高。技术实现要素：本发明所解决的技术问题在于提供一种变速风机叶片的多工况气动优化的改进方法。实现本发明目的的技术解决方案为：一种变速风机叶片的多工况气动优化的改进方法，包括以下步骤：步骤1、对控制点进行初始化，利用贝塞尔曲线拟合初始叶片的弦长和扭角曲线，得到控制点的初始值；同时构建初始表格，该表格为最优叶尖速比-设计叶尖速比-权重表，表格包含上述三种参数；步骤2、利用贝塞尔曲线的控制点生成叶片气动外形，该外形包括弦长和扭角分布；步骤3、判断叶片外形是否满足所有约束条件，若满足，执行步骤4，否则，将目标函数值设为零，执行步骤7；步骤4、调用bladed软件计算风能利用系数和叶尖速比的关系曲线以及最优叶尖速比λopt；所述风能利用系数和叶尖速比的关系曲线为cp-λ曲线；步骤5、根据最优叶尖速比λopt更新目标函数，确定需要考虑的工况以及对应工况的权重系数μj；步骤6、计算目标函数值，得到平均风能捕获效率pfavg；步骤7、判断是否满足终止条件，若不满足，通过调用遗传算法产生新一代个体，转向步骤2；若满足，则输出最优叶片的气动外形。本发明与现有技术相比，其显著优点为：1)本发明考虑了风轮实际运行中受mppt动态过程影响，在设计优化阶段，提升多个工况下的平均气动性能；2)本发明提出了在优化过程中主动调整运行工况的方法，主动调整运行工况使其集中分布于更容易提升cp的叶尖速比处，可以进一步提升风力机的平均气动性能；3)本发明无需限制优化风轮的λope分布特性与初始风轮，使得所得优化风轮的运行工况集中分布于更容易提升风能利用系数cp的叶尖速比处，能够进一步提高风力机的平均风能捕获效率。下面结合附图对本发明作进一步详细描述。附图说明图1为本发明的改进的多工况气动设计方法流程图。图2为本发明优化后的风轮与其他两种优化方法获得的风轮以及初始风轮的气动外形对比图，其中图(a)为弦长对比图，图(b)为扭角对比图。其中，通过集总参数λopt的多工况改进设计方法标记为mpvwblade，固定优化风轮λope分布的多工况设计方法标记为mpcwblade，多工况数值计算的方法标记为mpncblade，以及初始叶片标记为originalblade。图3为本发明优化后的风轮与其他两种优化方法得到的风轮分别与初始风轮进行来流风能比例分布和风能利用系数的对比图，其中图(a)为mpvwblade与originalblade的对比图，图(b)为mpcwblade与originalblade的对比图，图(c)为mpncblade与originalblade的对比图。具体实施方式结合附图，本发明的一种变速风机叶片的多工况气动优化的改进方法，包括以下步骤：步骤1、对控制点进行初始化，利用贝塞尔曲线拟合初始叶片的弦长和扭角曲线，得到控制点的初始值；同时构建初始表格，该表格为最优叶尖速比-设计叶尖速比-权重表，表格包含上述三种参数；步骤2、利用贝塞尔曲线的控制点生成叶片气动外形，该外形包括弦长和扭角分布；步骤3、判断叶片外形是否满足所有约束条件，若满足，执行步骤4，否则，将目标函数值设为零，执行步骤7；所述约束条件为：1)弦长变化不超过初始叶片弦长的1.05倍；2)弦长变化不低于初始叶片弦长的0.8倍；3)叶片弦长和扭角曲线平滑过渡，从最大值位置开始直到叶尖位置的弦长和扭角均单调递减。步骤4、调用bladed软件计算风能利用系数和叶尖速比的关系曲线以及最优叶尖速比λopt；所述风能利用系数和叶尖速比的关系曲线为cp-λ曲线；步骤5、根据最优叶尖速比λopt更新目标函数，确定需要考虑的工况以及对应工况的权重系数μj；所述目标函数为：式中，k是设计叶尖速比λj的个数，是在λj处的风能利用系数，μj是在处的权重系数，并且满足0＜μj＜1和是在λj处的来流风能比例，pfavg是平均风能捕获效率；根据最优叶尖速比λopt更新目标函数包括以下步骤：步骤5-1、将最优叶尖速比λopt与步骤1中的最优叶尖速比-设计叶尖速比-权重表进行比对，若最优叶尖速比λopt已经存在，则直接进行查表，更新权重；否则，进行步骤5-2；步骤5-2、将步骤4中计算得到的最佳叶尖速比λopt的新的风轮外形，基于风力机闭环系统进行动态仿真，得到来流风能分布从而确定设计叶尖速比λj和权重μj，具体步骤如下：步骤5-2-1、按照来流风能分布依次递减的顺序选择设计叶尖速比区间直到满足条件其中rtot为90％；步骤5-2-2、将步骤5-2-1中确定的区间中点作为对应设计叶尖速比λj，相应的权重系数μj由下式确定：步骤6、计算目标函数值，得到平均风能捕获效率pfavg；步骤7、判断是否满足终止条件，若不满足，通过调用遗传算法产生新一代个体，转向步骤2；若满足，则输出最优叶片的气动外形。所述终止条件是指：1)遗传代数不超过500代；或者，2)遗传算法连续计算300代，目标函数的值未变化。本发明无需限制优化风轮的λope分布特性与初始风轮，使得所得优化风轮的运行工况集中分布于更容易提升风能利用系数cp的叶尖速比处，能够进一步提高风力机的平均风能捕获效率。下面结合实施例对本发明做进一步详细的描述：实施例利用英国garradhassanpartnerslimited公司(简称gh公司)开发的通过gl(德国劳埃船级社)认证的用于风电机组设计的专业软件bladed进行仿真计算。风力机选用美国国家能源部可再生能源实验室(nationalrenewableenergylaboratory，nrel)开发的1.5mw机型，基准叶片参数见表1。表11.5mw风力机叶片参数计算pfavg所用湍流风速参数如表2所示。表2湍流风速参数参数数值平均风速(m/s)5定义速度的高度(m)84风速序列长度(sec.)3600时间步长(sec.)0.05湍流模型ieckaimal湍流强度a积分尺度150按照传统的多工况设计方法(mpcw)，根据初始风轮的λope的分布特性来确定目标函数，并且保持目标函数不变，以此来优化风轮的气动外形，并且计算出mpcw方法得到优化后风轮的设计叶尖速比λopt和平均风能捕获效率pfavg。然后，按照本发明提出的改进的多工况设计方法(mpvw)，不固定优化风轮的λope分布特性与初始风轮相同，通过风轮的集总参数λopt确定目标函数中不同的叶尖速比及其权重。计算出mpvw方法优化后风轮的设计叶尖速比λopt和平均风能捕获效率pfavg。此外，还做了一组不通过集总参数λopt，每一次都计算来流风能比例分布并且按照要求确定目标函数的仿真，也就是通过数值计算的多工况优化方法(mpnc)，也计算出优化后风轮的设计叶尖速比λopt和平均风能捕获效率pfavg。三种方法的优化过程耗时记录如表3，仿真结果记录如表4：表3优化仿真耗时优化方法mpcwmpvwmpnc仿真时间(h)6.323.5368.1表4仿真结果从表3和表4中可以看出：1)本发明提出的mpvw方法，相比于初始风轮，pfavg提升1.97％，是三种优化方法中最高的。2)未改进前的mpcw方法，虽然相比于初始风轮，pfavg虽然提升1.43％，但是本发明中改进以后的方法，能够提升更多的效率。3)mpnc方法由于未引入集总参数λopt直接进行数值计算，虽然目标函数的确定每一次都是精准的，pfavg相比于初始风轮提升1.89％，但也没有本发明中提出的mpvw方法优化出来的叶片提升效率高。此外，mpnc方法的仿真耗时远远超过mpvw方法，显然优化过程过长，时间成本过高。图2(a)是三种优化方法与初始叶片的弦长分布对比图，从图中可以看出，除了靠近叶根的部分，mpvw方法叶片的弦长比初始叶片小，节省了叶片的材料成本。mpvw和mpnc方法设计得到的叶片弦长在叶片中部(30％到70％处)很接近，这部分也是捕获功率最主要的部分。图2(b)是三种优化方法与初始叶片的扭角分布对比图，从图中可以看出，三种优化方法得到的叶片扭角很接近，都比初始叶片的扭角大。图3(a)是mpvw叶片与初始叶片的性能比较，从图中可以看出，由于两种叶片的cp-λ曲线相交于叶尖速比为6.7的位置，因此将运行叶尖速比的分布范围划分为两个区间u(5.0,6.7)和u(6.7,10.0)。尽管在区间u(5.0,6.7)中mpvw叶片的cp较初始叶片小，但是对应的来流风能分布的比例也较小，而在区间u(6.7,10.0)中来流风能分布得更多，具体的分布比例p(u(6.7,10.0))＝60.8％，相应的cp值也更高，这解释了mpvw叶片效率提升的主要原因。图3(b)是mpcw叶片与初始叶片的性能比较，从图中可以看出，两种叶片的来流风能分布是一样的，但是mpcw叶片的cp-λ曲线较初始叶片有了普遍的提升，因而效率更高。图3(c)是mpnc叶片与初始叶片的性能比较，与mpvw的分析类似，在来流风能分布更集中的设计叶尖速比区间提升风能利用效率，从而提升平均风能捕获效率。由上述实施例，可以验证本发明提出的mpvw优化方法，相比于传统的多工况优化方法，在湍流风况下能够进一步改善风力机的风能捕获效率，此外，相较于每一次优化都采用直接数值计算的方法，本发明引入风轮集总参数λopt，不仅在风能捕获效率上有了一定的提升，更重要的是减少了很多优化时间，节省了时间成本。当前第1页12