CN104093971B - 风轮机叶片控制方法 - Google Patents

风轮机叶片控制方法 Download PDF

Info

Publication number
CN104093971B
CN104093971B CN201280061666.7A CN201280061666A CN104093971B CN 104093971 B CN104093971 B CN 104093971B CN 201280061666 A CN201280061666 A CN 201280061666A CN 104093971 B CN104093971 B CN 104093971B
Authority
CN
China
Prior art keywords
blade
mode
shape
lifting device
active lifting
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.)
Active
Application number
CN201280061666.7A
Other languages
English (en)
Other versions
CN104093971A (zh
Inventor
P.比克
C.F.安德森
M.O.斯洛特
C.斯科拜
S.B.博杰森
M.拉夫恩
M.克里特加亚尔德
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
LM Wind Power AS
Original Assignee
LM Wind Power AS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by LM Wind Power AS filed Critical LM Wind Power AS
Publication of CN104093971A publication Critical patent/CN104093971A/zh
Application granted granted Critical
Publication of CN104093971B publication Critical patent/CN104093971B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D7/00—Controlling wind motors 
    • F03D7/02—Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/022—Adjusting aerodynamic properties of the blades
    • F03D7/0232—Adjusting aerodynamic properties of the blades with flaps or slats
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D7/00—Controlling wind motors 
    • F03D7/02—Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/0296—Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor to prevent, counteract or reduce noise emissions
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D7/00—Controlling wind motors 
    • F03D7/02—Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/04—Automatic control; Regulation
    • F03D7/042—Automatic control; Regulation by means of an electrical or electronic controller
    • F03D7/043—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D7/00—Controlling wind motors 
    • F03D7/02—Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/04—Automatic control; Regulation
    • F03D7/042—Automatic control; Regulation by means of an electrical or electronic controller
    • F03D7/043—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic
    • F03D7/045—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic with model-based controls
    • 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
    • F05B2240/00—Components
    • F05B2240/20—Rotors
    • F05B2240/30—Characteristics of rotor blades, i.e. of any element transforming dynamic fluid energy to or from rotational energy and being attached to a rotor
    • F05B2240/31—Characteristics of rotor blades, i.e. of any element transforming dynamic fluid energy to or from rotational energy and being attached to a rotor of changeable form or shape
    • 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
    • F05B2270/00—Control
    • F05B2270/30—Control parameters, e.g. input parameters
    • F05B2270/334—Vibration measurements
    • 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
    • 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
    • Y02P70/00—Climate change mitigation technologies in the production process for final industrial or consumer products
    • Y02P70/50—Manufacturing or production processes characterised by the final manufactured product

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Wind Motors (AREA)

Abstract

描述了一种用于风轮机的控制方法，尤其是用于风轮机叶片的控制方法。该控制方法利用了叶片模态形状或者叶片的自然振动形状来检测叶片自然振动的激励水平，并且控制了在叶片上的主动提升装置，以便为了减小激励水平，减小在叶片和总体风轮机结构中的负载。还提供了一种设计用于在此方法中使用的风轮机叶片的方法。

Description

风轮机叶片控制方法
技术领域
本发明涉及一种用于风轮机叶片的控制方法，以及设计风轮机叶片的方法，以及根据该方法所设计的风轮机叶片。
背景技术
目前研究的主动空气动力装置在风轮机叶片中使用，以便减小结构上随时间变化的负载。主动空气动力装置(AAD)是可以以可变的和受控的方式操纵叶片区段的空气动力的装置的广义用语。主动空气动力装置的实例是后缘襟翼，其很类似于飞机机翼的外侧副翼。通过将此类装置安装在叶片上并且适当控制，由湍流引起的负载可被部分地消除。
不断的挑战在于设计出在一些位置上具有此类AAD以提供最高性能的用于风轮机的叶片，以及控制此类AAD的合适方法。
本发明的一个目的在于提供一种用于使用AAD的风轮机叶片的新控制方法，以及将此类AAD适当地定位在叶片上的新设计方法。
发明内容
因此，提供了一种在风轮机的操作期间控制风轮机叶片以减小叶片的根部动量的方法，该叶片具有末梢端和根部端，该方法包括以下步骤：
识别风轮机叶片的激励模态形状；以及
调整设在所述叶片上的至少一个主动提升装置，以在所述风轮机的操作期间减小所述叶片的模态振动，其中所述调整步骤是基于所述激励模态形状。
通过使风轮机叶片的控制基于叶片模态，提出了优于现有技术的若干优点。基于叶片模态的控制可提供与标准的根部动量控制系统相比更快的控制响应，以及由风轮机所经受的疲劳负载的减小。另外，基于叶片模态的控制提出了在与基于先进模型模拟的控制进行比较时明显减少的计算要求。
优选地，所述识别激励模态形状的步骤包括：
测量或导出所述风轮机叶片的偏转，以及
将所述测得的偏转与所述风轮机叶片的至少一个已知的模态形状相比较，以确定所述至少一个模态形状的激励水平，
其中所述调整步骤包括基于所述激励水平来促动所述至少一个主动提升装置，以减小所述激励模态形状的大小(magnitude)，以减小所述叶片的模态振动。
测量叶片的模态形状的激励水平，并且基于所述激励水平来操作主动提升装置的促动水平。将理解的是，测量或导出叶片的偏转的步骤可基于例如加速计、叶片动量传感器、应变仪、光学偏转传感器、位置传感器(例如GPS传感器)等任何适合的传感器或传感器阵列的输出。
优选地，所述叶片的模态振动基于所述叶片的模态1形状的激励，并且其中该方法包括以下步骤：
将至少第一主动提升装置设在朝向所述叶片末梢的位置处；以及
促动所述至少一个第一主动提升装置来控制所述模态1形状的激励水平，使用所述模态1形状的激励坐标来作为所述至少一个第一主动提升装置的输入。
将理解的是，模态1形状被认作是风轮机叶片的第一自然频率或模态。
将理解的是，激励坐标是风轮机叶片的模态的瞬时激励水平。
激励坐标用作主动提升装置的反馈控制器的输入。这提供了简单的控制系统，其可容易地调节来优化独立叶片设计的性能。
优选地，所述叶片的模态振动基本上由模态1形状的激励和所述叶片的模态2形状的激励的组合提供，并且其中该方法包括以下步骤：
将第一主动提升装置和第二主动提升装置设在叶片上，所述第一主动提升装置和所述第二主动提升装置能有选择地促动来控制所述模态1形状和模态2形状，其中叶片上的所述第一主动提升装置和所述第二主动提升装置的位置被选择成在主动提升装置被促动时提供在模态1形状与模态2形状之间的最小干扰。
使用主动提升装置控制第一叶片模态和第二叶片模态导致相对简单的二阶控制系统，其可容易地适应独立叶片设计的特性。通过最大限度地减小在主动提升装置之间的干涉或串扰，这从叶片模态1和2中的各个的最佳控制提供，而不同时地激励另一个叶片模态。
将理解的是，叶片的模态2形状被认作是风轮机叶片的第二自然频率或模态。
优选地，该方法包括以下步骤：
将第一主动提升装置设在沿着所述模态2形状的结点与叶片的末梢端之间的叶片长度的位置处；以及
将第二主动提升装置设在沿着所述模态2形状的结点与叶片的根部端之间的叶片长度的位置处。
将主动提升装置定位在叶片模态2形状的结点的任一侧上，这提供了叶片模态2形状的最大平衡控制，同时提供了相对主导的叶片模态1形状的有效控制。
将理解的是，叶片模态2形状的结点是指叶片模态形状与法向轴线交叉或相交的点。
优选地，所述比较步骤包括分析所述测得的偏转来确定叶片的模态1形状的激励值和叶片的模态2形状的激励值，并且其中所述促动步骤包括基于叶片的模态1形状和模态2形状的激励值的组合来促动所述第一主动提升装置和所述第二主动提升装置。
通过执行第一模态形状和第二模态形状的分析和基于这些模态形状控制主动提升装置，叶片可被调节来提供叶片和较大的风轮机结构中的疲劳负载和动量的最大控制。
优选地，该方法包括以下步骤：
将叶片的模态1形状的激励值与阈值相比较，以及
当所述模态1激励值超过所述阈值时，执行模态1控制操作，以及
当所述模态1激励值小于所述阈值时，执行模态1控制操作和模态2控制操作。
通过将激励水平与阈值相比较，有可能将控制系统的响应定制成集中于在模态1相当占主导时减小模态1激励。有时，当模态1不太占主导时，即当模态1的激励低于阈值时，则控制系统试图减小模态1和模态2两者。
优选地，执行模态1控制操作的所述步骤包括促动所述主动提升装置来减小激励的叶片模态1形状的大小。
优选地，执行模态2控制操作的所述步骤包括促动所述主动提升装置来减小激励的叶片模态2形状的大小。
优选地，该方法包括以下步骤：
通过在相同方向上促动所述第一主动提升装置和所述第二主动提升装置来控制所述叶片的模态1形状；以及
通过在相反方向上促动所述第一主动提升装置和所述第二主动提升装置来控制所述叶片的模态2形状。
还提供了一种具有控制器的风轮机，该控制器具有储存在计算机可读存储器上的计算机可读指令，这些指令在读取时能操作成实施以上方法步骤中的任何一个。
还提供了一种设计风轮机叶片的方法，其包括：
提供风轮机叶片；
执行所述风轮机叶片的模态分析来确定风轮机叶片的模态1形状和模态2形状；以及
基于至少一个确定的模态形状将至少一个主动提升装置定位在所述风轮机叶片上，使得所述至少一个确定的模态形状的激励水平能由所述主动提升装置的促动来控制。
通过分析叶片模态形状，有可能将主动提升装置定位成用于所述模态形状的最大控制，导致在叶片操作期间减小的叶片疲劳负载。优选地，该设计方法被用于与上述控制方法一起使用来设计叶片。
优选地，所述定位步骤包括在沿着与所述模态2形状的结点对应的叶片长度的位置处将主动提升装置设在所述风轮机叶片上，使得所述主动提升装置能操作成控制所述模态1形状的激励水平，而不会影响所述模态2形状的激励水平。
将主动提升装置设在模态2形状的结点处意味着模态1形状的控制可提供成没有与叶片模态2形状的串扰或干扰，即促动结点处的装置将不影响叶片的模态2形状的激励。通过将主动装置设在该结点处，将理解的是，该装置沿着叶片的长度定位，在模态2形状的结点处定心(centre)。
备选地，所述定位步骤包括将第一主动提升装置和第二主动提升装置设在所述风轮机叶片上，
其中所述第一主动提升装置设在沿着叶片的长度的第一位置处，并且所述第二主动提升装置设在沿着叶片的长度的第二位置处，
其中叶片上的所述第一主动提升装置和所述第二主动提升装置的位置被选择成在主动提升装置被促动时提供模态1形状与模态2形状之间的最小串扰(cross-talk)/干扰(interference)。
优选地，所述第一主动提升装置设在沿着在所述模态2形状的结点与叶片的末梢端之间的叶片长度的位置处，并且其中第二主动提升装置设在沿着在所述模态2形状的结点与叶片的根部端之间的叶片长度的位置处，使得主动提升装置能操作成控制所述模态1形状和所述模态2形状的激励水平。
优选地，所述至少一个主动提升装置包括设在风轮机叶片的后缘处的空气动力襟翼。
所使用的标准命名：
D 阻尼矩阵
K 刚度矩阵
M 质量矩阵
c 翼弦长度[m]
C1翼型件的分段升力系数，Cl=l/(1/2ρU2c)，其中ρ为流体密度。
Faero外部空气动力
Mz围绕平面外轴线的弯矩(襟翼方向动量)
U 流速
x 转子径向位置
y 转子平面外方向的位移
a 翼型件的冲角
ΔC1 升力系数的变化，C1在AAD被触动时
AAD 主动空气动力装置，例如，也可称为主动提升装置的襟翼、阻流板、凸片、流体注入装置
BMS 弯矩传感器，例如，校准的应变仪
Xc 变量X的中点
X∞ 变量X的自由流性质
Xel 变量X的单束元件
附图说明
现在仅通过举例的方式参照附图来描述本发明的实施例，在附图中：
图1示出了风轮机；
图2示出了根据本发明的风轮机叶片的示意图；
图3示出了图2的叶片的翼型轮廓的示意图；
图4为具有主动空气动力装置或主动提升装置的风轮机叶片的示范性视图；
图5为使用的坐标系的图示；
图6为风轮机叶片中的主导模态形状的图表；
图7为襟翼偏转的阶跃响应的一系列图表，示出了弯矩和模态形状的坐标；
图8为AAD的不同叶片模态的可控性的图表；
图9为在叶片的模态1和模态2之间的可控性的一系列图表；
图10为襟翼促动期间对应于模态1和模态2的可控性格兰姆行列式HSV的图表；
图11示出了风轮机叶片上的两个襟翼的优化问题的图表；
图12为由使用的不同控制系统引起的风轮机的不同构件测得的疲劳负载的图表；以及
图13为根据本发明的控制系统的示意图。
将使用相同的参考标号来表示在不同实施例之间的共同元件。
具体实施方式
图1示出了根据所谓的”丹麦构想”的常规的现代逆风风轮机，其具有塔架4、机舱6和具有基本上水平的转子轴的转子。转子包括转毂8和从转毂8沿着径向延伸的三个叶片10，每一个具有最接近转毂的叶片根部端16和最远离转毂8的叶片末梢端14。转子具有表示为R的半径。
图2示出了根据本发明的实施例的风轮机叶片10的第一实施例的示意图。风轮机叶片10具有常规的风轮机叶片形状，并且包括最接近转毂的根部区30、最远离转毂的成型(profiled)区或翼型区34，以及在根部区30与翼型区34之间的过渡区32。当叶片安装在转鼓上时，叶片10包括面朝叶片10的旋转方向的前缘18，以及面朝前缘18的相对方向的后缘20。
翼型区34(也称为成型区)具有关于产生提升的理想的或者几乎理想的叶片形状，而根部区30由于结构考虑具有基本上圆形或椭圆形的截面，其例如使得将叶片10更容易且更安全地安装在转毂上。根部区30的直径(或翼弦)沿着整个根部区域30通常为恒定的。过渡区32具有从根部区30的圆形或椭圆形40逐渐地变至翼型区34的翼型轮廓50的过渡轮廓42。过渡区32的弦长通常随着离转毂的距离r增大而基本上线性地增大。
过渡区34具有翼型轮廓50，其中翼弦在叶片10的前缘18与后缘20之间延伸。翼弦的宽度随着离转毂的距离r增大而减小。
应当注意的是，叶片的不同区段的翼弦一般位于公共平面中，因为叶片可扭转和/或弯曲(即预先弯曲)，从而向翼弦平面提供了对应的扭转和/或弯曲路线，这是最常见的情况，以便补偿取决于离转毂的半径的叶片的局部速度。
图3示出了以各种参数绘出的风轮机的典型叶片的翼型轮廓50的示意图，其通常用于限定翼型件的几何形状。翼型轮廓50具有压力侧52和吸入侧54，在使用期间，即在转子旋转期间，一般分别面朝向上风(或逆风)侧和下风(或顺风)侧。翼型件50具有翼弦60，翼弦60具有弦长c，其在叶片的前缘56与后缘58之间延伸。翼型件50具有厚度t，其限定为压力侧52与吸入侧54之间的距离。翼型件的厚度t沿着翼弦60变化。与对称轮廓的偏离由弧线62给出，弧线62为穿过翼型轮廓50的中线。中线可通过绘出从前缘56到后缘58的内切圆来找到。中线沿着这些内切圆的中心，并且与翼弦60的偏离或距离称为弧f。对称还可通过使用称为上弧和下弧的参数来限定，其分别限定为离翼弦60和吸入侧54和压力侧52的距离。
翼型轮廓通常特征为以下参数：弦长c、最大弧f、最大弧f的位置df、最大翼型厚度t，其为沿着中间弧线62的内切圆的最大直径，最大厚度t的位置dt，以及鼻部半径(未示出)。这些参数通常限定为与弦长c的比。
图4示出了具有被偏转-10度的后缘襟翼11的风轮机转子叶片10的实例。
1.介绍
开发出了用于叶片的简单模型。该模型为叶片的襟翼方向偏转的有限元模型(基本上从转子平面向外的方向)。未考虑与边缘方向(转子平面中)的偏转或与扭转的联接。已经忽略了转子空气动力(引起的尾迹和末梢效应)，但已经包括由叶片的动态偏转引起的空气动力阻尼。为了简单起见，还忽略了陀螺交感和离心刚化。主动空气动力装置的空气动力可被控制。已经使用了逆风5MW参考涡轮叶片来作为实例，其为61.5米长(细节见Jonkman,J.,Butterfield,S.,Musial W.,Scott,G.,2009,Definition of a5-MW Reference WindTurbine for Offshore System Development Definition of a5-MW Reference WindTurbine for Offshore System Development.Tech.Rep.二月,NREL)。
Andersen对AAD的早期研究(Andersen,P.B.,2010的Risoe博士报告AdvancedLoad Alleviation for Wind Turbines using Adaptive Trailing Edge Flaps:Sensoring and Control,丹麦技术大学)示出了应变仪传感器，其用于控制襟翼，应当置于外侧较远(离根部大约28m)，以便控制襟翼。
提出了使用叶片的模态形状来控制AAD。大多数早期研究使用了叶片的应变的局部测量或局部位移测量。模态形状为叶片的固有自由度，并且提出了在试图阻尼叶片的振动时必须使用它们。
2.模型
梁模型由N联接伯努利欧拉梁元件构成，其为标准立方体梁元件。各个梁元件均具有4个自由度，y方向上的两个位移和围绕平面外轴线θ的2个旋转。
例如，Rao S.S.,2004Mechanical Vibration.Prentice Hall的作为元件质量和刚度矩阵的教科书实例。
参看图5，并且示出了使用的坐标系的实例。叶片为离散化成N个结点。Y是转子平面外位移，且θi为梁的局部斜率。θi为图5中所示的情况中的负值。弯曲黑线指出了偏转状态，而直线红线指出了未偏转状态。
元件质量和刚度矩阵
第i个元件已经具有以下刚度矩阵：
其中EI为元件的刚度，且L为长度。Y和θ为元件的结点(端点)处的位移和斜率。
其中M为每元件长度的质量。
阻尼矩阵通过瑞利阻尼引入，其中调节的系数给予叶片本征频率下大约3％的对数衰减。
Dstruc，el=0.188.5Mel+0.0048Kel
2.1空气动力阻尼
通过使用梁元件的立方偏转形状，可计算元件的中心的偏转。元件的中心处的偏转为
yc，i=0.5(yi+yi+1)+L/8(θ1-θ2)
分段升力系数C1基本上为C1＝2πα，其中α为以弧度为单元的冲角。在该简单线性模态中，平均冲角设置为0，且因此使用叶片的相对速度
以及中心点的平移速率，由元件的运动引起的元件上空气动力为
元件的空气动力阻尼矩阵变为：
2.2外力(AAD)
叶片上的外力仅来自AAD，其可将局部升力系数改变一定量ΔC1。如果AAD仅置于第i个结点处，则可将常数对角矩阵公式化
其应当乘以列向量ΔC1，其含有从主动装置对升力系数的贡献。ΔC1将为状态空间模型中的控制输入u。
2.3.状态空间模型
系统的运动方程为
其中q=[y2，θ2，···，yj，θi，···，yN+l，θN+1]T包含叶片的位移和旋转。该形式的常微分方程系统
可对于此特定情况形成:
其中M,K和D为组合的质量、刚度和阻尼矩阵。
系统组件
如果Kel写作
则组合的刚度矩阵写作
同样对于质量和阻尼矩阵完成。
夹紧梁的边界条件为壁处的偏转和斜率为零。这意味着ql＝0。这可通过从所有矩阵除去与第一结点(y1和θ1)相关联的行和列来实现，将其减小至2Nx2N系统。
2.3.1.输出方程1：根部弯矩
如果方程对于系统输出公式化，则我们利用单个输入u(对于该系统ΔC1)和单个输出y来使状态空间公式化。输出方程具有该形式
y=Cx+Du
相关的输出可为叶片的根部弯矩，其通过最内侧区段处的曲率乘以刚度来得出。通过使用立方形状函数，曲率得出为以下
且动量限定为
Mz(x)=y(x)"EI(x) (1)
如果根部动量为期望输出，则C向量变为
2.3.2.输出方程2：模态形状坐标
代替使用给定区段处的弯矩来控制叶片，还尝试使用叶片的模态形状来控制振动。
模态形状已经通过解决本征值问题来计算
(K-λM)·Φ＝0
其给出具有向量中的对应本征值λ的矩阵的模态形状Φ。Φ的列向量按本征值的大小来分类，且最低频率模态在第一列中。
这可用于改变C矩阵，以便输出y为给定的模态形状坐标。模态形状坐标为模态形状的瞬时激励水平。由于模态形状通常在叶片的末梢处标准化为1，故模态形状坐标还可看作是对末梢偏转的贡献。如果第一模态形状的坐标为期望的控制参数/输出，则C矩阵采取该形式：
Cmode1=[1 0 0 ... 0][Φ-1 0]
或对于模态2：
Cmode2=[0 1 0 ... 0][Φ-1 0]
3.结果-促动器放置
最先的3个模态形状在图6中给出，其示出了3个主导模态形状。所有模态形状已被标准化为末梢处的1m偏转。重要的是注意第二模态形状的正偏转将由于根部处的负曲率(见方程1)而给予根部弯矩负贡献。这将证明是系统控制的问题。
3.1.叶片上的单个AAD
3.1.1.阶跃响应
最近的研究已经使用气动弹性代码来发现促动器和传感器的最佳放置(Andersen,P.B.,Henriksen,L.C,Gaunaa,M.,Bak,C,Buhl,T.,2010.Deformable trailingedge flaps for modern megawatt wind turbine controllers using strain gaugesensors.Wind Energy(2009年12月),193-206)。结论是当6.3m襟翼置于离根部大约50m时，应变仪传感器应当放置在离根部大约28m。该结果通过使用重气动弹性代码(HAWC2)的精细参数研究获得。利用到手的新模型，看起来真实的是，外侧传感器位置归因于根部弯矩的非最小相位行为。
在图7中，示出了置于离根部50m处(黑实线)和离根部58m处(红虚线)的6.3m宽襟翼的在t＝0秒处的襟翼偏转的阶跃响应(利用状态空间模型上的Matlab阶跃命令作出)。绘出了六个传感器，即根部弯矩，在28m和40m处的弯矩，以及前三个襟翼方向模态的大小。三个上部的图表示出了3个位置处弯矩；根部、28m处和40m处。三个下部的图表示出了前三个襟翼方向模态的模态形状坐标。
如所见到的，3个模态形状的激励取决于AAD的所处的位置。如果AAD置于离根部(黑色)50m，则第二模态形状几乎没有任何激励，因为AAD置于第二模态形状的结点中(见图6)。
取决于弯矩传感器所处的位置，看到了不同的响应。如果弯矩传感器(BMS)置于根部中，则响应为清楚的非最小相位，尤其是如果襟翼置于外侧更远。这归因于第二模态形状的激励，其具有根部附近的模态形状中的负曲率(见图6)。
如来自方程1，弯矩由叶片的曲率给出。如果MBS置于离根部28m，则看不到来自第二模态形状的响应，因为第二模态形状的曲率在28m处为零。然而，响应仍略微为非最小相位，因为在弯矩上看到3模态形状。如果BMS置于40m处，则3模态形状的曲率为零，并且看不到来自3模态形状的响应。然而，第二模态形状的曲率在40m处不为零，但至少其是正的。这意味着其指向与第一模态形状形同的方向，并且因此响应不再是最小相位。当襟翼置于离根部50m处时，第二模态形状未受激励，因为这是第二模态形状的结点的位置。
第一模态和第二模态的阶跃响应看起来类似简单二阶系统，其对于调节控制器是理想的。因此，模态在以下区段中用作传感器。
3.1.2.参数研究
以下研究尝试使用状态空间模型来分许给定AAD的最佳位置。该分析基于系统的可控性，这意味着系统(具有AAD的叶片)在有限时间内达到给定状态(偏转)的能力。如果系统是可控的，则在有限时间内达到给定状态的成本可使用格兰姆理论来评估。AAD位置和大小的给定构造的可控性成本越小，控制叶片的振动可越好。然而，可控性成本可在第一本征型上较低，意味着AAD可激励第一本征型，但如果第二模态也由相同AAD移动来激励，则叶片的总振动可受到不利影响。
可控性成本ε可由ε＝(1/Λ)估计，其中Λ为已经由Laub的方法计算的矩阵q的汉克尔奇异值(HSV)
Ωq=RTQqR
其中R是分解的解
Pij=RTR
并且Pij为使用B向量的系统的可控性格兰姆行列式，其改变为具有位置i和宽度j处的襟翼的性质。下标q为人们期望研究其可控性的模态形状，且将为以Qq的参数，其为观察模态形状q的系统的可观察性的格兰姆行列式。可控性的计算是很计算有效的(使用Matlab中的lyap函数)，且参数研究可很快进行，且该方法对于在需要许多目标函数估计的因子优化器中使用很可行。
在下文中，用语可控性将表示Λ，其与可控性成本成反比。可控性成本越低，则可控性越高，因此较高的可控性较好。
图8示出了AAD位置和宽度的各种组合的模态1(黑色)、模态2(红色)和模态3(颜色，虚线)的可控性Λ。”不可行区域”标记了AAD延伸到叶片末梢外侧的区域。
在图8中，模态1、模态2和模态3的可控性示为襟翼位置和宽度的函数。将清楚的是，模态1比模态2高一个数量级，并且两个的数量级高于模态3。模态3预计在叶片动力中的重要性很低。第一模态形状的蛤壳形轮廓示出了襟翼越大，则获得的可控性越高。当襟翼较小时，最佳位置在末梢附近，并且在襟翼大小增大时更向内侧移动(且有效地脱离其不贡献的叶片的末梢)。当襟翼较小且置于末梢附近或者襟翼较大且置于第二模态形状的”谷部”中(图6中接近x＝40m)时，第二模态形状受到激励。
有趣的是看到当襟翼很小时，将襟翼放在模态1最可控的位置可能不太明智，因为模态2也将受控/激励。因此，例如在10m的襟翼宽度处，襟翼将可能不置于其给予第一模态的最大可控性的55m处，而是在50m处，在该处，其仍具有75％的可控性，但不同样多地激励模态2。模态1对模态2的较高比率是被期望的。
3.1.3.最佳化
参数研究并未公开C1分布在襟翼上的另一个分布是否可减小模态1与模态2之间的串扰。因子优化器(来自Matlab)用于增大模态1的可控性和减小模态2的可控性，并且作为最大限度减小襟翼宽度的另一个目的。设计变量为襟翼的端部位置以及襟翼的各个端部处的最大ΔC1。ΔC1在襟翼上线性地变化。图9中示出了优化的结果。结果为帕累托前沿，并且因此示出了给定目标函数值的最佳解。
图9示出了以下：
上左：模态1与模态2的可控性Λ之间的帕累托前沿。
下左：模态1的可控性与襟翼宽度之间的帕累托前沿。
上右：模态2的可控性与襟翼宽度之间的帕累托前沿。
下右：对应于为10的模态1的可控性的解。
在上左图中看到，模态1和模态2的可控性遵循达到模态2的可控性在模态1的为10的值周围的宽度增大的点的线性关系。这对应于上右图中的20m的襟翼宽度。在下左图中，示出了模态1的为10的可控性的ΔC1的解。
3.2.叶片上的两个AAD的放置
5m宽度的两个襟翼被各自置于变化的两个独立位置处的叶片上。
两个襟翼的给定组合的可控性又计算为矩阵q的第一(最大)汉克尔奇异值(HSV)，但现在B已经改变成包括位置i和j处的两个襟翼之后Pij为可控性格兰姆行列式。B可改变为使两个襟翼同相(在相同方向上的襟翼)或异相(在相反方向上的襟翼)工作。下标q为人们期望研究其可控性的给定模态形状，并且将为以Qq的参数，其为观察模态形状q的系统的可观察性的格兰姆行列式。
3.2.1.襟翼位置的参数研究
它们又经历一个输入ΔC1，但襟翼可在相同方向上或在相反方向上移动。对于5m宽度的两个襟翼的情况，结果在图10中示出。
图10示出了在2个襟翼被在相同方向上(下三角形)和相反方向(上三角形)促动时对应于模态1(黑色实线)和模态2(红色虚线)的可控性格兰姆行列式HSV。由于两个襟翼位置可互换，故轮廓图围绕中心蓝色对角线对称，并且两个图表已经被组合。
图中的下三角形示出了在襟翼在相同方向上移动时发生什么，并且由于两个襟翼位置可互换，故轮廓图围绕蓝色中线对称。将看到的是，第1模态形状通常比第二模态形状更可控。用于控制第1模态形状的最佳位置在襟翼离根部50m和55m(彼此紧邻)的情况下在更远侧。用于控制第二模态的最佳位置是通过将两个促动器放置成在大约35m和40m处彼此紧邻。然而，当将促动器置于第一模态被最好促动的位置处时，第二模态也被促动。
图10的上三角形示出了两个襟翼被在相反方向上促动的情况。图片现在变化，并且第二模态现在比整个第一模态更可控。当一个襟翼置于很远且另一个尽可能置于内侧时，在该处其不会给予更大的力且实际上停用时，第一模态最可控。当一个襟翼置于较远外侧且一个置于图6中的第二模态形状的局部最小值处的40m时，第二模态很可控。两个襟翼的最佳位置为两个襟翼被在相同方向上促动时模态一很可控的位置，以及当它们被在相反方向上促动时模态2很可控的位置。因子优化器用于示出此权衡。
3.2.2.两个襟翼的优化
优化问题被建立，其中两个目标为：
1.当襟翼在相同方向上(也共同表示)移动时最大限度地增大模态1的可控性
2.当襟翼被在相反方向上促动时最大限度地增大模态2的可控性
用于优化的约束在于襟翼不必重叠，不必在叶片上向内侧延伸40m，不必进一步向外侧延伸61m，并且每一个襟翼必须在转子半径的3％至20％之间。
应当优化的两个目标函数为
其代表具有最适合的控制器(共同或相反)的给定模态的可控性的剩余。当总襟翼长度约束为10m时的解在图11中给出。
图11示出了具有转子半径的20％的总宽度的两个襟翼的优化。
4.最佳传感器位置
直到现在，模态形状已经通过知道沿着叶片的位移并且使用模态对形状的认识来找出模态形状坐标来确定。另一个方法通过使用弯矩传感器，即校准的应变仪。如果确定最先的三个模态形状，并且可忽略所有高阶模态形状，则其将至少需要三个BMS。模态形状坐标Gxi，其中i为简单地通过求解以下得到的模态形状数目
当动量Mz已经利用在位置x1,x2和x3中的BMS测得时。Φi’’(x)为位置x处的模态形状数目i的曲率，并且EI(x)为x处的刚度。
鉴于段落3.1.1，用于负载控制的最佳传感器位置可变窄。例如，如果还可忽视第三模态形状，则仅模态形状1相关来控制。因此，单个传感器可置于第二模态形状的零曲率点上，xp＝28m，以给出模态形状1的偏转。
通过了解第一较少模态形状坐标，可从模态形状发现任何区段处的瞬时弯矩或位移。
5.FLEX5中的AAD的实施
为了评定是否值得控制涡轮上的第1和第2模态两者，使用了气动弹性工具FLEX5。FLEX5为在以时间域工作且广泛用于风轮机行业的气动弹性模拟工具。其由丹麦技术大学的Stig开发。更多信息见(Oeye,S.1996.FLEX4-Simulation of Wind TurbineDynamics.In:Proceedings of the28th IEA Meeting of Experts State of the Art ofAeroelastic Codes for Wind Turbine Calculations.4月11-12.Denmark.Lyngby,丹麦技术大学，第71–76页)。
在下文中，将模拟逆风5MW参考风轮机。涡轮已经使用IEC Ed.3风级1B来模拟。地基在地面水平处是刚性的。在此研究中已经仅考虑了正常运算，正常湍流模型。对于5m/s到25m/s风速的每0.5m/s，将对于疲劳模拟和分析10分钟时间序列。疲劳负载通过IEC Ed.3的标准程序来计算，即通过计算各个风速的各个传感器的负载的雨流。给予的疲劳负载为20年寿命的给定传感器的破坏等效负载。对于叶片负载，使用了m＝10的系数，但对于塔底部，传动系和动量，使用了m＝3的值。极端负载基于相同量的负载情况，仅利用极端断流模型。
两个AAD在翼展位置x＝47m和x＝54m处安装到叶片上，各个都具有6.3m的宽度。这对应于在70％至80％的转子半径处的襟翼，以及从80％到90％的襟翼。10％c的襟翼用于NACA64618上，其中偏转限制为±10°：襟翼数据使用CFD求解器Ellipsys2D计算。
更多细节见：Baek,P.,Gaunaa,M.,Unsteady Wind Tunnel Results for aMiniflap and a Trailing Edge Flap on a Wind Turbine Airfoil.Wind Energy,1-9；Michelsen,J.A.,1992.Basis3D-a Platform for Development of Multiblock PDESolvers-AFM92-05.Tech.rep.,丹麦技术大学；Michelsen,J.A.,1994.Block structuredMultigrid Solution of2D and3D Elliptic PDEs-AFM94-06.Tech.rep.,丹麦技术大学；Soerensen,N.N.,1995.General Purpose Flow Solver applied to Flow over Hills-Risoe-R-827(EN).Tech.rep.,RisoNational Laboratory,丹麦技术大学。
襟翼可将基底极线在α＝5°处改变ΔC1＝0:6。
对于所有情况，使用了相同的控制器，但是具有不同的调节参数。对于所有情况，SISO Plead控制器调节成具有7rad/s的交叉频率，其中相位容限为60°。P-lead控制器是有效的相同比例导数控制。
对于单独的情况计算：
基本负载设置。没有主动控制。仅普通共同节距变速控制器。
在x＝1:5m处。SISO控制向两个襟翼给出相同输入。对SISO控制器的输入为根部弯矩。
。SISO控制器向两个襟翼给出相同的输入。SISO控制器的输入为28m处的弯矩动量。这几乎是Andersen的控制策略(Andersen P.B.r.,2010.Advanced Load Alleviationfor Wind Turbines using Adaptive Trailing Edge Flaps:Sensoring and Control-博士论文)。
模态1控制器。模态1形状坐标输入到促动两个襟翼的SISO控制。两个AAD作为一个移动，其中反馈控制器监测离根部x＝0m处的叶片襟翼方向弯矩。
模态1+2控制器。模态1形状坐标输入至在相同方向上促动两个襟翼的SISO控制。模态2形状坐标输入至促动两个襟翼但在相反方向上的第二SISO控制器。模态2形状控制器具有15rad/s的交叉频率，其中相角为60°。
不同控制器在图12中比较，其示出了不同控制器关于基本情况的各种涡轮构件的疲劳负载(上)和极端负载(下)。对于疲劳负载，看到四个控制器给予襟翼方向负载的类似减小，例外是MBSx＝28m控制器。对于其它风轮机构件，模态形状控制器具有优点，尤其是具有较高带宽并且可减小破坏钢构件的低循环疲劳的模态形状1+2控制器。这清楚地示出了提出的新控制方法用于在涡轮操作寿命期间减小疲劳和极端负载。
6.结论
使用具有主动空气动力装置的逆风5MW参考涡轮叶片的有限元梁模型，探究了促动器位置和传感器位置对系统行为的影响。AAD被认为是阻尼和减小例如通过风轮机叶片上的湍流引起的波动负载。为了研究该系统，调查了多少给定的AAD位置可激励叶片，合理的是，如果系统可被激励，则其还可被主动地阻尼。
关于促动器位置，发现了单个AAD置于越外侧，则第一模态形状的可控性就更好。然而，将AAD置于更外侧将激励高阶模态形状，其中第二模态形状最重要。对第二模态形状的影响通过将AAD置于第二模态形状的结点中减小，但代价是减小了第1模态形状的可控性。示出了当AAD置于更外侧时，如果系统的输出是根部弯矩，则第二模态形状的影响将产生系统最小相位。AAD的目的在于控制根部弯矩，其将负载传递至涡轮的其余部分。
模型验证
有限元模型可相对于教科书实例[9]验证。对于具有恒定刚度、长度和单位质量的梁，则自然频率在表1中给出，并且看到了具有20元件的有限元梁的优异一致。
模态 ω，理论[rad/s] ω，具有21结点的模型(rad/s)
1 3.516 3.516
2 22.034 22.035
3 61.7010 61.698
4 120.912 120.910
表1：具有刚度、长度和单位质量的梁。
通过适当分析和监测叶片模态形状，有可能高效地且准确地控制操作期间的叶片振动和负载。
参看图13，示出了风轮机叶片10的控制方法的示意图。叶片10联接到控制器(未示出)上，控制器包括处理装置和计算机可读储存介质。控制器联接到沿着叶片10的长度设置的多个传感器100上，传感器能操作成检测叶片的偏转。传感器可包括任何适合阵列的偏转检测元件，例如加速计、叶片动量传感器、应变仪、光学偏转传感器、位置传感器(例如，GPS)传感器等。
在风轮机叶片10在风轮机上操作时，叶片10将经历基于旋转期间遇到的力的偏转。传感器100能操作成输出该偏转102的测量结果。
控制器接收偏转测量结果102，并且执行偏转形状的数学分析，以从偏转102确定叶片104的模态1形状的激励水平和叶片106的模态2形状的激励水平。
基于不同模态激励水平或值102,104，控制器开始不同的控制方案。模态1形状102的激励水平用于调节模态1控制系统108，而模态2形状106的激励水平用于调节模态2控制系统110。
控制系统108,110的输出用作在叶片后缘处设在叶片10上的一对主动提升装置112,114的输入。第一主动提升装置112朝向叶片10的末梢端部定位，而第二主动提升装置114定位在第一主动装置112与叶片的根部端之间，优选地在叶片10的模态2形状的结点的位置的相对侧上。主动提升装置112,114可包括能够调整主动装置的位置处的升力特性的任何装置，例如襟翼、阻流板、凸片、流体注入装置等。
第一主动提升装置112基于模态1控制器108和模态2控制器110的输出的总和116来受控制。第二主动提升装置114基于模态1控制器108和模态2控制器110的输出的负加成118来受控制。主动提升装置112,114因此在反馈系统中受控以减小叶片的振动的大小，并且因此减小叶片和风轮机中的负载。
为了控制叶片10的模态1形状，主动提升装置112,114在相同方向上被促动。为了控制叶片10的模态2形状，主动提升装置112,114在相反方向上被促动。
基于叶片的偏转且具体是叶片模态的激励水平来改变主动提升装置112,114提供了使有效性与减小的计算复杂性组合的控制叶片振动的方法。以此方式在模态2形状的结点的任一侧上控制主动提升装置112,114确保了叶片模态之间的最小串扰或干扰，这用于减小叶片振动的大小。
在本发明的另一个增强方案中，应用的控制方案可基于叶片中检测到振动的大小来调节。对于相对较低的振动，可能足以控制模态1形状，同时如果振动相对较高，则其将优选为执行模态1形状和模态2形状的控制。
例如，当叶片偏转102低于阈值时，其可足以促动主动提升装置112,114来仅控制叶片的模态1形状。在图13中所示的系统中，这可通过在叶片偏转102低于预先限定的阈值偏转时在模态2控制器110停用一定周期，并且仅操作模态1控制器108用于低于阈值的偏转来实现。如果叶片偏转102超过阈值，则控制系统能操作成启动模态2控制器110，从而操作控制方案来减小叶片10的模态1形状和模态2形状两者。
尽管以上实例示出了叶片的模态1激励和模态2激励(这些模态形状为叶片激励和振动期间的主导模态形状)，但是将理解的是，可将类似系统用于仅控制主导模态1形状。此系统可能需要控制仅单个的主动提升装置，其可能定位在模态2形状的结点处，以确保模态1形状的主动提升装置的任何促动不会单独地激励模态2形状。
将理解的是，控制系统可对于结合叶片、传感器和/或主动提升装置的任何适合构造变化。在一些实施例中，使用叶片模态形状的坐标值作为对主动提升装置的输入来控制该特定的叶片模态可能就足够了。
本发明不限于本文所述的实施例，并且可在不脱离本发明的范围的情况下改变或改动。

Claims (17)

1.一种在风轮机的操作期间控制风轮机叶片以减小所述叶片的根部动量的方法，所述叶片具有末梢端和根部端，所述方法包括以下步骤：
识别所述风轮机叶片的激励模态形状；调整设在所述叶片上的至少一个主动提升装置，以在所述风轮机的操作期间减小所述叶片的模态振动，其中所述调整步骤是基于所述激励模态形状，并且其中所述叶片的模态振动至少基于所述叶片的模态1形状的激励，并且
其中所述方法包括以下步骤：
将至少第一主动提升装置设在朝向所述叶片末梢的位置处；以及
促动所述至少第一主动提升装置来控制所述模态1形状的激励水平，使用所述模态1形状的激励坐标作为所述至少一个第一主动提升装置的输入。
2.根据权利要求1所述的方法，其中所述识别激励模态形状的步骤包括：
测量所述风轮机叶片的偏转，以及
将所述测得的偏转与所述风轮机叶片的至少一个已知的模态形状相比较，以确定所述至少一个模态形状的激励水平，
其中所述调整步骤包括基于所述激励水平促动所述至少一个主动提升装置，以减小所述激励模态形状的大小，以减小所述叶片的模态振动。
3.根据权利要求1所述的方法，其中所述叶片的模态振动至少基于所述叶片的模态1形状的激励，并且其中所述方法包括以下步骤：
将至少第一主动提升装置设在朝向所述叶片末梢的位置处；以及
促动所述至少一个第一主动提升装置来控制所述模态1形状的激励水平，使用所述模态1形状的激励坐标作为所述至少一个第一主动提升装置的输入。
4.根据权利要求1-3中任一项所述的方法，其中所述叶片的模态振动基本上通过所述叶片的模态1形状的激励和模态2形状的激励的组合来提供，并且其中所述方法包括以下步骤：
将第一主动提升装置和第二主动提升装置设在所述叶片上，所述第一主动提升装置和第二主动提升装置能有选择地促动来控制所述模态1形状和模态2形状，其中在所述叶片上的所述第一主动提升装置和第二主动提升装置的位置被选择成在这些主动提升装置被促动时提供在所述模态1形状与模态2形状之间的最小干扰。
5.根据权利要求4所述的方法，其中所述方法包括以下步骤：
将第一主动提升装置设在沿着在所述模态2形状的结点与所述叶片的末梢端之间的所述叶片的长度的位置处；以及
将第二主动提升装置设在沿着在所述模态2形状的结点与所述叶片的根部端之间的所述叶片的长度的位置处。
6.根据权利要求4所述的方法，其中所述比较步骤包括分析所述测得的偏转来确定所述叶片的模态1形状的激励值和所述叶片的模态2形状的激励值，并且其中所述促动步骤包括基于所述叶片的模态1形状的所述激励值和所述叶片的模态2形状的所述激励值的组合来促动所述第一主动提升装置和所述第二主动提升装置。
7.根据权利要求5所述的方法，其中所述比较步骤包括分析所述测得的偏转来确定所述叶片的模态1形状的激励值和所述叶片的模态2形状的激励值，并且其中所述促动步骤包括基于所述叶片的模态1形状的所述激励值和所述叶片的模态2形状的所述激励值的组合来促动所述第一主动提升装置和所述第二主动提升装置。
8.根据权利要求4所述的方法，其中所述方法包括以下步骤：
将所述叶片的模态1形状的激励值与阈值相比较，以及
当所述模态1激励值超过所述阈值时，执行模态1控制操作，以及
当所述模态1激励值小于所述阈值时，执行模态1控制操作和模态2控制操作。
9.根据权利要求8所述的方法，其中所述执行模态1控制操作的步骤包括促动所述这些主动提升装置以减小所述激励叶片模态1形状的大小。
10.根据权利要求8所述的方法，其中所述执行模态2控制操作的步骤包括促动所述这些主动提升装置以减小所述激励叶片模态2形状的大小。
11.根据权利要求4所述的方法，其中所述方法包括以下步骤：
通过在相同方向上促动所述第一主动提升装置和所述第二主动提升装置来控制所述叶片的模态1形状；以及
通过在相反方向上促动所述第一主动提升装置和所述第二主动提升装置来控制所述叶片的模态2形状。
12.一种风轮机，其包括至少一个风轮机叶片和具有计算机可读存储器的控制器，所述计算机可读存储器具有储存在其上的指令，计算机可读指令在读取时能操作成实施权利要求1-11中任一项所述的方法。
13.一种设计风轮机叶片的方法，其包括：
提供风轮机叶片；
执行所述风轮机叶片的模态分析来确定所述风轮机叶片的模态1形状和模态2形状；以及
基于所述模态1形状和/或所述模态2形状将至少一个主动提升装置定位在所述风轮机叶片上，使得所述模态1形状和/或所述模态2形状的激励水平能通过所述主动提升装置的促动来进行控制。
14.根据权利要求13所述的方法，其中所述定位步骤包括在沿着与所述模态2形状的结点相对应的叶片长度的位置处将主动提升装置设在所述风轮机叶片上，使得所述主动提升装置能操作成控制所述模态1形状的激励水平，而不会影响所述模态2形状的激励水平。
15.根据权利要求13所述的方法，其中所述定位步骤包括将第一主动提升装置和第二主动提升装置设在所述风轮机叶片上，
其中所述第一主动提升装置设在沿着所述叶片的长度的第一位置处，并且所述第二主动提升装置设在沿着所述叶片的长度的第二位置处，
其中在所述叶片上的所述第一主动提升装置和所述第二主动提升装置的位置被选择成在这些主动提升装置被促动时提供在所述模态1形状与所述模态2形状之间的最小串扰。
16.根据权利要求15所述的方法，其中所述第一主动提升装置设在沿着所述模态2形状的结点与所述叶片的末梢端之间的叶片长度的位置处，并且其中第二主动提升装置设在沿着所述模态2形状的结点与所述叶片的根部端之间的叶片长度的位置处，使得这些主动提升装置能操作成控制所述模态1形状的激励水平和所述模态2形状的激励水平。
17.一种根据权利要求13-16中任一项所述的方法所设计的风轮机叶片。
CN201280061666.7A 2011-12-15 2012-12-05 风轮机叶片控制方法 Active CN104093971B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
GB1121590.2 2011-12-15
GB201121590A GB201121590D0 (en) 2011-12-15 2011-12-15 A wind turbine blade control method
PCT/EP2012/074442 WO2013087468A1 (en) 2011-12-15 2012-12-05 A wind turbine blade control method

Publications (2)

Publication Number Publication Date
CN104093971A CN104093971A (zh) 2014-10-08
CN104093971B true CN104093971B (zh) 2017-10-13

Family

ID=45560521

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280061666.7A Active CN104093971B (zh) 2011-12-15 2012-12-05 风轮机叶片控制方法

Country Status (6)

Country Link
US (1) US9909562B2 (zh)
EP (1) EP2791502B1 (zh)
CN (1) CN104093971B (zh)
DK (1) DK2791502T3 (zh)
GB (1) GB201121590D0 (zh)
WO (1) WO2013087468A1 (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10385826B2 (en) * 2014-09-12 2019-08-20 Ge Infrastructure Technology, Llc Wind turbine air deflector system control
US9926911B2 (en) * 2014-09-12 2018-03-27 Ge Infrastructure Technology, Llc Wind turbine air deflector system control
CN107784143B (zh) * 2016-08-30 2018-10-30 北京金风科创风电设备有限公司 风力发电机组叶片气弹稳定性的确定方法及装置
DK179472B1 (en) * 2017-05-22 2018-11-27 Envision Energy (Denmark) Aps Blade with pre-deflection for downwind type wind turbine
US10539119B2 (en) * 2017-07-10 2020-01-21 WindESCo, Inc. System and method for augmenting control of a wind turbine assembly
CN109519341A (zh) * 2018-10-16 2019-03-26 安徽科技学院 一种风力机叶片及其风力机
EP3667065A1 (en) * 2018-12-13 2020-06-17 Siemens Gamesa Renewable Energy A/S Removing vibrations in wind turbine blades
EP3667074A1 (en) * 2018-12-13 2020-06-17 Siemens Gamesa Renewable Energy A/S Device and method of damping front and backward movements of a tower of a wind turbine
EP3667066A1 (en) * 2018-12-13 2020-06-17 Siemens Gamesa Renewable Energy A/S Reduction of oscillating wind turbine rotor blade movements
EP3667064A1 (en) * 2018-12-13 2020-06-17 Siemens Gamesa Renewable Energy A/S Damping vibrations in a wind turbine
CN109946029B (zh) * 2019-02-22 2020-07-10 北京理工大学 一种涡轮增压器的多频疲劳试验方法
CN110069822B (zh) * 2019-03-22 2021-02-26 西安交通大学 一种叶片动应变测量的传感器布置方法
US11409249B1 (en) * 2020-01-30 2022-08-09 The Mathworks, Inc. Simulating transverse motion response of a flexible rotor based on a parameter dependent eigenmodes
US11462204B2 (en) * 2020-09-08 2022-10-04 Siemens Gamesa Renewable Energy A/S Wind turbine and method for noise reduction for a wind turbine
EP4050206A1 (en) * 2021-02-25 2022-08-31 Siemens Gamesa Renewable Energy A/S Controlling a wind turbine with respect to dynamic stability
CN113051679B (zh) * 2021-03-18 2021-10-26 湖南南方宇航高精传动有限公司 一种风电齿轮箱主轴承的载荷处理方法
EP4180655A1 (en) * 2021-11-15 2023-05-17 Siemens Gamesa Renewable Energy A/S Actively damping vibrations in wind turbines

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011000628A2 (en) * 2009-06-30 2011-01-06 Vestas Wind Systems A/S Improved wind turbine blade control
WO2011057633A2 (en) * 2009-11-11 2011-05-19 Vestas Wind Systems A/S Improved control of wind turbine blade lift regulating means
EP2350452A2 (en) * 2008-10-14 2011-08-03 Vestas Wind Systems A/S Wind turbine blade with device for changing the aerodynamic surface or shape
WO2011101079A1 (de) * 2010-02-22 2011-08-25 Repower Systems Ag Verfahren zum betrieb einer windenergieanlage

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2521045C (en) 2003-03-31 2010-03-23 Dan Christian Bak Control of power, loads and/or stability of a horizontal axis wind turbine by use of variable blade geometry control
GB2469854A (en) * 2009-04-30 2010-11-03 Vestas Wind Sys As Wind turbine rotor blade
US8321062B2 (en) * 2009-11-05 2012-11-27 General Electric Company Systems and method for operating a wind turbine having active flow control
DE102010013429A1 (de) 2010-01-25 2011-07-28 Smolka, Peter P., Dr., 48161 Leistungsstarkes Rotorblatt

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2350452A2 (en) * 2008-10-14 2011-08-03 Vestas Wind Systems A/S Wind turbine blade with device for changing the aerodynamic surface or shape
WO2011000628A2 (en) * 2009-06-30 2011-01-06 Vestas Wind Systems A/S Improved wind turbine blade control
WO2011057633A2 (en) * 2009-11-11 2011-05-19 Vestas Wind Systems A/S Improved control of wind turbine blade lift regulating means
WO2011101079A1 (de) * 2010-02-22 2011-08-25 Repower Systems Ag Verfahren zum betrieb einer windenergieanlage

Also Published As

Publication number Publication date
GB201121590D0 (en) 2012-01-25
WO2013087468A1 (en) 2013-06-20
US9909562B2 (en) 2018-03-06
EP2791502A1 (en) 2014-10-22
CN104093971A (zh) 2014-10-08
EP2791502B1 (en) 2016-06-15
DK2791502T3 (en) 2016-09-26
US20140334928A1 (en) 2014-11-13

Similar Documents

Publication Publication Date Title
CN104093971B (zh) 风轮机叶片控制方法
CA2758840C (en) Pressure based load measurement
Barlas et al. State of the art and prospectives of smart rotor control for wind turbines
RU2485013C2 (ru) Способ активной деформации аэродинамического профиля
Pechlivanoglou Passive and active flow control solutions for wind turbine blades
Castaignet et al. Full‐scale test of trailing edge flaps on a Vestas V27 wind turbine: active load reduction and system identification
CN103890383B (zh) 控制风力涡轮机的方法及相关系统
Rasmussen et al. Present status of aeroelasticity of wind turbines
Andersen Advanced load alleviation for wind turbines using adaptive trailing edge flaps: sensoring and control
CA2724311A1 (en) Monitoring of wind turbines
EP2275671A1 (en) System and method for designing airfoils
Andersen et al. Load alleviation on wind turbine blades using variable airfoil geometry
KR20140113515A (ko) 블레이드의 균형유지 방법
DK2715122T3 (en) A WAY TO OPERATE A WINDMILL
Berg et al. Design, fabrication, assembly and initial testing of a SMART rotor
Madsen et al. The potentials of the controllable rubber trailing edge flap (CRTEF)
Marten et al. Numerical and experimental investigation of trailing edge flap performance on a model wind turbine
CN112283051B (zh) 一种基于升力线模型的振动信号特征优化方法及系统
Niether et al. Development of a fluidic actuator for adaptive flow control on a thick wind turbine airfoil
Barlas et al. Aerodynamic and load control performance testing of a morphing trailing edge flap system on an outdoor rotating test rig
Abdelrahman et al. Development of a wind turbine test rig and rotor for trailing edge flap investigation: Static flap angles case
Luhmann et al. Coupled aero-elastic multi-body simulation of two-bladed wind turbines in wind farm arrays
AU2010236656B2 (en) Pressure based load measurement
Avalos et al. Aerodynamic techniques to mitigate the 3D loss in the power coefficient of vertical axis wind turbines
CN116718343A (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