CN105888971B - 一种大型风力机叶片主动降载控制***及方法 - Google Patents
一种大型风力机叶片主动降载控制***及方法 Download PDFInfo
- Publication number
- CN105888971B CN105888971B CN201610274672.0A CN201610274672A CN105888971B CN 105888971 B CN105888971 B CN 105888971B CN 201610274672 A CN201610274672 A CN 201610274672A CN 105888971 B CN105888971 B CN 105888971B
- Authority
- CN
- China
- Prior art keywords
- signal
- control
- wing flap
- wind
- blade
- 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
- 238000000034 method Methods 0.000 title claims abstract description 22
- 239000000835 fiber Substances 0.000 claims abstract description 30
- 238000006243 chemical reaction Methods 0.000 claims abstract description 11
- 238000005457 optimization Methods 0.000 claims abstract description 9
- 238000012544 monitoring process Methods 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 9
- 230000003044 adaptive effect Effects 0.000 claims description 8
- 238000002474 experimental method Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 230000009897 systematic effect Effects 0.000 claims description 3
- 230000001276 controlling effect Effects 0.000 claims 1
- 239000013307 optical fiber Substances 0.000 claims 1
- 230000001869 rapid Effects 0.000 claims 1
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 abstract 1
- 230000001419 dependent effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
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/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/044—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic with PID control
-
- 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/32—Wind speeds
-
- 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/70—Type of control algorithm
- F05B2270/706—Type of control algorithm proportional-integral-differential
-
- 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/80—Devices generating input signals, e.g. transducers, sensors, cameras or strain gauges
- F05B2270/808—Strain gauges; Load cells
-
- 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
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
本发明涉及一种大型风力机叶片主动降载控制***及方法。光纤应变传感器测量叶片根部的应力值,传输至控制单元,控制单元包括使用教与学算法进行参数寻优的PID控制器和对叶片根部的应变力信号进行处理的模糊控制器,通过模糊控制器与PID控制器对襟翼摆角进行切换控制,综合了模糊控制与PID控制各自的优点;同时,机前风场信号传感设备测量风力机前风速并传输至控制单元,控制单元中的前馈控制器通过实时监测机前风速变化,计算出因降低随机风或湍流风引起的不均匀载荷所需的控制量;控制单元将两部分控制量进行耦合,完成对襟翼摆角的控制。本发明有效降低了叶片根部的应变力,延长了叶片的使用寿命,降低了风力机组的使用成本。
Description
技术领域
本发明属于风力发电技术领域,特别涉及一种大型风力机叶片主动降载控制***及方法。
背景技术
随着风力发电的高速发展,风力机也逐渐向着离岸化、大型化的趋势发展,这就对风力机的核心部件—风机叶片提出了更高的要求。风机大型化意味着增加了风机的载荷和***的质量,使疲劳载荷和极限载荷急剧增加,严重降低了机组的使用寿命,同时也增加了硬件和工程成本以及后期维护费用。
为了降低风机的疲劳载荷和极限载荷,提高风机寿命、发电量和出力的稳定性,近年来产生了一种含有主动尾缘襟翼的大型风力机叶片,其硬件结构和控制方法与传统的叶片不尽相同。由于该叶片结构的复杂性,建立其精确的数学模型是非常困难的,并且随机风和湍流风对风机叶片的载荷影响不可忽视,因此常规的单级PID控制很难满足叶片的降载需求。
发明内容
针对现有技术不足,本发明提供了一种大型风力机叶片主动降载控制***及方法。
一种大型风力机叶片主动降载控制***,所述***包括信号采集模块、控制模块和电气伺服模块;
所述信号采集模块包括光纤应变传感器、机前风场信号传感设备和光纤应变信号处理设备;所述控制模块包括5路低通滤波器,5路模数转换器,4路控制单元,PLC,4路数模转换器和4路信号隔离器;所述电气伺服模块包括4路襟翼作动器驱动电路和4路襟翼作动器;
所述光纤应变传感器安装在风力机叶片根部,并与光纤应变信号处理设备连接;光纤应变信号处理设备的1路信号输出和机前风场信号传感设备对应4路襟翼的4路信号输出分别对应1路低通滤波器和1路模数转换器顺次连接;与光纤应变信号处理设备对应的1路模数转换器分别连接至4路控制单元,与机前风场信号传感设备对应的4路模数转换器一一对应连接至4路控制单元,每路控制单元通过PLC分别对应1路数模转换器、1路信号隔离器、1路襟翼作动器驱动电路和1路襟翼作动器顺次连接;
所述光纤应变传感器用于采集叶片根部的应变力信号,所述光纤应变信号处理设备用于将光纤应变传感器采集的信号转换为电压信号;所述机前风场信号传感设备用于测量风力机前风速;光纤应变信号处理设备和机前风场信号传感设备分别将信号传送至对应的低通滤波器,所述低通滤波器用于滤掉高频干扰信号;所述模数转换器用于将模拟信号转换成数字信号;所述控制单元包括对机前风场信号进行处理的前馈控制器,使用教与学算法进行参数寻优的PID控制器和对叶片根部的应变力信号进行处理的模糊控制器,用于减小叶片根部应力的控制运算;所述PLC用于切换PID控制器和模糊控制器的输出信号,并耦合来自前馈控制器的控制信号;所述数模转换器用于将数字信号转换成模拟信号;所述信号隔离器用于将控制***输出信号和电气伺服模块隔离开;所述襟翼作动器驱动电路产生驱动襟翼作动器的电信号;所述襟翼作动器根据襟翼作动器驱动电路输出信号调节襟翼使其产生不同的摆角。
所述前馈控制器是在风力发电机组运行时,通过实时监测机前风速变化,计算出因降低随机风或湍流风引起的不均匀载荷所需的控制量。
上述一种大型风力机叶片主动降载控制***的控制方法,具体包括以下步骤:
步骤1:对模糊控制器FCi、前馈控制器FBi和PID控制器PIDi进行初始化,i=1,2,3,4;
步骤2:读取当前风速和叶片根部的应力值,将得到的叶片根部应力值y(k)与叶片根部额定应力值r(k)进行差值运算,得到应力偏差e(k)及偏差变化率ec(k),其中叶片额定应力值r(k)是由风机出厂前实验测得;
步骤21:将步骤2得到的应力偏差e(k)、偏差变化率ec(k)作为模糊控制器FCi的输入变量;
步骤22:选择隶属度函数进行模糊化,并依据模糊规则得到襟翼的控制量,反模糊化后求得模糊控制器FCi的输出变量,此输出为襟翼控制期望角θ1i;
步骤23:将步骤21的应力偏差e(k)作为PID控制器PIDi的输入信号,利用教与学寻优算法对PIDi参数KPi,KIi,KDi进行在线自整定,所述PID控制器PIDi的输出变量为襟翼控制期望角θ2i;
步骤24:采集风力发电机组的机前风速v(t),将其作为自变量,将襟翼摆动角θi作为因变量,对襟翼角度-风速进行拟合,建立襟翼角度-风速的前馈控制器FBi模型:θi(v)=a0+a1v+a2v2+L+anvn,采用最小二乘法确定各项系数;将风速作为前馈控制器FBi的输入信号,则前馈控制器FBi的输出变量为襟翼控制期望角θ3i;
步骤3:将步骤22与步骤23得到的襟翼期望角控制信号分别送往PLC进行处理,在PLC中设置切换算法:当叶根应力误差大于设定值时使用来自模糊控制器FCi的信号θ1i,来自PID控制器PIDi的控制信号θ2i将不起作用;当叶根应力误差小于设定值时使用来自PID控制器PIDi的控制信号θ2i,来自模糊控制器FCi的信号θ1i将不起作用;最后,PLC再耦合来自前馈控制器FBi的控制信号θ3i,并将这些信号传输给对应的襟翼作动器驱动电路;
步骤4:4个襟翼作动器分别接受来自4个襟翼作动器驱动电路的信号,执行襟翼摆动动作以减小叶片根部应力;
上述步骤2-4反复运行,使叶片根部应力值y(k)等于或接近风力机出厂前实验测得的叶片根部额定应力值r(k),直至完成控制任务。
用于PID参数的教与学寻优算法包括以下步骤:
步骤1):设置初始参数,搜索区域范围定义为X=(x1,x2,…,xd)∈[L,U],L=(L1,L2,…,Ld)是空间下界,U=(U1,U2,…,Ud)是空间上界,d为优化问题的维数,d维空间第i个学员定义为学员规模为N,最大迭代次数为maxgen;
步骤2):教师的教学阶段:
步骤21):计算每个学生的适应值,选择最好个体作为老师Xteacher,计算个体平均值然后根据学员与个体平均水平的差异进行学习,如下式:
TFi=2-gen/maxgen(2)
式中:和分别表示第i个学员学习前和学习后的值;w1=1-gen/maxgen为自适应权系数;ri为0-1之间的随机数;TFi为1-2之间的某个数,其值随迭代次数的变化而变化;ri和TFi用于调整学习速率;gen与maxgen分别为当前迭代次数与最大迭代次数;
步骤22):学员更新:
如果的适应值比的适应值好,那么用代替否则,继续使用
步骤3):学员之间相互学习阶段:
步骤31):每个学员Xi在班级中随机选取一个学习对象Xj(j≠i),Xi通过分析自己和学员Xj之间的差异进行学习调整,如下式:
若Xi优于Xj,
若Xj优于Xi,
式中:w2=1-gen/maxgen为自适应权系数;ri为0-1之间的随机数;
步骤32):学员更新:
如果的适应值比的适应值好,那么用代替否则,继续使用
步骤4):根据适应度函数计算每个学员的适应度值,其公式如下:
式中,e(t)为***误差,根据适应度函数更新学员的全局最优值,当计算所得到最优值达到设定值或算法达到最大迭代次数时,退出教与学寻优算法,否则返回步骤2)。
本发明的有益效果为:
1.相对于传统的单级PID控制***,本***中通过引入前馈控制器可有效降低由于风的随机性波动对叶片根部应力产生的影响。2.使用模糊控制与PID控制在设定值处进行切换控制,既利用了模糊控制器不需要建立精确地数学模型、能适应被控对象非线性和时变性的优点,又利用PID控制器算法简单、稳定性好的特点。3.使用教与学算法寻优速度快,求解精度高的特点,在线优化PID控制器的三个参数,从而有效降低了叶片根部的应变力,延长了叶片的使用寿命,降低了风力机组的使用成本。4.本发明设计简单,应用方便,控制更加准确可靠,非常适合大型风力机叶片的建模和控制。
附图说明
图1为一种大型风力机叶片主动降载控制***的结构框图。
图2是一种大型风力机叶片主动降载控制***的控制原理框图。
图3为对PID控制器参数进行在线自整定的教与学寻优算法流程图。
图4为本发明一种大型风力机叶片主动降载控制方法与现有方法的对比结果。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
如图1所示一种大型风力机叶片主动降载控制***,所述***包括信号采集模块、控制模块和电气伺服模块;
所述信号采集模块包括光纤应变传感器、机前风场信号传感设备和光纤应变信号处理设备;所述控制模块包括5路低通滤波器,5路模数转换器,4路控制单元,PLC,4路数模转换器和4路信号隔离器;所述电气伺服模块包括4路襟翼作动器驱动电路和4路襟翼作动器;
所述光纤应变传感器安装在风力机叶片根部,并与光纤应变信号处理设备连接;光纤应变信号处理设备的1路信号输出和机前风场信号传感设备对应4路襟翼的4路信号输出分别对应1路低通滤波器和1路模数转换器顺次连接;与光纤应变信号处理设备对应的1路模数转换器分别连接至4路控制单元,与机前风场信号传感设备对应的4路模数转换器一一对应连接至4路控制单元,每路控制单元通过PLC分别对应1路数模转换器、1路信号隔离器、1路襟翼作动器驱动电路和1路襟翼作动器顺次连接;
所述光纤应变传感器用于采集叶片根部的应变力信号,所述光纤应变信号处理设备用于将光纤应变传感器采集的信号转换为电压信号;所述机前风场信号传感设备用于测量风力机前风速;光纤应变信号处理设备和机前风场信号传感设备分别将信号传送至对应的低通滤波器,所述低通滤波器用于滤掉高频干扰信号;所述模数转换器用于将模拟信号转换成数字信号;所述控制单元用于减小叶片根部应力的控制运算,其包括对机前风场信号进行处理的前馈控制器,所述前馈控制器是在风力发电机组运行时,通过实时监测机前风速变化,计算出因降低随机风或湍流风引起的不均匀载荷所需的控制量;使用教与学算法进行参数寻优的PID控制器;以及对叶片根部的应变力信号进行处理的模糊控制器;所述PLC用于切换PID控制器和模糊控制器的输出信号,并耦合来自前馈控制器的控制信号;所述数模转换器用于将数字信号转换成模拟信号;所述信号隔离器用于将控制***输出信号和电气伺服模块隔离开;所述襟翼作动器驱动电路产生驱动襟翼作动器的电信号;所述襟翼作动器根据襟翼作动器驱动电路输出信号调节襟翼使其产生不同的摆角。
如图2所示一种大型风力机叶片主动降载控制方法,具体包括以下步骤:
步骤1:对模糊控制器FCi、前馈控制器FBi和PID控制器PIDi进行初始化,i=1,2,3,4。
步骤2:读取当前风速和叶片根部的应力值,将得到的叶片根部应力值y(k)与叶片根部额定应力值r(k)进行差值运算,得到应力偏差e(k)及偏差变化率ec(k),其中叶片额定应力值r(k)是由风机出厂前实验测得;
步骤21:将步骤2得到的应力偏差e(k)、偏差变化率ec(k)作为模糊控制器FCi的输入变量;
步骤22:选择隶属度函数进行模糊化,并依据模糊规则得到襟翼的控制量,反模糊化后求得模糊控制器FCi的输出变量,此输出为襟翼控制期望角θ1i;
步骤23:将步骤21的应力偏差e(k)作为PID控制器PIDi的输入信号,利用教与学寻优算法对PIDi参数KPi,KIi,KDi进行在线自整定,所述PID控制器PIDi的输出变量为襟翼控制期望角θ2i;
如图3所示,用于PID参数的教与学寻优算法包括以下步骤:
步骤1):设置初始参数,搜索区域范围定义为X=(x1,x2,…,xd)∈[L,U],L=(L1,L2,…,Ld)是空间下界,U=(U1,U2,…,Ud)是空间上界,d为优化问题的维数,d维空间第i个学员定义为学员规模为N,最大迭代次数为maxgen;
步骤2):教师的教学阶段:
步骤21):计算每个学生的适应值,选择最好个体作为老师Xteacher,计算个体平均值然后根据学员与个体平均水平的差异进行学习,如下式:
TFi=2-gen/maxgen (2)
式中:和分别表示第i个学员学习前和学习后的值;w1=1-gen/maxgen为自适应权系数;ri为0-1之间的随机数;TFi为1-2之间的某个数,其值随迭代次数的变化而变化;ri和TFi用于调整学习速率;gen与maxgen分别为当前迭代次数与最大迭代次数;
步骤22):学员更新:
如果的适应值比的适应值好,那么用代替否则,继续使用
步骤3):学员之间相互学习阶段:
步骤31):每个学员Xi在班级中随机选取一个学习对象Xj(j≠i),Xi通过分析自己和学员Xj之间的差异进行学习调整,如下式:
若Xi优于Xj,
若Xj优于Xi,
式中:w2=1-gen/maxgen为自适应权系数;ri为0-1之间的随机数;
步骤32):学员更新:
如果的适应值比的适应值好,那么用代替否则,继续使用
步骤4):根据适应度函数计算每个学员的适应度值,其公式如下:
式中,e(t)为***误差,根据适应度函数更新学员的全局最优值,当计算所得到最优值达到设定值或算法达到最大迭代次数时,退出教与学寻优算法,否则返回步骤2);
步骤24:采集风力发电机组的机前风速v(t),将其作为自变量,将襟翼摆动角θi作为因变量,对襟翼角度-风速进行拟合,建立襟翼角度-风速的前馈控制器FBi模型:θi(v)=a0+a1v+a2v2+L+anvn,采用最小二乘法确定各项系数;将风速作为前馈控制器FBi的输入信号,则前馈控制器FBi的输出变量为襟翼控制期望角θ3i。
步骤3:将步骤22与步骤23得到的襟翼期望角控制信号分别送往PLC进行处理,在PLC中设置切换算法:当叶根应力误差大于设定值时使用来自模糊控制器FCi的信号θ1i,来自PID控制器PIDi的控制信号θ2i将不起作用;当叶根应力误差小于设定值时使用来自PID控制器PIDi的控制信号θ2i,来自模糊控制器FCi的信号θ1i将不起作用;最后,PLC再耦合来自前馈控制器FBi的控制信号θ3i,并将这些信号传输给对应的襟翼作动器驱动电路。
步骤4:4个襟翼作动器分别接受来自4个襟翼作动器驱动电路的信号,执行襟翼摆动动作以减小叶片根部应力。
上述步骤2-4反复运行,使叶片根部应力值y(k)等于或接近风力机出厂前实验测得的叶片根部额定应力值r(k),直至完成控制任务。
5WM参考风力机在风况为11.4m/s的湍流风下,以叶片的叶根弯矩作为风力机降载目标,采用本发明一种大型风力机叶片主动降载控制方法与现有方法的对比结果如图4所示,可见采用本发明的方法有效降低了叶片根部的应变力。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何形式上的限制,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的简单修改、等同变化,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (4)
1.一种大型风力机叶片主动降载控制***,其特征在于,所述***包括信号采集模块、控制模块和电气伺服模块;
所述信号采集模块包括光纤应变传感器、机前风场信号传感设备和光纤应变信号处理设备;所述控制模块包括5路低通滤波器,5路模数转换器,4路控制单元,PLC,4路数模转换器和4路信号隔离器;所述电气伺服模块包括4路襟翼作动器驱动电路和4路襟翼作动器;
所述光纤应变传感器安装在风力机叶片根部,并与光纤应变信号处理设备连接;光纤应变信号处理设备的1路信号输出和机前风场信号传感设备对应4路襟翼的4路信号输出分别对应1路低通滤波器和1路模数转换器顺次连接;与光纤应变信号处理设备对应的1路模数转换器分别连接至4路控制单元,与机前风场信号传感设备对应的4路模数转换器一一对应连接至4路控制单元,每路控制单元通过PLC分别对应1路数模转换器、1路信号隔离器、1路襟翼作动器驱动电路和1路襟翼作动器顺次连接;
所述光纤应变传感器用于采集叶片根部的应变力信号,所述光纤应变信号处理设备用于将光纤应变传感器采集的信号转换为电压信号;所述机前风场信号传感设备用于测量风力机前风速;光纤应变信号处理设备和机前风场信号传感设备分别将信号传送至对应的低通滤波器,所述低通滤波器用于滤掉高频干扰信号;所述模数转换器用于将模拟信号转换成数字信号;所述控制单元包括对机前风场信号进行处理的前馈控制器,使用教与学算法进行参数寻优的PID控制器和对叶片根部的应变力信号进行处理的模糊控制器,用于减小叶片根部应力的控制运算;所述PLC用于切换PID控制器和模糊控制器的输出信号,并耦合来自前馈控制器的控制信号;所述数模转换器用于将数字信号转换成模拟信号;所述信号隔离器用于将控制***输出信号和电气伺服模块隔离开;所述襟翼作动器驱动电路产生驱动襟翼作动器的电信号;所述襟翼作动器根据襟翼作动器驱动电路输出信号调节襟翼使其产生不同的摆角。
2.根据权利要求1所述一种大型风力机叶片主动降载控制***,其特征在于,所述前馈控制器是在风力发电机组运行时,通过实时监测机前风速变化,计算出因降低随机风或湍流风引起的不均匀载荷所需的控制量。
3.权利要求1-2任一权利要求所述一种大型风力机叶片主动降载控制***的控制方法,其特征在于,具体包括以下步骤:
步骤1:对模糊控制器FCi、前馈控制器FBi和PID控制器PIDi进行初始化,i=1,2,3,4;
步骤2:读取当前风速和叶片根部的应力值,将得到的叶片根部应力值y(k)与叶片根部额定应力值r(k)进行差值运算,得到应力偏差e(k)及偏差变化率ec(k),其中叶片额定应力值r(k)是由风机出厂前实验测得;
步骤21:将步骤2得到的应力偏差e(k)、偏差变化率ec(k)作为模糊控制器FCi的输入变量;
步骤22:选择隶属度函数进行模糊化,并依据模糊规则得到襟翼的控制量,反模糊化后求得模糊控制器FCi的输出变量,此输出为襟翼控制期望角θ1i;
步骤23:将步骤21的应力偏差e(k)作为PID控制器PIDi的输入信号,利用教与学寻优算法对PIDi参数KPi,KIi,KDi进行在线自整定,所述PID控制器PIDi的输出变量为襟翼控制期望角θ2i;
步骤24:采集风力发电机组的机前风速v(t),将其作为自变量,将襟翼摆动角θi作为因变量,对襟翼角度-风速进行拟合,建立襟翼角度-风速的前馈控制器FBi模型:θi(v)=a0+a1v+a2v2+L+anvn,采用最小二乘法确定各项系数;将风速作为前馈控制器FBi的输入信号,则前馈控制器FBi的输出变量为襟翼控制期望角θ3i;
步骤3:将步骤22与步骤23得到的襟翼期望角控制信号分别送往PLC进行处理,在PLC中设置切换算法:当叶根应力误差大于设定值时使用来自模糊控制器FCi的信号θ1i,来自PID控制器PIDi的控制信号θ2i将不起作用;当叶根应力误差小于设定值时使用来自PID控制器PIDi的控制信号θ2i,来自模糊控制器FCi的信号θ1i将不起作用;最后,PLC再耦合来自前馈控制器FBi的控制信号θ3i,并将这些信号传输给对应的襟翼作动器驱动电路;
步骤4:4个襟翼作动器分别接受来自4个襟翼作动器驱动电路的信号,执行襟翼摆动动作以减小叶片根部应力;
上述步骤2-4反复运行,使叶片根部应力值y(k)等于或接近风力机出厂前实验测得的叶片根部额定应力值r(k),直至完成控制任务。
4.根据权利要求3所述一种控制方法,其特征在于,用于PID参数的教与学寻优算法包括以下步骤:
步骤1):设置初始参数,搜索区域范围定义为X=(x1,x2,…,xd)∈[L,U],L=(L1,L2,…,Ld)是空间下界,U=(U1,U2,…,Ud)是空间上界,d为优化问题的维数,d维空间第i个学员定义为学员规模为N,最大迭代次数为maxgen;
步骤2):教师的教学阶段:
步骤21):计算每个学生的适应值,选择最好个体作为老师Xteacher,计算个体平均值然后根据学员与个体平均水平的差异进行学习,如下式:
TFi=2-gen/maxgen (2)
式中:和分别表示第i个学员学习前和学习后的值;w1=1-gen/maxgen为自适应权系数;ri为0-1之间的随机数;TFi为1-2之间的某个数,其值随迭代次数的变化而变化;ri和TFi用于调整学习速率;gen与maxgen代次数与最大迭代次数;
步骤22):学员更新:
如果的适应值比的适应值好,那么用代替否则,继续使用
步骤3):学员之间相互学习阶段:
步骤31):每个学员Xi在班级中随机选取一个学习对象Xj(j≠i),Xi通过分析自己和学员Xj之间的差异进行学习调整,如下式:
若Xi优于Xj,
若Xj优于Xi,
式中:w2=1-gen/maxgen为自适应权系数;ri为0-1之间的随机数;
步骤32):学员更新:
如果的适应值比的适应值好,那么用代替否则,继续使用
步骤4):根据适应度函数计算每个学员的适应度值,其公式如下:
式中,e(t)为***误差,根据适应度函数更新学员的全局最优值,当计算所得到最优值达到设定值或算法达到最大迭代次数时,退出教与学寻优算法,否则返回步骤2)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610274672.0A CN105888971B (zh) | 2016-04-28 | 2016-04-28 | 一种大型风力机叶片主动降载控制***及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610274672.0A CN105888971B (zh) | 2016-04-28 | 2016-04-28 | 一种大型风力机叶片主动降载控制***及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105888971A CN105888971A (zh) | 2016-08-24 |
CN105888971B true CN105888971B (zh) | 2018-12-18 |
Family
ID=56702436
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610274672.0A Expired - Fee Related CN105888971B (zh) | 2016-04-28 | 2016-04-28 | 一种大型风力机叶片主动降载控制***及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105888971B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108241288A (zh) * | 2016-12-23 | 2018-07-03 | 重庆邮电大学 | 一种基于改进tlbo算法的pid控制器参数整定方法 |
CN107559143B (zh) * | 2017-07-27 | 2019-08-27 | 华北电力大学 | 一种大型风力机尾缘襟翼结构参数寻优及多目标襟翼优化控制方法 |
CN107762730B (zh) * | 2017-08-23 | 2019-06-18 | 华北电力大学 | 一种带有尾缘襟翼的大型变桨风力机控制***及控制方法 |
CN109958574B (zh) * | 2017-12-22 | 2020-05-05 | 北京金风科创风电设备有限公司 | 用于风力发电机组的降载控制方法和装置 |
CN109162878B (zh) * | 2018-08-01 | 2020-06-30 | 汕头大学 | 一种风力发电机智能叶片测控方法及其*** |
EP3667076A1 (en) * | 2018-12-13 | 2020-06-17 | Siemens Gamesa Renewable Energy A/S | Estimating wind speed |
CN111749845B (zh) * | 2019-03-29 | 2022-04-05 | 北京金风科创风电设备有限公司 | 用于风电机组的降载控制方法和装置 |
CN115242072B (zh) * | 2022-09-16 | 2022-12-23 | 深圳市首航新能源股份有限公司 | 一种降低逆变管应力的方法、装置及逆变器 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001195102A (ja) * | 2000-01-13 | 2001-07-19 | Seiko Seiki Co Ltd | 可変飽和機能を有する制御装置 |
CN102192084A (zh) * | 2011-06-22 | 2011-09-21 | 成都阜特科技有限公司 | 风力发电机组变桨用伺服驱动控制*** |
CN102878016A (zh) * | 2012-10-26 | 2013-01-16 | 华北电力大学 | 基于智能结构的大型风力机叶片摆振抑制***及控制方法 |
CN102900609A (zh) * | 2012-10-26 | 2013-01-30 | 华北电力大学 | 超磁致伸缩襟翼结构的风力机叶片减振***及控制方法 |
CN103912447A (zh) * | 2014-04-04 | 2014-07-09 | 上海电机学院 | 基于攻角权系数的模糊pid控制的独立变桨距控制*** |
CN104005909A (zh) * | 2014-04-22 | 2014-08-27 | 重庆邮电大学 | 非线性前馈与模糊pid结合的风力发电机组变桨距控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6600240B2 (en) * | 1997-08-08 | 2003-07-29 | General Electric Company | Variable speed wind turbine generator |
-
2016
- 2016-04-28 CN CN201610274672.0A patent/CN105888971B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001195102A (ja) * | 2000-01-13 | 2001-07-19 | Seiko Seiki Co Ltd | 可変飽和機能を有する制御装置 |
CN102192084A (zh) * | 2011-06-22 | 2011-09-21 | 成都阜特科技有限公司 | 风力发电机组变桨用伺服驱动控制*** |
CN102878016A (zh) * | 2012-10-26 | 2013-01-16 | 华北电力大学 | 基于智能结构的大型风力机叶片摆振抑制***及控制方法 |
CN102900609A (zh) * | 2012-10-26 | 2013-01-30 | 华北电力大学 | 超磁致伸缩襟翼结构的风力机叶片减振***及控制方法 |
CN103912447A (zh) * | 2014-04-04 | 2014-07-09 | 上海电机学院 | 基于攻角权系数的模糊pid控制的独立变桨距控制*** |
CN104005909A (zh) * | 2014-04-22 | 2014-08-27 | 重庆邮电大学 | 非线性前馈与模糊pid结合的风力发电机组变桨距控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105888971A (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105888971B (zh) | 一种大型风力机叶片主动降载控制***及方法 | |
CN102354107B (zh) | 一种交流位置伺服***模型参数在线辨识和控制方法 | |
Bottasso et al. | Aero-servo-elastic modeling and control of wind turbines using finite-element multibody procedures | |
CN104595106B (zh) | 基于强化学习补偿的风力发电变桨距控制方法 | |
CN107203137B (zh) | 非线性抽水蓄能机组调速***增益启发式自适应pid控制方法 | |
Butt et al. | Multi-variable integral sliding mode control of a two degrees of freedom helicopter | |
CN107505841B (zh) | 一种基于干扰估计器的机械臂姿态鲁棒控制方法 | |
Iyasere et al. | Optimum seeking-based non-linear controller to maximise energy capture in a variable speed wind turbine | |
CN102411305A (zh) | 单旋翼直升机/涡轴发动机综合抗扰控制***设计方法 | |
CN101917150A (zh) | 永磁同步电机模糊神经网络广义逆鲁棒控制器及构造方法 | |
CN112696319A (zh) | 利用准确在线模型的基于风力涡轮模型的控制和估计 | |
WO2018115423A1 (en) | Fatigue load minimization in an operation of a wind farm | |
CN111781835A (zh) | 一种镇定二阶惯性加纯滞后***的线性自抗扰控制器设计方法 | |
CN102410138B (zh) | 一种获取风力发电机组最优控制输入的方法 | |
CN112523944A (zh) | 一种风力发电机变桨***自适应动态面控制方法 | |
CN109611274B (zh) | 一种风力发电机组高风速区lqg优化控制方法 | |
Maheedhar et al. | A behavioral study of different controllers and algorithms in real-time applications | |
CN109884881A (zh) | 一种基于非线性pid控制技术的稳瞄伺服控制器的设计 | |
CN111538232B (zh) | 基于自适应神经模糊控制的无人行车防摇定位方法及*** | |
CN110018695B (zh) | 超大跨度悬索桥颤振气动翼板主动控制方法 | |
CN106499583B (zh) | 基于rbf神经网络技术的风力发电机组***辨识方法 | |
CN108155836A (zh) | 基于全局在线启发式动态规划永磁同步电机矢量控制方法 | |
Chai et al. | Design of crane anti-swing controller based on differential flat and linear active disturbance rejection control | |
CN112947606A (zh) | Bp神经网络pid预测控制的锅炉液位控制***及方法 | |
Fu et al. | Design of adrc for unmanned surface vehicle heading based on immune particle swarm optimization |
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: 20181218 |
|
CF01 | Termination of patent right due to non-payment of annual fee |