CN109409019A - 一种风力机叶片气动外形的预弯优化方法 - Google Patents
一种风力机叶片气动外形的预弯优化方法 Download PDFInfo
- Publication number
- CN109409019A CN109409019A CN201811537228.9A CN201811537228A CN109409019A CN 109409019 A CN109409019 A CN 109409019A CN 201811537228 A CN201811537228 A CN 201811537228A CN 109409019 A CN109409019 A CN 109409019A
- Authority
- CN
- China
- Prior art keywords
- bending
- blade
- pneumatic equipment
- function
- equipment bladess
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000005452 bending Methods 0.000 title claims abstract description 58
- 238000005457 optimization Methods 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012549 training Methods 0.000 claims abstract description 29
- 230000004044 response Effects 0.000 claims abstract description 18
- 238000013461 design Methods 0.000 claims abstract description 15
- 238000006243 chemical reaction Methods 0.000 claims abstract description 11
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 238000013507 mapping Methods 0.000 claims abstract description 8
- 238000003050 experimental design method Methods 0.000 claims abstract description 6
- 239000003795 chemical substances by application Substances 0.000 claims description 15
- 238000009826 distribution Methods 0.000 claims description 11
- 238000005314 correlation function Methods 0.000 claims description 9
- 238000013401 experimental design Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000009792 diffusion process Methods 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000031068 symbiosis, encompassing mutualism through parasitism Effects 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 239000012141 concentrate Substances 0.000 abstract 1
- 238000005070 sampling Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000003321 amplification Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003365 glass fiber Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Abstract
本发明公开了一种风力机叶片气动外形的预弯优化方法,属于优化方式领域,包含以下步骤,根据风力机叶片的优化目标,确定气动外形参数的优化变量及其范围,采用最优化拉丁超立方试验设计方法进行训练样本点的采样;利用数值模拟方法求解训练样本点对应的响应值;利用Kriging代理模型为近似模型方法,求得模型的超参数θ和设计变量与响应的映射关系;求解期望函数的最大值,获得下一个采样点的位置;求解新样本点的响应值,并加入到初始样本点集中,最终求得期望函数的最优解满足收敛条件。基于该方法优化预弯得到的叶片,即可以有效增大叶尖与塔架之间的间隙,又使得风力机输出功率得到明显的提升。
Description
技术领域
本发明涉及优化方法领域,更具体的,涉及一种风力机叶片气动外形的预弯优化方法。
背景技术
风力机额定功率不断增大,风力机叶片的尺寸也随之增大。现代大型风力机叶片制造材料均采用玻璃纤维,导致叶片刚度明显降低,促使其柔性进一步增大。为了避免风力机在运行过程中,特别在强风作用下发生叶片弯曲变形与塔架碰撞,有公司提出一种避免叶片尖部与塔架碰撞的策略,即叶片向迎风方向提前进行预先弯曲处理,可以有效增大叶尖与塔架之间的间隙。叶片的柔性预弯处理,可以降低局部应力载荷,减少生产叶片的材料,从而降低叶片的制造成本。但是有学者指出,预弯叶片的输出功率小于原有叶片的输出功率,使得风力机整体的功率降低。
发明内容
为了克服现有技术的缺陷,本发明所要解决的技术问题在于提出提出一种风力机叶片气动外形的预弯优化方法,其能够克服目前叶片技术的不足,提供一种风力机叶片气动外形的预弯优化方法,使叶片预弯后的输出功率得到一定的提升。
为达此目的,本发明采用以下技术方案:
本发明提供了一种风力机叶片气动外形的预弯优化方法,包括以下步骤:
步骤1:确定风力机叶片的优化目标,以叶片扭矩为优化目标;
步骤2:分析影响叶片扭矩的因素,得出其主要受预弯展向位置r/R和预弯距离Δz影响,以预弯展向位置r/R和预弯距离Δz为设计因子,即叶片的优化参数,利用最优化拉丁超立方试验设计方法采集优化参数的训练样本点;
步骤3:采用数值模拟方法获得训练样本点对应的响应值;
步骤4:采用模式算法与遗传算法的混合算法改进Kriging代理模型,求得模型的超参数θ的初始解;
步骤5:构建Kriging代理模型,求得模型的超参数θ以及设计因子与响应值的映射关系;
步骤6:采用期望函数判断映射关系是否收敛,求出期望函数的最大值,判断其是否满足收敛条件,若满足收敛条件,即优化成功,如不满足收敛条件,求解期望函数的最大值,以期望函数的最大值所对应的设计因子为下一训练样本点;
步骤7:将下一样本点加入到步骤2的求得的初始样本点集中,回到步骤3进行循环,直到期望函数的最优解满足收敛条件。
在本发明较佳地技术方案中,所述步骤1中,风力机的输出功率为P=2π×Ω0×MTorque×B/60;其中,Ω0为叶轮转速,单位为r/min,MTorque,为叶片扭矩,单位为,N·m,B为叶片数;运行工况不变情况下,所述输出功率只受所述叶片扭矩影响,则以所述叶片扭矩为优化目标。
在本发明较佳地技术方案中,所述步骤2中,以叶片尖部到预弯展向位置的距离得到预弯展向位置r/R,采用多项式插值得到截面翼型的预弯距离Δz。拉丁超立方试验设计可以按照以下公式选取:其中,1≤j≤n,1≤i≤l,U是0到1之间的随机数值,π为0,1,…,l-1的独立随机排列,共生成20个训练样本点。
在本发明较佳地技术方案中,所述步骤5中Kriging代理模型由多项式和随机分布组成,具体模型为:其中,βp×1是回归系数向量,f(x)p×1 T是已知的回归模型,为多项式函数,可以为常数、一阶或二阶函数,f(x)p×1 Tβp×1是已知部分,对试验设计空间全局近似。
在本发明较佳地技术方案中,z(x)是随机误差,服从正态分布,协方差非零,具有统计特性。E[z(x)]=0
其中,xi和xj是两个任意的训练样本点,θ是各向异性相关参数,R(θ,xi,xj)是训练样本点的相关函数。
在本发明较佳地技术方案中,R为相关函数矩阵,其表达式为:
其中m=20。已知的回归函数模型通常选取一阶线性多项式模型。
在本发明较佳地技术方案中,所述一阶线性多项式的形式为:
随机函数z(x)服从正态分布,的对数似然函数为:
问题可以转化求解一个非线性的无约束全局优化问题min[Ψ(θ)]≡σa(θ)2|R|。
在本发明较佳地技术方案中,所述步骤6中定义fmax=max(y1(x),y2(x),…,yn(x))为训练样本点中的最大响应值定义fmax=max(y1(x),y2(x),…,yn(x))为训练样本点中的最大响应值,而期望函数的定义为:
其中,为预测点x的响应值,s为预测点x的标准差,Φ为标准正太分布函数,Ψ为标准正太分布的概率密度函数。
在本发明较佳地技术方案中,所述步骤3中数值模拟采用的计算几何模型为1.5MW水平轴风力机叶片,采用非定常的雷诺三维不可压N-S方程,湍流模型应用SST k-ω湍流模型。对于风力机叶片三维流场的数值模拟,取得每个所述训练样本单所对应的响应值,该响应值为所述训练样本点对应的叶片扭矩。
在本发明较佳地技术方案中,所述步骤3中,对叶片表面以及叶片附近区域进行了网格加密处理,计算域的进口边界给定为速度入口条件,出口边界给定为压力出口条件,并给出平均静压,叶片的表面采用无滑移边界条件。速度方程和压力方程的耦合求解应用SIMPLEC算法,对流项应用二阶迎风格式,扩散项应用中心差分格式。
本发明的有益效果为:
本发明提供的一种风力机叶片气动外形的预弯优化方法,将叶片气动性能预测方法与Kriging代理模型相结合,以及以预弯展向位置r/R和预弯距离Δz为设计因子,经过模型的筛选映射,最终输出最优解,基于该方法优化预弯得到的叶片,在叶片可以弯曲设置的同时,叶片扭矩也会得到提升,在工况不变的情况下提高风力机的输出功率。
附图说明
图1是本发明具体实施方式提供的风力机叶片气动外形的预弯优化方法的建模流程图;
图2是本发明具体实施方式提供的一种风力机叶片气动外形预弯优化方法的优化预弯叶片与原叶片气动外形对比图;
图3是本发明具体实施方式提供的一种风力机叶片气动外形预弯优化方法的扭矩系数分布对比图;
图4是本发明具体实施方式提供的一种风力机叶片气动外形预弯优化方法的推力系数分布对比图。
具体实施方式
下面结合附图并通过具体实施方式来进一步说明本发明的技术方案。
如图1所示,实施例中提供了一种风力机叶片气动外形的预弯优化方法,所述建模方法包括以下步骤:
步骤1:确定风力机叶片的优化目标,以叶片扭矩为优化目标。
风力机的输出功率表达式为P=2π×Ω0×MTorque×B/60
其中,Ω0为叶轮转速,单位为r/min,MTorque为叶片扭矩,单位为N·m,B为叶片数。在风力机的运行工况不变的情况下,即叶轮转速以及叶片数不变的情况下,影响风力机输出功率的因素只有叶片扭矩,因此可以确定叶片扭矩为优化目标。
步骤2:分析影响叶片扭矩的因素,得出其主要受预弯展向位置r/R和预弯距离Δz影响,以预弯展向位置r/R和预弯距离Δz为设计因子,即叶片的优化参数,利用最优化拉丁超立方试验设计方法采集优化参数的训练样本点。本小节以叶片的扭矩为优化目标,即响应值。叶片的尺寸保持不变的情况下,以预弯展向位置r/R和预弯距离Δz为设计因子。在风力机叶片的预弯优化过程中,设计变量与目标函数之间是隐函数关系,并且是高度非线性。根据最优化拉丁超立方试验设计方法,共生成20个训练样本点。以叶片尖部到预弯展向位置的距离得到预弯展向位置r/R,采用多项式插值方法得到截面翼型的预弯距离Δz。
根据以上的介绍,本小节风力机叶片的预弯优化过程,是以叶片扭矩为优化目标,约束条件为1m≤Δz≤2m、0.25≤r/R≤0.75和R=35m,这个约束条件即为预测点。则优化问题可以用如下形式表示:
拉丁超立方试验设计一般可以按照以下公式选取:
其中,1≤j≤n,1≤i≤l,U是0到1之间的随机数值,π为0,1,…,l-1的独立随机排列。该方法通过选取较少的样本点,从而反映整个试验设计的空间属性。
步骤3:采用数值模拟方法获得训练样本点对应的响应值;
采用的计算几何模型是美国可再生能源实验室WindPACT工程设计的1.5MW水平轴风力机叶片。采用非定常的雷诺三维不可压N-S方程,湍流模型应用SSTk-ω湍流模型。对于风力机叶片三维流场的数值模拟,主要研究对象为风力机叶片,故对叶片表面以及叶片附近区域进行了网格加密处理,可以确保得到准确的计算结果。计算域的进口边界给定为速度入口条件,出口边界给定为压力出口条件,并给出平均静压,叶片的表面采用无滑移边界条件。速度方程和压力方程的耦合求解应用SIMPLEC算法,对流项应用二阶迎风格式,扩散项应用中心差分格式。
步骤4:采用模式算法与遗传算法的混合算法改进Kriging代理模型,求得模型的超参数θ的初始解;
基于Kriging代理模型和改进的Kriging代理模型,分别拟合测试函数F1(x)、F2(x)和F3(x),利用最优化拉丁超立方试验设计方法分别对测试函数F1(x)、F2(x)和F3(x)选取了33个、60个和60个训练样本点。表1对比了Kriging代理模型和改进的Kriging代理模型的均方根误差。改进的Kriging代理模型均能准确的预测测试函数的分布,三个函数的均方根误差均有所下降,其中,Schaffer函数分布简单,其均方根误差降低幅度最大。
(1)Branin函数
其中,x1∈[-5,10],x2∈[0,15]。
(2)Schaffer函数
其中,x1∈[-2,2],x2∈[-2,2]。
(3)Ackley函数
其中,x1∈[0,2],x2∈[0,2]。
表1
步骤5:构建Kriging代理模型,求得模型的超参数θ和设计变量与响应的映射关系;
Kriging代理模型由多项式和随机分布组成,具体模型为:
其中,βp×1是回归系数向量;f(x)p×1 T是已知的回归模型,多为多项式函数,可以为常数、一阶或二阶函数;f(x)p×1 Tβp×1是已知部分,对试验设计空间全局近似;z(x)是随机误差,服从正态分布,协方差非零,具有统计特性。
E[z(x)]=0
其中,xi和xj是两个任意的训练样本点,θ是各向异性相关参数,R(θ,xi,xj)是训练样本点的相关函数,R是相关函数矩阵,其表达式为:
其中m=20,已知的回归函数模型通常选取常数、一阶和二阶线性多项式模型。本发明中,选择一阶线性多项式模型。
相关函数的形式为:
当确定了回归函数模型和相关函数,相关函数矩阵R、广义最小二乘估计β*和预测方差的最大似然估计σa 2均与各向异性相关参数θ有关。随机函数z(x)服从正态分布,则的对数似然函数为:
问题可以转化求解一个非线性的无约束全局优化问题:
min[Ψ(θ)]≡σa(θ)2|R|
步骤6:采用期望函数判断映射关系是否收敛,求出期望函数的最大值,判断其是否满足收敛条件,若满足收敛条件,即优化成功,如不满足收敛条件,求解期望函数的最大值,以期望函数的最大值所对应的设计因子为下一训练样本点。定义fmax=max(y1(x),y2(x),…,yn(x))为训练样本点中的最大响应值。EI(x)函数定义为:
其中,为预测点x的响应值,s为预测点x的标准差,Φ为标准正太分布函数,Ψ为标准正太分布的概率密度函数。这些参数均根据映射关系得到。当。求出EI(x)的最大值EImax,然后将EImax与收敛条件Δs=0.001比较,如果EImax小于Δs=0.001,则期望函数收敛,优化过程收敛,优化成功,这部分训练样本点达到优化效果。如果EImax大于于Δs=0.001,期望函数不收敛,则需要以EImax所对应的设计因子为下一训练样本点。
步骤7:将下一样本点加入到步骤2的求得的初始样本点集中,回到步骤3进行循环,直到期望函数的最优解满足收敛条件。
如图2以及表2所示,优化后,叶片向迎风侧预弯约为1.25m,叶片在r/R=0.35处开始向迎风方向进行了弯曲。右侧为美国可再生能源实验室的WindPACT叶片,即优化前的叶片,左侧为预弯叶片,即优化后的叶片。由于预弯叶片向迎风方向弯曲,叶尖到旋转轴线的距离减小了约0.03m。优化后叶轮的扭矩提升了约为0.33%,即预弯叶片在发电能力方面略微高于直叶片的发电能力,同时预弯叶片的功率系数大于原直叶片的功率系数。因此,预弯叶片向迎风方向进行一定的弯曲,可以增大叶尖与塔架之间的间隙,确保的风力机在强风作用下避免叶片尖部与塔架碰撞,同时叶轮的扭矩得到增加。
表2
如图3所示,弯叶片和原直叶片的扭矩系数和推力系数均呈现先上升后下降的趋势。在r/R=0.2处,两种叶片截面的扭矩系数为最大值。优化后叶片截面的当地扭矩系数在展向0.2<r/R<0.3和r/R>0.55均有一定的增大,当地扭矩的最大增幅处在r/R=0.20处,增幅比例为7.10%.在叶片展向位置的中部,优化后的扭矩系数略有减小。
如图4所示,弯叶片的原直叶片的推力系数均呈现先上升后下降的趋势,在r/R=0.45处两种叶片,两种叶片的推力达到最大值。弯叶片的推力一直大于原直叶片,因此可以体现出优化后的叶片是具有增大功率的效果。
本发明是通过优选实施例进行描述的,本领域技术人员知悉,在不脱离本发明的精神和范围的情况下,可以对这些特征和实施例进行各种改变或等效替换。本发明不受此处所公开的具体实施例的限制,其他落入本申请的权利要求内的实施例都属于本发明保护的范围。
Claims (10)
1.一种风力机叶片气动外形的预弯优化方法,其特征在于:
包括以下步骤:
步骤1:确定风力机叶片的优化目标,以叶片扭矩为优化目标;
步骤2:分析影响叶片扭矩的因素,得出其主要受预弯展向位置r/R和预弯距离Δz影响,以预弯展向位置r/R和预弯距离Δz为设计因子,即风力机叶片的优化参数,利用最优化拉丁超立方试验设计方法采集优化参数的训练样本点;
步骤3:采用数值模拟方法获得训练样本点对应的响应值;
步骤4:采用模式算法与遗传算法的混合算法改进Kriging代理模型,求得模型的超参数θ的初始解;
步骤5:构建Kriging代理模型,求得模型的超参数θ以及设计因子与响应值的映射关系;
步骤6:采用期望函数判断映射关系是否收敛,求出期望函数的最大值,判断其是否满足收敛条件,若满足收敛条件,即优化成功,如不满足收敛条件,求解期望函数的最大值,以期望函数的最大值所对应的设计因子为下一训练样本点;
步骤7:将下一样本点加入到步骤2的求得的初始样本点集中,回到步骤3进行循环,直到期望函数的最优解满足收敛条件。
2.根据权利要求1所述的一种风力机叶片气动外形的预弯优化方法,其特征在于:
所述步骤1中,风力机的输出功率为P=2π×Ω0×MTorque×B/60;其中,Ω0为叶轮转速,单位为r/min,MTorque,为叶片扭矩,单位为,N·m,B为叶片数;运行工况不变情况下,所述输出功率只受所述叶片扭矩影响,则以所述叶片扭矩为优化目标。
3.根据权利要求1所述的一种风力机叶片气动外形的预弯优化方法,其特征在于:
所述步骤2中,以叶片尖部到预弯展向位置的距离得到预弯展向位置r/R,采用多项式插值得到截面翼型的预弯距离Δz;
拉丁超立方试验设计可以按照以下公式选取:其中,1≤j≤n,1≤i≤l,U是0到1之间的随机数值,π为0,1,…,l-1的独立随机排列,共生成20个训练样本点。
4.根据权利要求1所述的一种风力机叶片气动外形的预弯优化方法,其特征在于:
所述步骤5中Kriging代理模型由多项式和随机分布组成,具体模型为:
其中,βp×1是回归系数向量;
f(x)p×1 T是已知的回归模型,为多项式函数,可以为常数、一阶或二阶函数;
f(x)p×1 Tβp×1是已知部分,对试验设计空间全局近似。
5.根据权利要求4所述的一种风力机叶片气动外形的预弯优化方法,其特征在于:
z(x)是随机误差,服从正态分布,协方差非零,具有统计特性;
E[z(x)]=0;
其中,xi和xj是两个任意的训练样本点,θ是各向异性相关参数,R(θ,xi,xj)是训练样本点的相关函数。
6.根据权利要求5所述的一种风力机叶片气动外形的预弯优化方法,其特征在于:
R为相关函数矩阵,其表达式为:
其中m=20;
已知的回归函数模型通常选取一阶线性多项式模型。
7.根据权利要求6所述的一种风力机叶片气动外形的预弯优化方法,其特征在于:
所述一阶线性多项式的形式为:
随机函数z(x)服从正态分布,的对数似然函数为:
问题可以转化求解一个非线性的无约束全局优化问题min[Ψ(θ)]≡σa(θ)2|R|。
8.根据权利要求7所述的一种风力机叶片气动外形的预弯优化方法,其特征在于:
所述步骤6中定义fmax=max(y1(x),y2(x),…,yn(x))为训练样本点中的最大响应值,而期望函数的定义为:
其中,为预测点x的响应值,s为预测点x的标准差,Φ为标准正太分布函数,Ψ为标准正太分布的概率密度函数。
9.根据权利要求1所述的一种风力机叶片气动外形的预弯优化方法,其特征在于:
所述步骤3中数值模拟采用的计算几何模型为1.5MW水平轴风力机叶片,采用非定常的雷诺三维不可压N-S方程,湍流模型应用SST k-ω湍流模型,对于风力机叶片三维流场的数值模拟,取得每个所述训练样本单所对应的响应值,该响应值为所述训练样本点对应的叶片扭矩。
10.根据权利要求9所述的一种风力机叶片气动外形的预弯优化方法,其特征在于:
所述步骤3中,对叶片表面以及叶片附近区域进行了网格加密处理,计算域的进口边界给定为速度入口条件,出口边界给定为压力出口条件,并给出平均静压,叶片的表面采用无滑移边界条件,速度方程和压力方程的耦合求解应用SIMPLEC算法,对流项应用二阶迎风格式,扩散项应用中心差分格式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811537228.9A CN109409019A (zh) | 2018-12-15 | 2018-12-15 | 一种风力机叶片气动外形的预弯优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811537228.9A CN109409019A (zh) | 2018-12-15 | 2018-12-15 | 一种风力机叶片气动外形的预弯优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109409019A true CN109409019A (zh) | 2019-03-01 |
Family
ID=65459747
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811537228.9A Pending CN109409019A (zh) | 2018-12-15 | 2018-12-15 | 一种风力机叶片气动外形的预弯优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109409019A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112765743A (zh) * | 2021-01-13 | 2021-05-07 | 国电联合动力技术有限公司 | 一种风电机组整机一体化设计多目标协同优化方法及装置 |
CN113378330A (zh) * | 2021-07-15 | 2021-09-10 | 河海大学 | 用于超长柔性下风向风力机叶片有效扫风半径的计算方法 |
CN114065572A (zh) * | 2021-10-27 | 2022-02-18 | 西北工业大学 | 一种厚度偏差与叶片厚度分布相关的叶片模型构建方法 |
CN116776600A (zh) * | 2023-06-21 | 2023-09-19 | 安徽工业大学 | 基于自适应代理模型的风力机叶片优化设计方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2253836A1 (en) * | 2009-05-18 | 2010-11-24 | Lm Glasfiber A/S | Wind turbine blade |
CN102750410A (zh) * | 2012-06-12 | 2012-10-24 | 中国科学院工程热物理研究所 | 一种水平轴风力机叶片铺层的优化设计方法 |
CN102867083A (zh) * | 2012-08-30 | 2013-01-09 | 浙江大学 | 考虑不确定性的压力机滑块机构高刚度轻量化设计方法 |
CN108549773A (zh) * | 2018-04-16 | 2018-09-18 | 西北工业大学 | 网格参数化方法及基于该网格参数化方法的涡轮叶片多学科可靠性设计优化方法 |
-
2018
- 2018-12-15 CN CN201811537228.9A patent/CN109409019A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2253836A1 (en) * | 2009-05-18 | 2010-11-24 | Lm Glasfiber A/S | Wind turbine blade |
CN102750410A (zh) * | 2012-06-12 | 2012-10-24 | 中国科学院工程热物理研究所 | 一种水平轴风力机叶片铺层的优化设计方法 |
CN102867083A (zh) * | 2012-08-30 | 2013-01-09 | 浙江大学 | 考虑不确定性的压力机滑块机构高刚度轻量化设计方法 |
CN108549773A (zh) * | 2018-04-16 | 2018-09-18 | 西北工业大学 | 网格参数化方法及基于该网格参数化方法的涡轮叶片多学科可靠性设计优化方法 |
Non-Patent Citations (1)
Title |
---|
王海鹏: "风力机非定常气动特性及优化设计研究", 《中国博士学位论文全文数据库(电子期刊) 工程科技Ⅱ辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112765743A (zh) * | 2021-01-13 | 2021-05-07 | 国电联合动力技术有限公司 | 一种风电机组整机一体化设计多目标协同优化方法及装置 |
CN113378330A (zh) * | 2021-07-15 | 2021-09-10 | 河海大学 | 用于超长柔性下风向风力机叶片有效扫风半径的计算方法 |
CN113378330B (zh) * | 2021-07-15 | 2023-09-08 | 河海大学 | 用于超长柔性下风向风力机叶片有效扫风半径的计算方法 |
CN114065572A (zh) * | 2021-10-27 | 2022-02-18 | 西北工业大学 | 一种厚度偏差与叶片厚度分布相关的叶片模型构建方法 |
CN116776600A (zh) * | 2023-06-21 | 2023-09-19 | 安徽工业大学 | 基于自适应代理模型的风力机叶片优化设计方法及系统 |
CN116776600B (zh) * | 2023-06-21 | 2024-04-12 | 安徽工业大学 | 基于自适应代理模型的风力机叶片优化设计方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109409019A (zh) | 一种风力机叶片气动外形的预弯优化方法 | |
CN109460566B (zh) | 一种风力机叶片内侧厚翼型的气动稳健优化设计方法 | |
Zanette et al. | A design methodology for cross flow water turbines | |
CN109190283A (zh) | 一种考虑高湍流自由来流效应的风力机翼型气动优化方法 | |
Maalawi et al. | A practical approach for selecting optimum wind rotors | |
CN109145506B (zh) | 一种高气动性能低噪声水平风力机外侧翼型的优化设计方法 | |
CN112507471B (zh) | 一种大攻角范围条件下垂直轴风力机翼型设计方法 | |
CN101923584B (zh) | 风力机专用翼型设计方法及风力机专用翼型 | |
CN101615216B (zh) | 对翼型尾缘进行流线形增厚的方法 | |
CN106126791B (zh) | 一种考虑几何不确定性的高超声速机翼气动力/热分析方法 | |
CN114091377B (zh) | 基于空间变化的动态双高斯风力机尾流风速的计算方法 | |
Bangga et al. | The effects of airfoil thickness on dynamic stall characteristics of high‐solidity vertical axis wind turbines | |
Liu et al. | Optimization model for rotor blades of horizontal axis wind turbines | |
CN106777567A (zh) | 一种风力机钝尾缘翼型优化设计方法 | |
CN115062563A (zh) | 基于空间变化的风力机三维尾流风速计算方法 | |
CN109783912B (zh) | 一种升力型垂直轴风力发电机叶片翼型优化设计方法 | |
CN106919730B (zh) | 一种采用风速衰减因子的风电场尾流确定方法 | |
Tanürün | Improvement of vertical axis wind turbine performance by using the optimized adaptive flap by the Taguchi method | |
Hao et al. | Aerodynamic and aeroacoustic optimization of wind turbine blade by a genetic algorithm | |
CN116108574A (zh) | 一种漂浮式风电叶片模型气动设计方法及模型 | |
Lee et al. | Lift correction model for local shear flow effect on wind turbine airfoils | |
CN113742861A (zh) | 一种适用于风力发电机风洞试验的叶片模型优化设计方法 | |
CN106227985A (zh) | 潮流能水平轴水轮机用叶片翼型族设计方法 | |
Simpson et al. | Influence of tower shadow on downwind flexible rotors: Field tests and simulations | |
CN112906166B (zh) | 一种考虑气动效率与气动载荷的风力机叶片优化设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190301 |