[发明专利]一种低噪声风力机叶片设计方法及低噪声风力机叶片

权利要求书

1.一种低噪声风力机叶片设计方法，其特征在于，包括以下步骤：

步骤1：构建叶片入流湍流噪声模型；

叶片入流湍流噪声频域包括高频和低频两个区域；

对于叶片高频域，叶片湍流声压级为：

其中，ρ为空气密度；c₀为声速；L为大气湍流长度因子；Δl为叶片翼段长度；为修正波动长度；Ma为来流的马赫数，是来流相对速度U与音速c₀的比值；I为湍流密度；为高频声音方向函数；r为观察者到声源距离；

湍流密度I是表面粗糙度z₀和叶片距地面高度z的方程，风力机叶片沿展向不同位置的湍流密度为：

式中γ为幂律系数，为γ＝0.24+0.096log₁₀z₀+0.016(log₁₀z₀)²；

长度因子L为表面粗糙度z₀和叶片距地面高度z的表达式：

对于低频域，叶片湍流声压级为：

其中，K_c为低频域修正因子；

步骤2：构建叶片翼型噪声模型；

包括翼型湍流边界层尾缘噪声模型及失速噪声模型；

翼型湍流边界层尾缘噪声是由湍流边界层在压力面产生的噪声S_PLp与在吸力面产生的噪声S_PLs之和：

其中，尾缘噪声是关于吸力面尾部边界层相对厚度和压力面尾部边界层相对厚度的函数，与翼型的攻角α和来流的雷诺数Re有关；S_t为斯特劳哈尔数，其中S_t1＝0.02Ma^-0.6，Ma为马赫数，S_tp、S_ts分别为压力面斯特劳哈尔数和吸力面斯特劳哈尔数；为高频声音方向函数；r为观察者到声源距离；A为频谱形状函数；W₁为振幅函数；ΔW₁为声压级修正函数；Δl为叶片翼段长度；

失速噪声为：

式中，W₂为振幅函数；B为频谱形状函数；

步骤3：构建风力机空气动力学模型；

风力机后尾流旋转时，来流风速在x、y方向的速度分量v_x、v_y为：

其中，v₀为风速，ω为风力机旋转角速度(rad/s)，r_b为叶素沿叶片展向位置，a、b为轴向诱导因子和周向诱导因子；

叶素处的入流角φ和攻角α分别为：

α＝φ-θ (11)

叶素处的合成入流速度v_rel为：

应用动量理论，推导出风力机的推力与扭矩为：

dM＝4πρωv₀bF(1-aF)r³dr (14)

式中N表示风力机叶片数，R为叶片长度；

应用叶素理论，推力和扭矩为：

其中C_n、C_t分别为法向力系数和切向力系数；式15中的F₁为法向力系数修正因子，式16中的F₁为切向力系数修正因子：

g＝exp[-0.125(Nλ-21)]+0.1 (18)

其中λ为叶尖速比；

联立(6)～(9)，求出轴向诱导因子a和周向诱导因子b为：

式中Y₂＝4Fsinφcosφ/(σC_tF₁)；σ＝Nc/(2πr)，c表示叶片翼段弦长；

步骤4：计算叶片气动噪声；

步骤5：构建叶片气动噪声优化数学模型；

步骤5.1：构建目标函数f(x)；

f(x)＝max(C_p/SPL_total) (21)

其中，C_p为风力机功率系数，SPL_total为叶片气动噪声；

风力机功率系数C_p为：

其中，P表示风力机功率；

基于动量叶素理论，经过微积分变换，风力机功率系数另一种表达式为：

叶片沿展向看作是由若干翼段组成，每个翼段噪声计算采用叶片入流湍流噪声计算公式及叶片翼型噪声计算公式预测；对于第i个翼段，其噪声计算公式为：

其中，j表示不同噪声源，包括叶片入流湍流噪声、翼型尾缘噪声及翼型失速噪声；K_A为加权过滤值(dB)；

则整个叶片噪声声压级预测由若干翼段噪声声压级或者功率级叠加：

步骤5.2：设计变量及约束条件；

选取叶片弦长及扭角沿叶片展向变化作为设计变量，其约束范围为：

	弦长c/m	扭角θ(°)
最大值	4.0	16.00
最小值	0.0	-6.00

对风力机叶根挥舞弯矩M_flap及扭矩M_T进行约束；

其中，M_flap,max表示风力机叶根挥舞弯矩的最大值，M_T,max表示风力机叶根挥舞扭矩的最大值；

步骤5.3：采用粒子群算法对步骤5.1中的目标函数进行优化。