一种运行工况下H型垂直轴风力机的钝尾缘风轮防冰优化设计方法技术

本发明专利技术公开了一种运行工况下H型垂直轴风力机的钝尾缘风轮防冰优化设计方法，包括以下步骤：采用中弧线

【技术实现步骤摘要】
一种运行工况下H型垂直轴风力机的钝尾缘风轮防冰优化设计方法


[0001]本专利技术属于风力机叶片的优化设计及改型
，尤其是涉及一种利用粒子群(PSO)算法耦合计算流体力学(CFD)方法进行垂直轴风力机的钝尾缘风轮防冰优化的方法。

技术介绍

[0002]随着离网式风力发电的迅速发展，垂直轴风力机的应用越来越广泛。然而，冻雨、雪等天气条件不可避免地令叶片表面出现结冰现象。由于垂直轴风力机通常在低叶尖速比下运行，较慢的旋转速度使液滴容易进入风轮内部并与叶片内侧充分接触，从而叶片的内、外两侧均会覆盖上冰。而且，当叶片位于上风区或即使位于下风区但其线速度高于来流风速时，前缘相较于尾缘更容易结冰。结冰会增加叶片表面的粗糙度，降低风力机的风能利用率；同时，也会引起叶片静动态负载增大和风轮不平衡性增加，加速叶片的疲劳破坏。然而，叶片受到长期的交变载荷和压力负载才能产生显著的结构危害，但结冰短期内会使叶片的气动性能明显下降。因此，研究合理、有效的防冰方法，改善运行工况下垂直轴风力机叶片结冰后的气动性能，具有重要的理论意义和工程应用价值。
[0003]垂直轴风力机叶片的防冰方法包括外部处理(热能防冰、涂料防冰等)和翼型设计等。但外部处理方法具有一定的局限性，且翼型的气动性能对风能利用率存在直接的影响。因此，进行结冰环境下翼型的专门设计，例如：翼型型线的直接优化和钝尾缘改型，是减小结冰对叶片性能影响的最佳解决方案。国内外学者针对由沙尘、结冰、风蚀以及昆虫污染等造成表面粗糙的翼型开展了直接优化和钝尾缘改型相关研究。Ram等针对由于前缘污染而造成表面粗糙度增大的水平轴风力机叶片，采用多目标遗传算法耦合XFOIL软件进行根部、中部以及轮毂区域翼型的优化设计。陈进等分别采用改进的多目标PSO算法耦合RFOIL软件和遗传算法耦合RFOIL软件，进行以翼型型线的控制参数为设计变量、干净与粗糙工况下翼型的升阻比最大为设计目标的优化设计。Sun等修正升、阻力系数变化率的权重系数，建立粗糙条件下翼型性能的评价指标，并采用NSGA
‑
II算法耦合XFOIL进行以干净与粗糙工况下翼型的升阻比最大、压力系数最小为设计目标的优化设计。Zhu等引入新的形状扰动函数进行翼型的参数化表达，利用BEM理论耦合XFOIL实现以干净与粗糙条件下的升力差最大为约束条件、风轮的实度及其干净与粗糙时的功率系数加权值最大为设计目标的翼型几何优化设计。另外，杨瑞等利用在前缘处叠加细折线模拟粗糙，并分析薄、钝尾缘翼型的气动特性。张旭等通过平移前缘指定位置坐标来添加凸台以模拟粗糙，采用坐标旋转变换和缩放横纵坐标系数相结合对具有粗糙度敏感位置的翼型进行钝尾缘改型，研究改型前后翼型的升阻比、压力系数和流动特性。上述研究结果表明，优化设计和尾缘改型均使得翼型的吸力面逆压梯度减小，流动分离延迟，粗糙位置明显后移，对粗糙位置的敏感性降低，升力系数和升阻比升高。鉴于两者的优势，张旭等采用指数混合函数和B
‑
样条函数构建非对称钝尾缘翼型的参数化控制方程组，应用CFD方法、神经网络和改进的量子粒子群算法相结合优化明冰条件下钝尾缘翼型的型线，使得干净及明冰条件下的气动性能均有所改善。
[0004]到目前为止，学者们主要是针对表面具有不连续粗糙的水平轴风力机叶片开展优化设计和钝尾缘改型研究的。然而，结冰通常会连续覆盖于叶片表面且形状更加复杂。而且，水平轴风力机的结构与运行原理不同与垂直轴风力机，其结冰分布规律并不适用于垂直轴风力机。因此，国内外针对垂直轴风力机结冰进行一系列研究。李岩等采用黏土作为附着物模拟垂直轴风力机叶片表面前缘的结冰情况，并通过风洞实验分析附着物对垂直轴风力机转速和功率的影响。然而，附着物在形状、质量等方面均与冰存在差异。之后，郭文峰、李岩等利用冰风洞实验，研究垂直轴风力机在低、高叶尖速比下旋转时的结冰特性。与冰风洞实验相比，数值计算具有成本低、风险小、效率高等优点，且能提供结冰冰形相对准确的估计，在结冰研究中具有一定优势。李岩等基于二维不可压缩稳态流中的N
‑
S方程，并采用离散相位模型计算不同风速和冷水流量条件下的静态垂直轴风力机各个攻角时的翼型结冰分布。由于运行工况下垂直轴风力机叶片的相对速度随方位角时刻变化，且叶片尾流与下游叶片的相互作用会产生较为复杂的流动，因此实现非定常结冰模拟具有极大的困难。但是，可以通过指定方位角和单次结冰时间进行结冰模拟，预测准稳态下垂直轴风力机叶片的结冰形状。Manatbayev等采用多参考坐标系模型来考虑液滴场的旋转效应，进而预测不同攻角下垂直轴风力机叶片的结冰形状。Baizhuma等利用滑移网格模型和多参考坐标系模型相结合进行垂直轴风力机准稳态结冰计算，并与实验结果比较分析。
[0005]上述研究表明，优化设计和钝尾缘改型可以提高风力机叶片在干净、污染以及结冰条件下的气动性能。然而，迄今而止，对运行工况下的垂直轴风力机进行钝尾缘风轮防冰优化设计的研究并没有涉及。

技术实现思路

[0006]本专利技术要解决的问题是提供一种运行工况下H型垂直轴风力机的钝尾缘风轮防冰优化设计方法，该方法通过建立H型垂直轴风力机钝尾缘风轮的几何外形控制方程组，明确钝尾缘风轮防冰优化模型的设计变量和设计目标；利用FLUENT
‑
Solution模块和FLUENT
‑
Icing模块在每间隔45
°
方位角上进行稳态条件下结冰，获得风轮旋转一周后的冰形；采用PSO算法耦合CFD方法进行优化模型求解，提出运行工况下H型垂直轴风力机的钝尾缘风轮防冰优化设计方法，实现结冰环境下风轮的尾流长度和区域的减小以及流动分离的减弱，很大程度上提升了风力机的风能利用率。
[0007]为解决上述技术问题，本专利技术采用的技术方案是：一种运行工况下H型垂直轴风力机的钝尾缘风轮防冰优化设计方法，其特征在于：包括以下步骤：
[0008]步骤(1)，钝尾缘风轮防冰优化模型构建：采用中弧线
‑
厚度函数、坐标的旋转
‑
缩放和气动设计理论，构建H型垂直轴风力机钝尾缘风轮的几何外形控制方程组：
[0009][0010]式中：α＝arcsin(il/c)，β＝
‑
arcsin((1
‑
i)l/c)；l为钝尾缘厚度；i为钝尾缘厚度分配比；p1，a1，b1，q1，g1，h1，p2，a2，b2，q2，g2，h2为翼型的形状控制因子；c为翼型弦长；x
k
为尖尾缘翼型某控制点的横坐标；(x
″
k
，y
″
k
)为钝尾缘翼型某控制点的坐标；k＝1、2分别代表翼型的上、下面；(x
kt
、y
kt
)为原翼型最大相对厚度处的坐标，(x
′
kt
、y
′
kt
)为坐标旋转
‑
缩放后翼型最大相对厚度处的坐标；(X
k
，Y
k
，Z...

【技术保护点】

【技术特征摘要】
1.一种运行工况下H型垂直轴风力机的钝尾缘风轮防冰优化设计方法，其特征在于：包括以下步骤：步骤(1)，钝尾缘风轮防冰优化模型构建：采用中弧线
‑
厚度函数、坐标的旋转
‑
缩放和气动设计理论，构建H型垂直轴风力机钝尾缘风轮的几何外形控制方程组：式中：α＝arcsin(il/c)，β＝
‑
arcsin((1
‑
i)l/c)；l为钝尾缘厚度；i为钝尾缘厚度分配比；p1，a1，b1，q1，g1，h1，p2，a2，b2，q2，g2，h2为翼型的形状控制因子；c为翼型弦长；x
k
为尖尾缘翼型某控制点的横坐标；(x
″
k
，y
″
k
)为钝尾缘翼型某控制点的坐标；k＝1、2分别代表翼型的上、下面；(x
kt
、y
kt
)为原翼型最大相对厚度处的坐标，(x
′
kt
、y
′
kt
)为坐标旋转
‑
缩放后翼型最大相对厚度处的坐标；(X
k
，Y
k
，Z
k
)为控制点的全局坐标，全局坐标系XYZ的原点选择在风轮的几何中心；R为风轮半径；H为叶片长度且P为额定功率，ρ为空气密度，U0为额定风速，C
P
为风能利用率；i
H
为高径...

【专利技术属性】
技术研发人员：张旭，崔冷双，李伟，
申请(专利权)人：天津工业大学，
类型：发明
国别省市：