CN115510381A - 一种海上风机多元相干效应风场载荷构建方法 - Google Patents
一种海上风机多元相干效应风场载荷构建方法 Download PDFInfo
- Publication number
- CN115510381A CN115510381A CN202211181493.4A CN202211181493A CN115510381A CN 115510381 A CN115510381 A CN 115510381A CN 202211181493 A CN202211181493 A CN 202211181493A CN 115510381 A CN115510381 A CN 115510381A
- Authority
- CN
- China
- Prior art keywords
- wind
- wind turbine
- load
- effect
- rotating
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Abstract
本发明涉及海上风场构建技术领域,具体涉及一种海上风机多元相干效应风场载荷构建方法,其包括以下步骤:一、检测海上风机的物理量和海上风机所在风场的物理量;二、进行时空等效;三、建立旋转叶片的频域旋转效应风谱;四、将海上风机所在风场的风场域划分为转动域、静止域和风浪作用域;五、进行等效变换,建立等效静态相干函数;六、构建互功率谱,生成多元相干效应联合功率谱矩阵;七、计算适应于海上风机的风速幅值和相位参数;八、通过分段评估进一步构建海上风机多元相干效应全域载荷,本申请具有能够结合复杂的海况来计算海上风机多元相干效应风场载荷的效果。
Description
技术领域
本发明涉及海上风场构建技术领域,具体涉及一种海上风机多元相干效应风场载荷构建方法。
背景技术
目前海上风能作为一种重要的战略资源已经引起了广泛的关注,成为不可或缺的能源组份。与陆地环境相比,海上的风能资源更加丰富,具有广阔的开发前景。随着海上风机高功率化的发展,海上风机叶片呈现出大细长比的特点,这导致风机叶片的旋转扰动效应更加明显。通常,海工结构依据《建筑结构荷载规范》,通过建立风振系数、风载荷体形系数、风压高度变化系数和基本风压函数关系的方式来计算垂直于海工结构的风载荷标准值。但这种风载荷计算方式仅适应于固定式海工结构,若按照固定式海工结构的脉动风进行气动载荷的求解将不能反映风机叶片的转动倍频效应。在海上环境载荷构建方面,风浪载荷常被分开模拟计算,这导致风场模型只包含不同脉动风速位置之间的空间相干性。然而海上风机的环境荷载是多域关联耦合的,叶片旋转域、塔架静止域以及相应的风浪作用域都具有相干性。因此在进行海上风机风场构建时,需要先构建适应风机叶片的旋转效应风谱,其次考虑多域载荷的空间相干性,进而才能准确评估海上风机的全域耦合风场。
现有一申请公开号为CN111985018A的发明专利,其保护了一种基于惯性力法和塔线分离法考虑塔线耦合影响的超高大跨越塔、线风载荷的计算方法,步骤为:搭建超高大跨越塔的塔线体系,获取塔线体系的物理参数;基于塔线耦合影响因子,求取杆塔等效阻尼系数、塔线体系超高大跨越塔的风振系数、塔线体系悬垂绝缘子串最大风偏角的风振系数、塔线体系风荷载脉动折减系数;对塔线体系超高大跨越塔的风振系数、风偏角的风振系数进行修正计算,得到塔线体系超高大跨越塔的修正风振系数和塔线体系输电线的修正风振系数;基于塔线分离法,在等效振动惯性力作用下计算超高大跨越塔线体系中超高输电塔的设计风荷载和大跨越输电线的设计风荷载。
上述技术方案中,由于海上尤其深海区域环境复杂多样,极端海况发生的概率更高,这导致风浪相互作用更加明显,载荷计算需要考虑更多的要素。
发明内容
为解决现有技术存在的不足,本发明提供了一种海上风机多元相干效应风场载荷构建方法。
本发明的技术方案为:
本发明提供了一种海上风机多元相干效应风场载荷构建方法,包括以下步骤:
一、检测海上风机的物理量和海上风机所在风场的物理量;
二、通过泰勒冻结原理,对海上风机叶片上不同空间位置和不同时刻的风速点进行时空等效;
三、结合相关性原理与狄拉克函数的性质建立旋转叶片的频域旋转效应风谱;
四、将海上风机所在风场的风场域划分为转动域、静止域和风浪作用域;
五、基于叶片转动建立转动域动态相干性函数,对转动坐标系下旋转风速与固定坐标系下脉动风速进行等效变换,建立海上风机叶片转动与塔架静止域的等效静态相干函数;
六、通过风浪相干函数构建塔架静止与波浪作用域、叶片转动与波浪作用域的互功率谱,结合旋转域的互功率谱、转动与静止域的互功率谱生成多元相干效应联合功率谱矩阵;
七、使用改进平方根法对多元相干效应联合功率谱矩阵进行分解,根据谐波叠加原理重构多元相干效应风速和波高时程,计算适应于海上风机的风速幅值和相位参数;
八、通过分段评估进一步构建海上风机多元相干效应全域载荷,海上风机多元相干效应全域载荷包括旋转域空气动力载荷、静止域风载荷以及波浪载荷,求解获得海上风机多元相干效应载荷。
本发明所达到的有益效果为:将海上风场划分为多元相干效应风场,在转动域引入相关函数原理提出具有物理表征意义的频域旋转效应风谱,准确识别叶片转动产生的1p-3p倍频效应,进而能准确反映海上风机结构动力响应的多模态振动信息,将转动坐标系下的旋转效应风谱转化为静止坐标系下的脉动风谱来实现动静转化技术,建立转动、静止以及风浪作用域的空间相干性函数,引入多元相干效应联合功率谱矩阵准确构建海上风机全域耦合风场。通过对风场域的多元的载荷计算,获得准确的海上风机多元相干效应载荷。
进一步,所述通过泰勒冻结原理,对海上风机叶片上不同空间位置和不同时刻的风速点进行时空等效还包括:
当风速点i和j分别位于静止位置与转动位置时,将静止位置风速点间的互相关函数与脉动风互功率谱构成傅里叶变换对,旋转位置风速点间的互相关函数与旋转风互功率谱构成傅里叶变换对:
其中Sii(f)和Sjj(f)表示点i和j的脉动风自功率谱,Coh(i,j;f)是i和j之间的相干函数,θij表示相位延迟参数。
进一步,所述结合相关性原理与狄拉克函数的性质建立旋转叶片的频域旋转效应风谱还包括:
选定f0作为叶片的转动频率,则相干函数可以表示为展开系数An(i,j;f)和f0的表达式,即:
其中An(i,j;f)为相干函数的傅里叶展开系数,φ0表示叶片的分布角;
基于傅里叶展开系数,旋转效应风谱可进一步展开为:
设置狄拉克函数δ为:
通过旋转效应风谱和狄拉克函数δ获得旋转效应风谱与转动效应模态、源谱平移效应和转动相位延迟的相关函数:
取相位延迟参数为0,并令θ=2πf0τ,则适配海上风机叶片的旋转效应自功率谱与互功率谱表示为:
Gauss-Lejeune积分公式由高斯系数和高斯点所代表的函数之和近似表达为:
式中,tj为高斯点,Aj表示高斯系数,通过上述公式求解旋转风速。
通过上述方案,从相关函数的角度构建了适配海上风机叶片的频域旋转效应风谱模型,模型具有物理表征意义,可根据叶片旋转频率重新分配1p~3p 频率处的波峰能量,Gauss-Lejeune方法的引入进一步解决了旋转效应风谱的震荡积分收敛问题,有助于为随后多元相干效应风场构建提供理论前提。
进一步,所述将海上风机所在风场的风场域划分为转动域、静止域和风浪作用域还包括:
将海上风机最上面承受的海上环境载荷划分为旋转载荷,将海上风机中间承受的海上环境载荷划分为静风载荷,将海上风机下面承受的海上环境载荷划分为波浪载荷,将海上风机所在风场的风场域对应旋转载荷所在区域划分为转动域,将海上风机所在风场的风场域对应静风载荷所在区域划分为静止域,将海上风机所在风场的风场域对应波浪载荷所在区域划分为风浪作用域。
通过上述方案,根据海上风机中间承受的海上环境载荷的性质来划分风场域,实现了全域载荷划分,使得最终计算结果能够考虑所有载荷的影响。
进一步,所述步骤五还包括:
使用Davenport相干函数表示海上风机叶片转动时叶片上不同点i和j的动态相干性:
其中,ri和rj分别表示叶片上i和j点到轮毂中心的距离;B1→B1和B1→B2分别表示i和j点位于相同和不同叶片上;
海上风机叶片o点和塔架p点的跨区域风速互功率谱用旋转效应风谱和静止域相干函数表示为:
通过上述方案,在旋转域建立了叶片转动频率相关的动态相干性函数,通过动静转换原则建立了旋转与静止域的等效静态相干函数,将海上风机波面以上的区域关联起来,并以此构建跨区域风速互功率谱。
进一步,所述步骤六还包括:
选用塔架p点和海平面q点建立塔架静止与波浪作用域的互功率谱:
式中,Tp表示谱峰周期且fp=1/Tp,Hs是有效波高,σ是峰形系数,γ是谱峰升高因子;
根据FETCH实验数据可得风浪相干函数为:
通过动静转换原理和风浪相干函数获得叶片上o点和海平面q点的互功率谱:
式中,Cohww′(o,q;f)表示转动域o点和海平面q点的相干函数,即:
结合多域交叉互功率谱,多元相干效应联合功率谱矩阵可表示为:
通过上述方案,在计算中引入风浪相关性,风浪相干函数将海上风机风与波浪相关联,建立了适应于海上风机风场的多元相干效应联合功率谱矩阵,构造的多元相干效应联合功率谱矩阵具备了全域关联性,使得计算更加准确。
进一步,所述步骤七还包括:
多元相干效应联合功率谱矩阵通过改进平方根法分解为上、下三角矩阵和对角矩阵组合的形式:
式中,L(f)q×q为分解后的下三角矩阵,Dq×q为对角矩阵,*表示矩阵的转置;
基于谐波叠加原理构建海上风机多元相干效应全域风速和波高时程:
式中,Δw为频率分辨率,wN为上限截止频率。
通过上述方案,
进一步,所述步骤八还包括:
在同步重构的全域风速和波高时间序列的基础上,通过分段评估的方式构建海上风机多元相干效应全域载荷:
式中,Bwind、Twind和Twave分别表示风或波浪的作用域;a′和a″分别为轴向和切向速度的诱导因子;Ω表示叶片的转速;CN是法向系数,Cd表示风的阻力系数; Aw表示有效迎风面积;ub和ut分别是叶片和塔架上风速的时间序列;CM和CD分别表示波浪的质量系数和阻力系数;ρw是水的密度,Dt表示塔架的直径,zt是塔架分段的长度;和uf分别为流体粒子的加速度和速度;
通过分段求解的方式求解海上风机多元相干效应载荷,其中叶片区空气动力载荷基于叶素动量理论求解,塔架区风载荷和波浪载荷基于载荷作用面积同步求解。
通过上述方案,通过改进平方根法实现了半正定矩阵的稳定分解,耦合全域相干性的多元相干联合功率谱矩阵的分解异常情况被消除。多元相干联合功率谱矩阵被同步分解,保障了重构风速、波高时程和进一步分段求解的多元相干效应全域载荷的同步性和关联性。
本发明的一种海上风机多元相干效应风场载荷构建方法具有以下优点:
1.构建了具有物理表征意义的频域旋转效应风谱解析模型,可以有效表征叶片转动的1p-3p转动倍频,避免了时域旋转采样风速的离散,引入 Gauss-Lejeune方法解决了旋转效应风谱模型震荡求解的收敛性问题。
2.建立叶片转动域不同风速点的动态相干性,首次通过动静转换提出了等效静态相干函数,建立了转动与静止域间的相干性,避免了传统固定式结构上风速模拟空间相干性单一的问题,静止与波浪域、旋转与波浪域相干函数的引入使得风场具有全域相干效应,首创适应于海上风机的、准确完整的全域联合功率谱矩阵。
3.引入改进平方根法对多元相干效应联合功率谱矩阵进行分解,消除了半正定矩阵Cholesky分解异常的问题,多维矩阵被同步分解,进一步保证了风速、波高时间序列以及全域载荷构造的同步性和相干性。
具体实施方式
为便于本领域的技术人员理解本发明,下面说明本发明的具体实施方式。
在本申请的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本申请的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本申请中的具体含义。
本发明提供了一种海上风机多元相干效应风场载荷构建方法,具体步骤如下:
步骤一、检测海上风机的物理量和海上风机所在风场的物理量。确定风速点i和j。检测的物理量具体根据计算要求决定。
步骤二、通过泰勒冻结原理,对海上风机叶片上不同空间位置和不同时刻的风速点进行时空等效。
S200、当风速点i和j分别位于静止位置与转动位置时,将静止位置风速点间的互相关函数与脉动风互功率谱构成傅里叶变换对,旋转位置风速点间的互相关函数与旋转风互功率谱构成傅里叶变换对:
S201、Sij(f)可表示为:
其中Sii(f)和Sjj(f)表示点i和j的脉动风自功率谱,Coh(i,j;f)是i和j之间的相干函数,θij表示相位延迟参数。
步骤三、结合相关性原理与狄拉克函数的性质建立旋转叶片的频域旋转效应风谱。
S300、选定f0作为叶片的转动频率,则相干函数可以表示为展开系数An(i,j;f) 和f0的表达式,即:
S301、基于傅里叶展开系数,旋转效应风谱可进一步展开为:
S302、设置狄拉克函数δ为:
S303、通过旋转效应风谱和狄拉克函数δ获得旋转效应风谱与转动效应模态、源谱平移效应和转动相位延迟的相关函数:
由上式可知,平移后的源谱受转动效应模态和转动相位延迟组成的震荡衰减函数影响来反映叶片转动导致的倍频效应。
上式表示的旋转效应风谱可以表征叶片转动的1p、2p和3p转动倍频。在使用上式进行旋转风速求解时,需要解决震荡积分导致的收敛性问题。故 Gauss-Lejeune方法被用于旋转风速震荡积分的求解。
S306、Gauss-Lejeune积分公式由高斯系数和高斯点所代表的函数之和近似表达为:
式中,tj为高斯点,Aj表示高斯系数,通过上述公式求解旋转风速。
从相关函数的角度构建了适配海上风机叶片的频域旋转效应风谱模型,模型具有物理表征意义,可根据叶片旋转频率重新分配1p~3p频率处的波峰能量,Gauss-Lejeune方法的引入进一步解决了旋转效应风谱的震荡积分收敛问题,有助于为随后多元相干效应风场构建提供理论前提。
步骤四、将海上风机最上面承受的海上环境载荷划分为旋转载荷,将海上风机中间承受的海上环境载荷划分为静风载荷,将海上风机下面承受的海上环境载荷划分为波浪载荷,将海上风机所在风场的风场域对应旋转载荷所在区域划分为转动域,将海上风机所在风场的风场域对应静风载荷所在区域划分为静止域,将海上风机所在风场的风场域对应波浪载荷所在区域划分为风浪作用域。
步骤五、基于叶片转动建立转动域动态相干性函数,对转动坐标系下旋转风速与固定坐标系下脉动风速进行等效变换,建立海上风机叶片转动与塔架静止域的等效静态相干函数。
S500、使用Davenport相干函数表示海上风机叶片转动时叶片上不同点i 和j的动态相干性:
其中,ri和rj分别表示叶片上i和j点到轮毂中心的距离;B1→B1和B1→B2分别表示i和j点位于相同和不同叶片上。
S502、海上风机叶片o点和塔架p点的跨区域风速互功率谱用旋转效应风谱和静止域相干函数表示为:
S503、Cohbt(o,p;f)表示叶片o点和塔架p点的等效静态相干函数,为:
将海上风场划分为四个相干耦合区域,在旋转域建立了叶片转动频率相关的动态相干性函数,通过动静转换原则建立了旋转与静止域的等效静态相干函数,将海上风机波面以上的区域关联起来,并以此构建跨区域风速互功率谱。
步骤六、通过风浪相干函数构建塔架静止与波浪作用域、叶片转动与波浪作用域的互功率谱,结合旋转域的互功率谱、转动与静止域的互功率谱生成多元相干效应联合功率谱矩阵。
S600、选用塔架p点和海平面q点建立塔架静止与波浪作用域的互功率谱:
式中,Cohww(p,q;f)为p点和q点风浪相干函数。
式中,Tp表示谱峰周期且fp=1/Tp,Hs是有效波高,σ是峰形系数,γ是谱峰升高因子。
S602、根据FETCH实验数据可得风浪相干函数为:
S603、通过动静转换原理和风浪相干函数获得叶片上o点和海平面q点的互功率谱:
S604、上述公式中Cohww′(o,q;f)表示转动域o点和海平面q点的相干函数,即:
S605、结合多域交叉互功率谱,多元相干效应联合功率谱矩阵可表示为:
其中多元相干效应联合功率谱矩阵主对角线元素依次表示旋转效应风谱、脉动风谱以及波浪谱,非主对角线元素表示转动、静止以及风浪相干域的交叉互功率谱。静止域风浪相干性的引入使联合功率谱矩阵具备了描述风浪关联的功能。此外,转动域的风浪相干性进一步保证了多域载荷构建的完整性和准确性。
在动态和等效静态相干性的基础上引入风浪相关性,风浪相干函数将海上风机风与波浪相关联,首创适应于海上风机风场的多元相干效应联合功率谱矩阵,构造的多元相干效应联合功率谱矩阵具备了全域关联性。
步骤七、使用改进平方根法对多元相干效应联合功率谱矩阵进行分解,根据谐波叠加原理重构多元相干效应风速和波高时程,计算适应于海上风机的风速幅值和相位参数。
S700、多元相干效应联合功率谱矩阵通过改进平方根法分解为上、下三角矩阵和对角矩阵组合的形式:
式中,L(f)q×q为分解后的下三角矩阵,Dq×q为对角矩阵,*表示矩阵的转置。
S701、基于谐波叠加原理构建海上风机多元相干效应全域风速和波高时程:
S702、通过上述公式计算得:
式中,Δw为频率分辨率,wN为上限截止频率。
步骤八、通过分段评估进一步构建海上风机多元相干效应全域载荷,海上风机多元相干效应全域载荷包括旋转域空气动力载荷、静止域风载荷以及波浪载荷,求解获得海上风机多元相干效应载荷。
S800、在同步重构的全域风速和波高时间序列的基础上,通过分段评估的方式构建海上风机多元相干效应全域载荷,海上风机多元相干效应全域载荷主要包含三个组成部分:旋转域(法向)空气动力载荷、静止域风载荷以及波浪载荷:
式中,Bwind、Twind和Twave分别表示风或波浪的作用域;a′和a″分别为轴向和切向速度的诱导因子;Ω表示叶片的转速;CN是法向系数,Cd表示风的阻力系数;Aw表示有效迎风面积;ub和ut分别是叶片和塔架上风速的时间序列;CM和CD分别表示波浪的质量系数和阻力系数;ρw是水的密度,Dt表示塔架的直径, zt是塔架分段的长度;和uf分别为流体粒子的加速度和速度。
S801、通过分段求解的方式求解海上风机多元相干效应载荷,其中叶片区空气动力载荷基于叶素动量理论求解,塔架区风载荷和波浪载荷基于载荷作用面积同步求解。
本方法解决了由于波浪低频和高频能量低、风低频能量高所造成的联合功率谱矩阵Cholesky分解的非正定问题,通过改进平方根法实现了半正定矩阵的稳定分解,耦合全域相干性的多元相干联合功率谱矩阵的分解异常情况被消除。多元相干联合功率谱矩阵被同步分解,保障了重构风速、波高时程和进一步分段求解的多元相干效应全域载荷的同步性和关联性。
以上所述的本发明实施方式,并不构成对本发明保护范围的限定。任何在本发明的精神和原则之内所作的修改、等同替换和改进等,均应包含在本发明的权利要求保护范围之内。
Claims (8)
1.一种海上风机多元相干效应风场载荷构建方法,其特征在于,包括以下步骤:
一、检测海上风机的物理量和海上风机所在风场的物理量;
二、通过泰勒冻结原理,对海上风机叶片上不同空间位置和不同时刻的风速点进行时空等效;
三、结合相关性原理与狄拉克函数的性质建立旋转叶片的频域旋转效应风谱;
四、将海上风机所在风场的风场域划分为转动域、静止域和风浪作用域;
五、基于叶片转动建立转动域动态相干性函数,对转动坐标系下旋转风速与固定坐标系下脉动风速进行等效变换,建立海上风机叶片转动与塔架静止域的等效静态相干函数;
六、通过风浪相干函数构建塔架静止与波浪作用域、叶片转动与波浪作用域的互功率谱,结合旋转域的互功率谱、转动与静止域的互功率谱生成多元相干效应联合功率谱矩阵;
七、使用改进平方根法对多元相干效应联合功率谱矩阵进行分解,根据谐波叠加原理重构多元相干效应风速和波高时程,计算适应于海上风机的风速幅值和相位参数;
八、通过分段评估进一步构建海上风机多元相干效应全域载荷,海上风机多元相干效应全域载荷包括旋转域空气动力载荷、静止域风载荷以及波浪载荷,求解获得海上风机多元相干效应载荷。
2.根据权利要求1所述的一种海上风机多元相干效应风场载荷构建方法,其特征在于,所述通过泰勒冻结原理,对海上风机叶片上不同空间位置和不同时刻的风速点进行时空等效还包括:
当风速点i和j分别位于静止位置与转动位置时,将静止位置风速点间的互相关函数与脉动风互功率谱构成傅里叶变换对,旋转位置风速点间的互相关函数与旋转风互功率谱构成傅里叶变换对:
其中Sii(f)和Sjj(f)表示点i和j的脉动风自功率谱,Coh(i,j;f)是i和j之间的相干函数,θij表示相位延迟参数。
3.根据权利要求2所述的一种海上风机多元相干效应风场载荷构建方法,其特征在于,所述结合相关性原理与狄拉克函数的性质建立旋转叶片的频域旋转效应风谱还包括:
选定f0作为叶片的转动频率,则相干函数可以表示为展开系数An(i,j;f)和f0的表达式,即:
其中An(i,j;f)为相干函数的傅里叶展开系数,φ0表示叶片的分布角;
基于傅里叶展开系数,旋转效应风谱可进一步展开为:
设置狄拉克函数δ为:
通过旋转效应风谱和狄拉克函数δ获得旋转效应风谱与转动效应模态、源谱平移效应和转动相位延迟的相关函数:
Gauss-Lejeune积分公式由高斯系数和高斯点所代表的函数之和近似表达为:
式中,tj为高斯点,Aj表示高斯系数,通过上述公式求解旋转风速。
4.根据权利要求1所述的一种海上风机多元相干效应风场载荷构建方法,其特征在于,所述将海上风机所在风场的风场域划分为转动域、静止域和风浪作用域还包括:
将海上风机最上面承受的海上环境载荷划分为旋转载荷,将海上风机中间承受的海上环境载荷划分为静风载荷,将海上风机下面承受的海上环境载荷划分为波浪载荷,将海上风机所在风场的风场域对应旋转载荷所在区域划分为转动域,将海上风机所在风场的风场域对应静风载荷所在区域划分为静止域,将海上风机所在风场的风场域对应波浪载荷所在区域划分为风浪作用域。
5.根据权利要求3所述的一种海上风机多元相干效应风场载荷构建方法,其特征在于,所述步骤五还包括:
使用Davenport相干函数表示海上风机叶片转动时叶片上不同点i和j的动态相干性:
其中,ri和rj分别表示叶片上i和j点到轮毂中心的距离;B1→B1和B1→B2分别表示i和j点位于相同和不同叶片上;
海上风机叶片o点和塔架p点的跨区域风速互功率谱用旋转效应风谱和静止域相干函数表示为:
6.根据权利要求5所述的一种海上风机多元相干效应风场载荷构建方法,其特征在于,所述步骤六还包括:
选用塔架p点和海平面q点建立塔架静止与波浪作用域的互功率谱:
式中,Tp表示谱峰周期且fp=1/Tp,Hs是有效波高,σ是峰形系数,γ是谱峰升高因子;
根据FETCH实验数据可得风浪相干函数为:
通过动静转换原理和风浪相干函数获得叶片上o点和海平面q点的互功率谱:
式中,Cohww′(o,q;f)表示转动域o点和海平面q点的相干函数,即:
结合多域交叉互功率谱,多元相干效应联合功率谱矩阵可表示为:
8.根据权利要求7所述的一种海上风机多元相干效应风场载荷构建方法,其特征在于,所述步骤八还包括:
在同步重构的全域风速和波高时间序列的基础上,通过分段评估的方式构建海上风机多元相干效应全域载荷:
式中,Bwind、Twind和Twave分别表示风或波浪的作用域;a′和a″分别为轴向和切向速度的诱导因子;Ω表示叶片的转速;CN是法向系数,Cd表示风的阻力系数;Aw表示有效迎风面积;ub和ut分别是叶片和塔架上风速的时间序列;CM和CD分别表示波浪的质量系数和阻力系数;ρw是水的密度,Dt表示塔架的直径,zt是塔架分段的长度;和uf分别为流体粒子的加速度和速度;
通过分段求解的方式求解海上风机多元相干效应载荷,其中叶片区空气动力载荷基于叶素动量理论求解,塔架区风载荷和波浪载荷基于载荷作用面积同步求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211181493.4A CN115510381B (zh) | 2022-09-27 | 2022-09-27 | 一种海上风机多元相干效应风场载荷构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211181493.4A CN115510381B (zh) | 2022-09-27 | 2022-09-27 | 一种海上风机多元相干效应风场载荷构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115510381A true CN115510381A (zh) | 2022-12-23 |
CN115510381B CN115510381B (zh) | 2023-08-22 |
Family
ID=84505159
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211181493.4A Active CN115510381B (zh) | 2022-09-27 | 2022-09-27 | 一种海上风机多元相干效应风场载荷构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115510381B (zh) |
Citations (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120179376A1 (en) * | 2011-01-11 | 2012-07-12 | Ophir Corporation | Methods And Apparatus For Monitoring Complex Flow Fields For Wind Turbine Applications |
CN204783486U (zh) * | 2015-07-30 | 2015-11-18 | 中国大唐集团新能源股份有限公司 | 一种加长叶片的连接结构 |
CN105426594A (zh) * | 2015-11-06 | 2016-03-23 | 西南交通大学 | 一种基于时空场和条件插值的平稳均质风场快速模拟方法 |
CN106762394A (zh) * | 2017-04-01 | 2017-05-31 | 孟英志 | 一种带加长装置的风机风轮或叶片及风力机 |
CN107247687A (zh) * | 2017-05-18 | 2017-10-13 | 西南交通大学 | 一种基于特征正交分解的非平稳随机过程快速模拟方法 |
CN107701376A (zh) * | 2017-10-17 | 2018-02-16 | 西南交通大学 | 风机单叶片安装桨距调节方法 |
EP3324043A1 (en) * | 2016-11-21 | 2018-05-23 | LM WP Patent Holding A/S | Method for controlling a floating offshore wind turbine, wind turbine control system and floating offshore wind turbine |
CN108229043A (zh) * | 2018-01-12 | 2018-06-29 | 中国海洋大学 | 考虑涡激效应的深海spar型浮式风机疲劳损伤分析方法 |
CN108694277A (zh) * | 2018-05-06 | 2018-10-23 | 北京工业大学 | 基于多平台联合仿真的风力发电机动力响应分析方法 |
CN109653963A (zh) * | 2018-12-26 | 2019-04-19 | 三重能有限公司 | 叶片动平衡评估方法、装置及系统 |
CN110414102A (zh) * | 2019-07-16 | 2019-11-05 | 大连理工大学 | 海上风力发电结构疲劳分析方法 |
CN110847210A (zh) * | 2019-11-18 | 2020-02-28 | 中国海洋大学 | 一种海上风机基础结构以及轻量化设计方法 |
CN111396246A (zh) * | 2019-11-27 | 2020-07-10 | 浙江运达风电股份有限公司 | 基于叶轮等效风速修正的激光雷达辅助控制方法 |
EP3705719A1 (en) * | 2019-03-05 | 2020-09-09 | MHI Vestas Offshore Wind A/S | Designing towers of offshore wind turbines |
CN111796244A (zh) * | 2020-06-09 | 2020-10-20 | 哈尔滨工业大学 | 一种海上风电场对地波超视距雷达探测效能影响的评估方法 |
CN111985018A (zh) * | 2020-03-31 | 2020-11-24 | 重庆科技学院 | 基于惯性力法和塔线分离法考虑塔线耦合影响的超高大跨越塔、线设计风载荷的计算方法 |
CN112727710A (zh) * | 2020-12-15 | 2021-04-30 | 北京天泽智云科技有限公司 | 一种基于音频信号的风场落雷密度统计方法及系统 |
CN113392547A (zh) * | 2021-05-31 | 2021-09-14 | 东南大学 | 一种考虑空间相关性的频域分布随机动载荷识别方法 |
CN113468692A (zh) * | 2021-07-19 | 2021-10-01 | 大连理工大学 | 一种基于延迟效应的三维风场高效模拟方法 |
CN113609622A (zh) * | 2021-08-20 | 2021-11-05 | 浙江大学 | 风力发电机组的塔架载荷建模方法和装置 |
CN113685314A (zh) * | 2021-08-24 | 2021-11-23 | 浙江大学 | 变桨控制方法、系统和可读存储介质 |
CN113822565A (zh) * | 2021-09-16 | 2021-12-21 | 中国电建集团华东勘测设计研究院有限公司 | 一种风机监测数据时频特征分级细化分析的方法 |
CN114117663A (zh) * | 2021-11-09 | 2022-03-01 | 明阳智慧能源集团股份公司 | 一种海上固定式风机一体化载荷分析方法及系统 |
CN114781074A (zh) * | 2022-06-20 | 2022-07-22 | 交通运输部天津水运工程科学研究所 | 一种复杂环境下船舶仿真模拟及风险评估方法 |
-
2022
- 2022-09-27 CN CN202211181493.4A patent/CN115510381B/zh active Active
Patent Citations (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120179376A1 (en) * | 2011-01-11 | 2012-07-12 | Ophir Corporation | Methods And Apparatus For Monitoring Complex Flow Fields For Wind Turbine Applications |
CN204783486U (zh) * | 2015-07-30 | 2015-11-18 | 中国大唐集团新能源股份有限公司 | 一种加长叶片的连接结构 |
CN105426594A (zh) * | 2015-11-06 | 2016-03-23 | 西南交通大学 | 一种基于时空场和条件插值的平稳均质风场快速模拟方法 |
EP3324043A1 (en) * | 2016-11-21 | 2018-05-23 | LM WP Patent Holding A/S | Method for controlling a floating offshore wind turbine, wind turbine control system and floating offshore wind turbine |
CN106762394A (zh) * | 2017-04-01 | 2017-05-31 | 孟英志 | 一种带加长装置的风机风轮或叶片及风力机 |
CN107247687A (zh) * | 2017-05-18 | 2017-10-13 | 西南交通大学 | 一种基于特征正交分解的非平稳随机过程快速模拟方法 |
CN107701376A (zh) * | 2017-10-17 | 2018-02-16 | 西南交通大学 | 风机单叶片安装桨距调节方法 |
CN108229043A (zh) * | 2018-01-12 | 2018-06-29 | 中国海洋大学 | 考虑涡激效应的深海spar型浮式风机疲劳损伤分析方法 |
CN108694277A (zh) * | 2018-05-06 | 2018-10-23 | 北京工业大学 | 基于多平台联合仿真的风力发电机动力响应分析方法 |
CN109653963A (zh) * | 2018-12-26 | 2019-04-19 | 三重能有限公司 | 叶片动平衡评估方法、装置及系统 |
EP3705719A1 (en) * | 2019-03-05 | 2020-09-09 | MHI Vestas Offshore Wind A/S | Designing towers of offshore wind turbines |
CN110414102A (zh) * | 2019-07-16 | 2019-11-05 | 大连理工大学 | 海上风力发电结构疲劳分析方法 |
CN110847210A (zh) * | 2019-11-18 | 2020-02-28 | 中国海洋大学 | 一种海上风机基础结构以及轻量化设计方法 |
CN111396246A (zh) * | 2019-11-27 | 2020-07-10 | 浙江运达风电股份有限公司 | 基于叶轮等效风速修正的激光雷达辅助控制方法 |
CN111985018A (zh) * | 2020-03-31 | 2020-11-24 | 重庆科技学院 | 基于惯性力法和塔线分离法考虑塔线耦合影响的超高大跨越塔、线设计风载荷的计算方法 |
CN111796244A (zh) * | 2020-06-09 | 2020-10-20 | 哈尔滨工业大学 | 一种海上风电场对地波超视距雷达探测效能影响的评估方法 |
CN112727710A (zh) * | 2020-12-15 | 2021-04-30 | 北京天泽智云科技有限公司 | 一种基于音频信号的风场落雷密度统计方法及系统 |
CN113392547A (zh) * | 2021-05-31 | 2021-09-14 | 东南大学 | 一种考虑空间相关性的频域分布随机动载荷识别方法 |
CN113468692A (zh) * | 2021-07-19 | 2021-10-01 | 大连理工大学 | 一种基于延迟效应的三维风场高效模拟方法 |
CN113609622A (zh) * | 2021-08-20 | 2021-11-05 | 浙江大学 | 风力发电机组的塔架载荷建模方法和装置 |
CN113685314A (zh) * | 2021-08-24 | 2021-11-23 | 浙江大学 | 变桨控制方法、系统和可读存储介质 |
CN113822565A (zh) * | 2021-09-16 | 2021-12-21 | 中国电建集团华东勘测设计研究院有限公司 | 一种风机监测数据时频特征分级细化分析的方法 |
CN114117663A (zh) * | 2021-11-09 | 2022-03-01 | 明阳智慧能源集团股份公司 | 一种海上固定式风机一体化载荷分析方法及系统 |
CN114781074A (zh) * | 2022-06-20 | 2022-07-22 | 交通运输部天津水运工程科学研究所 | 一种复杂环境下船舶仿真模拟及风险评估方法 |
Non-Patent Citations (1)
Title |
---|
孟繁星: "风浪联合作用下的海上风机支撑结构时域疲劳分析" * |
Also Published As
Publication number | Publication date |
---|---|
CN115510381B (zh) | 2023-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Coulling et al. | Validation of a FAST semi-submersible floating wind turbine numerical model with DeepCwind test data | |
Veers | Modeling stochastic wind loads on vertical axis wind turbines | |
Sasaki | Numerical variational analysis formulated under the constraints as determined by longwave equations and a low-pass filter | |
CN104408295B (zh) | 一种大跨桥梁下部结构风‑浪耦合作用荷载数值模拟方法 | |
CN105354377A (zh) | 一种确定输电杆塔的脉动风振荷载的方法 | |
CN110794459A (zh) | 一种海底近断层地震动的拟合方法 | |
Abhinav et al. | Nonlinear dynamical behaviour of jacket supported offshore wind turbines in loose sand | |
CN103617308B (zh) | 一种风电场频域等效模型的构建方法 | |
CN104516771A (zh) | 一种非平稳随机过程高效模拟方法 | |
Halfpenny | Dynamic analysis of both on and offshore wind turbines in the frequency domain | |
Selvam | Individual pitch control for large scale wind turbines | |
CN110059286A (zh) | 一种基于fft的结构非平稳响应高效分析方法 | |
CN115114841A (zh) | 基于U-net框架的边缘海高时空分辨率次表层温度场重构方法 | |
CN115510381A (zh) | 一种海上风机多元相干效应风场载荷构建方法 | |
CN114813123A (zh) | 一种基于pso-vmd-mckd的滚动轴承微弱故障诊断方法 | |
CN112785052A (zh) | 基于粒子滤波算法的风速风向预测方法 | |
Liu | A CFD study of fluid-structure interaction problems for floating offshore wind turbines | |
Pegalajar-Jurado et al. | Accelerated hydrodynamic analysis for spar buoys with second-order wave excitation | |
CN112918635B (zh) | 一种基于电流信号的海上浮式风电机组基础运动监测方法 | |
Chen et al. | Extraction of natural frequencies and mode shapes of rotating beams by variational iteration method | |
Jiang et al. | Rolling element bearing fault feature extraction using an optimal chirplet | |
Li et al. | Structural response of offshore monopile foundations to ocean waves | |
Mu et al. | Simulation and Control Applications Based on PM Wave Spectra | |
CN111400968B (zh) | 一种曲线坐标系下可压缩两相流界面条件构建方法 | |
CN109782338B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |