CN105508147B - 风电叶片单点疲劳加载试验弯矩匹配方法 - Google Patents
风电叶片单点疲劳加载试验弯矩匹配方法 Download PDFInfo
- Publication number
- CN105508147B CN105508147B CN201511004114.4A CN201511004114A CN105508147B CN 105508147 B CN105508147 B CN 105508147B CN 201511004114 A CN201511004114 A CN 201511004114A CN 105508147 B CN105508147 B CN 105508147B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- blade
- moment
- balancing weight
- 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.)
- Expired - Fee Related
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2260/00—Function
- F05B2260/83—Testing, e.g. methods, components or tools therefor
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2260/00—Function
- F05B2260/84—Modelling or simulation
-
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/30—Wind power
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/20—Climate change mitigation technologies for sector-wide applications using renewable energy
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Wind Motors (AREA)
Abstract
本发明属于风机叶片的单点疲劳加载试验***,具体涉及一种风电叶片单点疲劳加载试验弯矩匹配方法,根据等效替代原则将叶片沿翼向离散成n个离散部分,确立实际总弯矩值模型,最后结合二分法原则,从叶片尖部开始,针对某一截面,将理论弯矩值与实际总弯矩值之差的最小值作为目标函数,其他截面的误差作为不等式约束条件,设计了基于目标函数及约束条件的优化求解算法,并通过建立数学模型优化出所添加配重块的个数、质量和位置且保证所添加的配重块个数最少。通过本发明能够在配重块个数添加最少的前提下使叶片的实际试验弯矩尽可能与理论弯矩相匹配,提高叶片的疲劳加载试验精度,其精度可严格达到工业要求的δ(目前通常为7%)的误差要求。
Description
技术领域
本发明涉及一种风电叶片单点疲劳加载试验弯矩匹配方法,尤其适用于兆瓦级风电叶片的疲劳加载试验,属于风机叶片的单点疲劳加载试验***。
背景技术
目前单点疲劳加载试验方法是国际上风电叶片疲劳测试的主流方法之一。单点疲劳加载试验方法通常在叶片翼向约70%处施加单点激励使其与叶片共振来完成疲劳测试。根据IEC 61400-23Full-Scale Structural Testing of Wind Turbine Blade标准,由该单点激励所引起的叶片各个端截面的实际弯矩应尽可能与叶片设计方提供的理论弯矩相匹配。国内外通常采用在叶片表面添加多个配重块来保证各截面的弯矩匹配误差δ控制在一定的误差范围内(目前该误差通常取值为7%)。但是由于缺乏有效的弯矩匹配方法,多数叶片测试厂家仅通过经验或简单计算,在叶片翼向的合适位置仅添加一个配重块来实现实际弯矩和理论弯矩的匹配,匹配误差较大,某些截面的相对误差远远超过7%,由此得出的疲劳测试数据精度不高,在一定程度上造成了单点疲劳加载试验测试结果失真,很难满足高精度的叶片疲劳测试要求。或者如张磊安的论文《MW级风电叶片加载***关键技术研究.》中,虽给出了弯矩计算模型,但模型过于简单且后续优化步骤含糊,没有给出具体的实现方法,由此得出的弯矩分布数据后人无法验证,可靠性不高。根据某叶片生产厂商提供的数据:设计服役期为20年的某型号叶片,虽然通过了单点疲劳加载试验,但实际寿命远达不到设计寿命。随着风力机逐渐向兆瓦级大功率方向发展,叶片尺寸随之增大,对叶片的强度和刚度的要求更加苛刻,高精度的叶片疲劳加载试验必将是叶片检测领域的研究重点之一。因此,提出一种有效的实际-理论弯矩匹配方法,可以提高试验中实际弯矩的分布精度,从而提高叶片的疲劳加载试验精度。
发明内容
根据以上现有技术中的不足,本发明要解决的技术问题是:提供一种能够使理论弯矩和实际弯矩两者之间的误差控制在δ以内、有效提高叶片的实际—理论弯矩匹配的风电叶片单点疲劳加载试验弯矩匹配方法。
本发明所述的风电叶片单点疲劳加载试验弯矩匹配方法,包括以下步骤:
(1)根据等效替代原则将叶片沿翼向分割为n个离散部分,得到(n+1)个截面,试验中通过疲劳加载驱动装置带动叶片振动,同时产生激振力;
(2)建立模型:
建立仅考虑叶片自重的实际弯矩模型,该模型表示为
T1k为该模型中叶片某一截面k(k=1,2,…,n+1)处的实际弯矩值;
建立添加配重块后、考虑驱动装置重量与配重块重量的实际弯矩模型,该模型表示为
T2k为该模型中叶片某一截面k(k=1,2,…,n+1)处的实际弯矩值;
由上述得到实际总弯矩模型,该模型表示为
Tk=T1k+T2k;[3]
Tk为该模型中叶片某一截面k(k=1,2,…,n+1)处的实际总弯矩值;
其中,i为截面编号,j为截面k右边的配重块编号(截面k右边如果没有配重块,j就不用编号了,即j不存在,则T2k=0),N为所有添加的配重块总数,p为截面k右边的配重块的总数(p≤N),ρi为各离散部分的线质量密度,bi为各离散部分的长度,Lki为端截面k与第i个离散部分ρi重心处的距离,f为叶片激振频率,yi为ρi所属离散部分重心处的振幅,g为重力加速度,tk为端截面k到叶根的距离,ymj为配重块mj重心处的振幅值,rk为疲劳加载驱动装置到端截面k的距离,M为疲劳加载驱动装置的等效质量,yM为疲劳加载驱动装置处的振幅值,以上各参数均为可测出的已知值;mj为添加的配重块质量,xj为添加的配重块距叶片根部的距离,mj与xj为未知值;
(3)要计算Tk,所添加的各个配重块的质量mj以及各个配重块到叶根的距离xj(即配重块的添加位置)均为未知值,而已知条件为:根据疲劳测试要求,对于叶片的任一截面k(1≤k≤n+1),其上的实际总弯矩值Tk与理论弯矩值Tk'的误差需要控制在一定范围内,即:
|T′k-Tk|/Tk≤δ(k=1,2…,n+1);[4]
其中,δ为原始设定值;
按照公式[4]进行优化求解,最终得出所添加配重块的个数和质量以及位置且所添加的配重块个数最少。在优化求解时,可以采用常规的数学软件,如Matlab、Lindo、Lingo等。
本发明中优选如下优化求解的方法,具体方法为:
1)计算出未添加配重块之前实际总弯矩值与理论弯矩值相对误差最大的截面,设该截面距离叶片根部的距离为从叶片尖部开始添加,先添加一个配重块(即令N=1),所添加的配重块距离该截面越远,所需要的配重块的质量越小,但最远距离不能超过满足误差要求的截面;找出实际总弯矩值与理论弯矩值相对误差最小的截面,设该截面距离叶片根部的距离为
①采用二分法确定配重块的初始位置:
根据二分法原则,令取x1 1为x1的初始值,任选一个截面a,a=1,2,…,n+1,将其上的理论弯矩值Ta'与实际总弯矩值Ta之差的最小值作为目标函数,将其他截面误差的差值作为该目标函数的不等式约束条件,建立目标优化数学模型如下:
minG=T′a-Ta;[5]
其中,公式[5]和公式[6]中的j=1;
求解m1的值,若m1有解,则令x1 1=x1,得到添加一个配重块时x1、m1的最优值,完成优化;
②若无解,根据二分法原则,再次令取x1 2为x1的初始值,任选一个截面c,截面c与①中的截面a可以是同一截面,也可以是不同截面,将其上的理论弯矩值Tc'与实际总弯矩值Tc之差的最小值作为目标函数,将其他截面的误差作为该目标函数的不等式约束条件,如①建立目标优化数学模型(将相应公式中的下标a改为c,p为截面c右侧配重块的个数),求解m1的值,若m1有解,则令x1 2=x1,得到添加一个配重块时x1、m1的最优值,完成优化;
③若无解,则再根据二分法原则,令如此循环,直到得到添加一个配重块时x1、m1的最优值;
④若x1 n趋近于(即ε根据误差大小给定)时仍然无解,即仍然不满足误差要求,则需要添加第二个配重块,并计算其添加位置,此时取x1 1、x1 2……x1 n中实际总弯矩值与理论弯矩值相对误差最小的位置点作为添加第一个配重块的位置初始值,并取该位置点在上述目标优化数学模型中求解的m1作为第一个配重块的质量,之后添加第二个配重块;
2)计算出添加完第一个配重块之后实际总弯矩值与理论弯矩值相对误差最大的截面,设该截面距离叶片根部的距离为同时找出实际总弯矩值与理论弯矩值相对误差最小的截面,设该截面距离叶片根部的距离为根据二分法原则,令然后取x2 1为x2的初始值,采用如同1)中的步骤①②③——由于采用的配重块个数为2,因此公式[5]和公式[6]中的N=2——得到添加两个配重块后的x1、m1、x2和m2的最优值;若x2 n趋近于(即且ε根据误差大小给定)时仍然无解,即仍然不满足误差要求,则需要添加第三个配重块并计算其添加位置,此时,取x2 1、x2 2……x2 n中实际总弯矩值与理论弯矩值相对误差最小的位置点作为第二个配重块的添加位置,将这两个配重块的添加位置代入上述建立的目标优化数学模型中(此时N=2),求解的m1作为第一个配重块的质量、m2作为第二个配重块的质量,之后添加第三个配重块;
3)添加第三个配重块的优化求解步骤如同2),若添加三个配重块仍无解,则需要添加第四个,……,直至添加第N’个配重块有解,则取已经求解出的N个配重块的添加位置点x1、x2……xN’,此时在上述建立的目标优化数学模型中(令N=N’)求解出各添加位置点所对应的配重块的质量;
按照上述优化求解的方法,即可得到满足公式[4]的所添加的配重块的个数、质量和位置,且所添加的配重块个数最少。
本发明与现有技术相比所具有的有益效果是:
1、本发明对叶片任一端截面,计算其实际弯矩分布的模型精度较高;
2、本发明针对所添加配重块的数量、位置和质量,给出了具体的计算和优化方法,并且满足所添加的配重块的数量最少;
3、本发明弯矩匹配的精度高,叶片任一端截面的理论弯矩和实际弯矩误差可严格控制在δ(目前通常为7%)以内。
附图说明
图1是仅考虑叶片自重的弯矩分布分析图;
图2是添加配重块后、考虑驱动装置重量与配重块重量的弯矩分布分析图;
图3是本实施例中某型叶片的理论弯矩值曲线图;
图4是本实施例中某型叶片的等效线质量密度的曲线图;
图5是本实施例中某型叶片的各截面的振幅值的曲线图;
图6是本实施例中某型叶片在激振力与自重共同作用下的各截面处弯矩值曲线图;
图7是本实施例中某型叶片的各截面的实际总弯矩值与理论弯矩值的比较图;
图8是本实施例中某型叶片各截面的实际总弯矩值与理论弯矩值匹配后的误差效果图。
具体实施方式
下面结合具体实施例对本发明做进一步描述:
某叶片测试中心提供了长度为40.3m的某型叶片,并根据等效替代原则将该叶片沿翼向离散成22个离散部分,得到23个截面,i为截面编号,每个离散部分的长度bi=1.5m,表1给出了该叶片疲劳加载试验中各截面的理论弯矩值T′k、每个离散部分的线质量密度ρi以及通过激光测试仪测出的各离散部分重心处的振幅yi;并通过疲劳加载驱动装置带动各离散部分旋转产生激振力,该疲劳加载驱动装置的质量M=700Kg,其重心处的振幅yM=0.4,其带动该叶片振动的振动频率f=0.78Hz。本实施例中采用数学软件Matlab进行优化求解。
表1
本实施例中的某型叶片的单点疲劳加载试验弯矩匹配方法如下所示:
(1)在疲劳加载驱动装置产生的激振力下,仅考虑叶片自重,将表1数据带入公式[1]:
对离散而形成的23个截面,其上的实际弯矩值可分别计算得出,计算过程如下所示:
(2)添加配重块后,考虑疲劳加载驱动装置的重量与配重块的重量,则所有截面的实际弯矩值计算如下:
添加的配重块质量mj和位置xj为未知变量,用矩阵表示为:
式中:N为添加配重块的数量。
将表1中的数据带入公式[2]:
其中:p为截面k右边的配重块的总数,p≤N;j为截面k右边的配重块编号;若截面k右边没有配重块,即j不存在,则j无需编号,T2k=0;
其上的实际弯矩值可分别计算得出,计算结果如下:
此时,[T21,T22,…,T223]中含未知量mj、xj。
(3)计算所有截面的实际总弯矩值,根据公式[3]:
Tk=T1k+T2k
可得:
此时,[T1,T2,…,T23]中含未知量mj、xj。
(4)根据疲劳测试要求,对于该叶片的任一截面k(1≤k≤23),其上的实际总弯矩值Tk与理论弯矩值Tk'的相对误差δ需要控制7%的范围内,即需要满足公式[4]:
|T′k-Tk|/Tk≤δ
即应满足:
按照公式[4]将上述带有未知量mj、xj的目标函数在Matlab中优化,其中优化步骤如下:
1)首先计算出未添加配重块之前实际总弯矩值与理论弯矩值相对误差最大的截面位置为所有截面都不满足误差要求,其中实际总弯矩值与理论弯矩值相对误差最小的截面位置为令根据二分法原则确定第一个配重块的初始值(设定时,ε=10-2,x1 n趋近于):
①首先取x1 1为x1的初始值,为方便计算,选截面k=1的实际总弯矩值与理论弯矩值之差的最小值作为目标函数,其他截面的误差的差值作为该目标函数的不等式约束条件,在Matlab中建立目标优化数学模型如下:
minG=T′1-T1 [5]
注:公式[5]、公式[6]中的N=1,p为截面k右边配重块的个数(p≤N);
计算结果显示无解,则继续步骤②。
②根据二分法原则重新配重块的初始位置,则x1 2=30.23m,如步骤①在Matlab中建立目标优化数学模型求解m1的值,结果仍显示无解,则继续步骤③。
③如同步骤②重新确定配重块的初始位置依次优化下去,结果均显示无解,则继续步骤④。
④当(即x1 n趋近于)时仍然不满足误差要求,则需要添加第二个配重块并计算其添加位置,此时,取x1 1、x1 2……x1 14中实际总弯矩值与理论弯矩值相对误差最小的位置点(x1 5=38.59m)作为添加第一个配重块的位置初始值,m1=417kg作为第一个配重块的质量,继续步骤2),添加第二个配重块。
2)计算出添加完第一个配重块之后实际总弯矩值与理论弯矩值相对误差最大的截面,其距离叶片根部的距离同时找出添加完第一个配重块之后实际总弯矩值与理论弯矩值相对误差最小的截面,其距离叶片根部的距离则 根据二分法原则,取x2 1为x2的初始值,采用如同1)中的步骤①②③(注:公式[5]、公式[6]中的N=2),结果均显示无解。当时仍然不满足误差要求,则需要添加第三个配重块并计算其添加位置。此时,取x2 1、x2 2……x2 10中实际总弯矩值与理论弯矩值相对误差最小的位置点x2 6=32.48m作为第二个配重块的添加位置,并取上述两个配重块的位置点x1=38.59m、x2=32.48m在Matlab中求解的m1=349.4kg作为第一个配重块的质量,m2=377kg作为第二个配重块的质量,继续3),添加第三个配重块。
3)添加第三个配重块的优化求解步骤如同2),求解出第三个配重块的添加位置点x3=21.5m,但仍不满足误差要求,则需要添加第四个配重块,步骤同2)求解出第四个配重块的添加位置点x4=10.5m,此时可求解出四个添加位置点所对应的配重块的质量如下:
x1=38.59m,m1=330kg
x2=32.48m,m2=336kg
x3=21.5m,m3=170kg
x4=10.5m,m4=90kg
即配重块个数最少为四个可以满足实际总弯矩值与理论弯矩值的相对误差在7%以内。
将以上数据带入公式[3],计算出添加配重块后各截面的实际总弯矩值,并将该实际总弯矩值与理论弯矩值进行比较,得出二者的误差百分比,具体如表2所示:
表2
由表2也可以看出,按照本发明所述叶片单点疲劳加载试验弯矩匹配方法所得的弯矩匹配的精度高,叶片任一截面的理论弯矩值和实际总弯矩值的误差可严格控制在所要求的7%以内。
Claims (1)
1.一种风电叶片单点疲劳加载试验弯矩匹配方法,其特征在于:包括以下步骤:
(1)根据等效替代原则将叶片沿翼向分割为n个离散部分,得到n+1个截面,试验中通过疲劳加载驱动装置带动叶片振动,同时产生激振力;
(2)建立模型:
建立仅考虑叶片自重的实际弯矩模型,该仅考虑叶片自重的实际弯矩模型表示为<math display = 'block'>
<mrow>
<msub>
<mi>T1</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mrow>
<msubsup>
<mo>&sum;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</msubsup>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>L</mi>
<mi>ki</mi>
</msub>
</mrow>
</mrow>
<mfenced open = '[' close = ']'>
<mrow>
<msup>
<mfenced open = '(' close = ')'>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
</mrow>
</mfenced>
<mn>2</mn>
</msup>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>&plus;</mo>
<mi>g</mi>
</mrow>
</mfenced>
</mrow>
</math>,T1k为该仅考虑叶片自重的实际弯矩模型中叶
片某一截面k处的实际弯矩值,其中k=1,2,…,n+1;
建立添加配重块后、考虑驱动装置重量与配重块重量的实际弯矩模型,该添加配重块后、考虑驱动装置重量与配重块重量的实际弯矩模型表示为
<math display = 'block'>
<mrow>
<msub>
<mi>T2</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mrow>
<msubsup>
<mo>&sum;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</msubsup>
<mrow>
<msub>
<mi>m</mi>
<mi>j</mi>
</msub>
<mfenced open = '(' close = ')'>
<mrow>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>&minus;</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
</mrow>
</mfenced>
</mrow>
</mrow>
<mfenced open = '[' close = ']'>
<mrow>
<msup>
<mfenced open = '(' close = ')'>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
</mrow>
</mfenced>
<mn>2</mn>
</msup>
<msub>
<mi>y</mi>
<mi>mj</mi>
</msub>
<mo>&plus;</mo>
<mi>g</mi>
</mrow>
</mfenced>
<mo>+</mo>
<mi>M</mi>
<msub>
<mi>r</mi>
<mi>k</mi>
</msub>
<mfenced open = '[' close = ']'>
<mrow>
<msup>
<mfenced open = '(' close = ')'>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
</mrow>
</mfenced>
<mn>2</mn>
</msup>
<msub>
<mi>y</mi>
<mi>M</mi>
</msub>
<mo>&plus;</mo>
<mi>g</mi>
</mrow>
</mfenced>
</mrow>
</math>,
T2k为该添加配重块后、考虑驱动装置重量与配重块重量的实际弯矩模型中叶片某一截面k处的实际弯矩值,其中,k=1,2,…,n+1;
实际总弯矩模型表示为Tk=T1k+T2k,Tk为该实际总弯矩模型中叶片某一截面k处的实际总弯矩值,其中k=1,2,…,n+1;
其中,i为截面编号,j为截面k右边的配重块编号,N为所有添加的配重块总数,p为截面k右边的配重块的总数,其中p≤N,ρi为各离散部分的线质量密度,bi为各离散部分的长度,Lki为端截面k与第i个离散部分ρi重心处的距离,f为叶片激振频率,yi为ρi所属离散部分重心处的振幅,g为重力加速度,tk为端截面k到叶根的距离,ymj为配重块mj重心处的振幅值,rk为疲劳加载驱动装置到端截面k的距离,M为疲劳加载驱动装置的等效质量,yM为疲劳加载驱动装置处的振幅值,以上ρi、bi、Lki、f、yi、tk、ymj、rk、M、yM均为可测出的已知值; mj为添加的配重块质量,xj为添加的配重块距叶片根部的距离,mj与xj为未知值;
(3)根据<math display = 'block'>
<mrow>
<mrow>
<apply>
<csymbol>
<msup>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>&#8242;</mo>
</msup>
<mo>&minus;</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
</csymbol>
</apply>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
</mrow>
<mi>&delta;</mi>
</mrow>
</math>进行优化求解,其中,<math display = 'block'>
<mrow>
<msup>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>&#8242;</mo>
</msup>
</mrow>
</math>表示理论弯矩值,k=1,2,…,n
+1,δ为原始设定值,最终得出所添加配重块的个数和质量以及位置且所添加的配重块个数
最少。 1
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511004114.4A CN105508147B (zh) | 2015-12-28 | 2015-12-28 | 风电叶片单点疲劳加载试验弯矩匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511004114.4A CN105508147B (zh) | 2015-12-28 | 2015-12-28 | 风电叶片单点疲劳加载试验弯矩匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105508147A CN105508147A (zh) | 2016-04-20 |
CN105508147B true CN105508147B (zh) | 2017-12-08 |
Family
ID=55716377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511004114.4A Expired - Fee Related CN105508147B (zh) | 2015-12-28 | 2015-12-28 | 风电叶片单点疲劳加载试验弯矩匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105508147B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102016117191A1 (de) * | 2016-09-13 | 2018-03-15 | fos4X GmbH | Verfahren und Vorrichtung zur Ermittlung von Belastungen auf einen Turm einer Windenergieanlage |
CN110296053B (zh) * | 2018-03-23 | 2020-06-05 | 北京金风慧能技术有限公司 | 风力发电机叶片的疲劳损伤监测方法和装置 |
CN110232238B (zh) * | 2019-06-11 | 2023-02-07 | 南昌工程学院 | 一种加载器布局方法、装置、计算机设备及存储介质 |
CN113447376B (zh) * | 2021-07-05 | 2023-03-28 | 山东理工大学 | 一种用于风电叶片双轴疲劳试验的弯矩匹配优化方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102004838A (zh) * | 2010-12-17 | 2011-04-06 | 中国航天空气动力技术研究院 | 一种基于有限差分法的风力机叶片结构确定方法 |
CN202768249U (zh) * | 2012-08-03 | 2013-03-06 | 国电联合动力技术有限公司 | 一种基于气动转矩计算模型的风电机组控制*** |
CN104732060A (zh) * | 2015-01-19 | 2015-06-24 | 湖南科技大学 | 一种大型风电机组叶片多重载荷在线识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014121800A1 (en) * | 2013-02-08 | 2014-08-14 | Vestas Wind Systems A/S | Model based controller for a wind turbine generator |
-
2015
- 2015-12-28 CN CN201511004114.4A patent/CN105508147B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102004838A (zh) * | 2010-12-17 | 2011-04-06 | 中国航天空气动力技术研究院 | 一种基于有限差分法的风力机叶片结构确定方法 |
CN202768249U (zh) * | 2012-08-03 | 2013-03-06 | 国电联合动力技术有限公司 | 一种基于气动转矩计算模型的风电机组控制*** |
CN104732060A (zh) * | 2015-01-19 | 2015-06-24 | 湖南科技大学 | 一种大型风电机组叶片多重载荷在线识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105508147A (zh) | 2016-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105508147B (zh) | 风电叶片单点疲劳加载试验弯矩匹配方法 | |
EP3791060B1 (en) | Wind turbine control method | |
CN102880764B (zh) | 一种风力发电机组前机架强度的计算方法 | |
CN103541862B (zh) | 降低风力涡轮机极端负载的桨距调节方法和*** | |
CN103970965A (zh) | 燃气涡轮发动机加速寿命试验试车方法 | |
CN102866637B (zh) | 一种基于二次降阶的带操纵面机翼非定常气动力模拟方法 | |
CN114638046B (zh) | 一种铁路桥墩数字孪生变截面仿真计算方法 | |
CN103678762A (zh) | 优化的复合材料机翼气动弹性风洞模型的缩比建模方法 | |
CN106777549B (zh) | 一种面向荷载试验的桥梁多尺度有限元模拟方法 | |
CN107633105A (zh) | 一种基于改进混合蛙跳算法的四旋翼无人机参数辨识方法 | |
CN105117539B (zh) | 风力机叶片模态频率及其双峰高斯分布拟合方法 | |
Campagnolo et al. | Design, manufacturing and characterization of aero-elastically scaled wind turbine blades for testing active and passive load alleviation techniques within a ABL wind tunnel | |
CN106446424A (zh) | 一种非定常气动力模型参数预测方法 | |
CN106286151B (zh) | 一种风电机组低速轴扭转载荷监测方法及载荷分析方法 | |
Bottasso et al. | Wind tunnel testing of wind turbines and farms | |
CN111680823A (zh) | 一种风向信息预测方法及*** | |
CN106097157A (zh) | 一种同步发电机励磁***模型参数辨识优化方法 | |
CN105956318B (zh) | 基于改进***h-k聚类方法的风电场机群划分方法 | |
CN107704679A (zh) | 直驱式发电机模型的修正设备和方法 | |
CN104318005A (zh) | 一种基于Virtools的多元数据融合虚拟试验显示*** | |
JP2021049841A (ja) | 制御装置、学習装置、制御方法、学習方法、及びプログラム | |
CN106786675A (zh) | 一种电力***稳定器及其实现方法 | |
CN104408223B (zh) | 一种风电机组的跟踪优化控制方法 | |
CN101794338A (zh) | 基于结构模态试验的矩阵型动力学模型修正方法 | |
CN116305589A (zh) | 一种直升机桨叶结构降阶分析方法、***、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171208 Termination date: 20191228 |
|
CF01 | Termination of patent right due to non-payment of annual fee |