CN109505742B - 一种确定非常规风力机推力系数的方法 - Google Patents

一种确定非常规风力机推力系数的方法 Download PDF

Info

Publication number
CN109505742B
CN109505742B CN201811569307.8A CN201811569307A CN109505742B CN 109505742 B CN109505742 B CN 109505742B CN 201811569307 A CN201811569307 A CN 201811569307A CN 109505742 B CN109505742 B CN 109505742B
Authority
CN
China
Prior art keywords
wind
speed
coefficient
wind turbine
wind speed
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
CN201811569307.8A
Other languages
English (en)
Other versions
CN109505742A (zh
Inventor
朱建勇
庄林学
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.)
West New Energy Technology Hebei Co ltd
Original Assignee
Shenyang Aerospace University
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 Shenyang Aerospace University filed Critical Shenyang Aerospace University
Priority to CN201811569307.8A priority Critical patent/CN109505742B/zh
Publication of CN109505742A publication Critical patent/CN109505742A/zh
Application granted granted Critical
Publication of CN109505742B publication Critical patent/CN109505742B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
    • 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
    • F03D3/00—Wind motors with rotation axis substantially perpendicular to the air flow entering the rotor 
    • F03D3/005—Wind motors with rotation axis substantially perpendicular to the air flow entering the rotor  the axis being vertical
    • 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
    • F03D3/00—Wind motors with rotation axis substantially perpendicular to the air flow entering the rotor 
    • F03D3/06—Rotors
    • F03D3/061—Rotors characterised by their aerodynamic shape, e.g. aerofoil profiles
    • 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
    • F03D3/00—Wind motors with rotation axis substantially perpendicular to the air flow entering the rotor 
    • F03D3/06—Rotors
    • F03D3/062—Rotors characterised by their construction elements
    • 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/74—Wind turbines with rotation axis perpendicular to the 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

一种确定非常规风力机推力系数的方法，步骤为：设定风力机启动风速、切入风速、额定风速及切出风速；在启动风速下，通过数值模拟获取不同尖速比下的扭矩系数和推力系数，将扭矩系数为零时所对应的推力系数作为启动风速下的推力系数；在切入风速与额定风速之间时，通过数值模拟获取额定风速时不同尖速比下的风能利用系数和推力系数，取最大风能利用系数对应尖速比下的推力系数作为切入风速与额定风速之间的推力系数；在额定风速与切出风速之间时，保持额定功率，通过推力系数随尖速比的变化曲线得到相应尖速比下的推力系数，同时通过流动控制手段诱导产生附加推力系数；当超过切出风速时，风力机停机，拟合出风力机在不同风速下的推力系数曲线。

Description

一种确定非常规风力机推力系数的方法
技术领域
本发明属于风力机气动特性及结构设计技术领域，特别是涉及一种确定非常规风力机推力系数的方法。
背景技术
环境和能源是人类赖以生存的基础，化石能源的巨量消耗会导致严重的环境污染问题，只有积极开发新能源和可再生能源，才能更好的缓解能源和环境问题。
风能是一种清洁无污染的可再生能源，而风力机则是将风能转换为机械能的设备，风力机气动性能的优劣直接影响风能转换效率，而风力机的推力系数又是表征风力机气动性能指标之一，其表示的是风力机所受推力的大小，推力系数在风力机的基础设计、结构设计等方面具有重要作用。
目前，对于常规的水平轴风力机来说，可以利用现有的动量-叶素理论来确定风力机的推力系数；对于常规的垂直轴达里厄风力机来说，可以利用现有的流管法确定风力机的推力系数。
然而，随着风能技术的日益发展，涌现出各种新颖型式的非常规风力机，现阶段，对于非常规风力机尚没有成熟的确定推力系数的方法。
发明内容
针对现有技术存在的问题，本发明提供一种确定非常规风力机推力系数的方法，适用于各种型号的非常规风力机，能够获取不同风速下的非常规风力机的推力系数，进而为风力机的结构设计提供必要的气动数据。
为了实现上述目的，本发明采用如下技术方案：一种确定非常规风力机推力系数的方法，包括如下步骤：
步骤一：根据风力机的设计指标，先对风力机的启动风速、切入风速、额定风速及切出风速进行设定；
步骤二：在设定的启动风速下，通过数值模拟方式获取不同尖速比下的扭矩系数和推力系数，并将扭矩系数为零时所对应的推力系数作为启动风速下的推力系数；
步骤三：当风速处于设定的切入风速与额定风速之间时，通过数值模拟方式获取额定风速时不同尖速比下的风能利用系数和推力系数，并取最大风能利用系数对应尖速比下的推力系数作为切入风速与额定风速之间的推力系数；
步骤四：当风速超出设定的额定风速且处于设定的额定风速与切出风速之间时，随着风速的增大，为了使风力机功率保持在额定功率，风力机的尖速比逐渐减小，通过推力系数随尖速比的变化曲线得到相应尖速比下的推力系数；同时，为了使风力机功率保持在额定功率，通过流动控制手段对风力机功率进行调节，且流动控制手段会诱导产生风力机的附加推力系数；
步骤五：当风速超过设定的切出风速时，风力机停机；
步骤六：根据步骤三至步骤五获取的推力系数数据，拟合出风力机在不同风速下的推力系数曲线。
采用的数值模拟方式为风力机计算流体力学数值模拟方法，计算域采用结构网格，并通过滑移网格技术模拟风力机的旋转运动；采用计算流体力学软件求解非定常不可压缩的时均N-S方程，湍流方程标准k-ε二方程模型，采用有限体积法离散上述方程，对流项采用一阶迎风格式，扩散项采用中心差分格式，通过SIMPLE算法耦合压力与速度；边界条件采用速度入口边界、压力出口以及壁面无滑移边界。
风力机的扭矩系数、风能利用系数、推力系数及尖速比的计算公式为
CM＝2M/ρAV2R
Cp＝2P/ρAV3
CT＝2FT/ρAV2
λ＝ωR/V
其中，P＝2πnM/60，式中，P为风力机功率，n为风力机转速，M为风力机扭矩，CM为扭矩系数，ρ为空气密度，A为风力机风轮扫掠面积，V为来流风速，R为风力机风轮旋转半径，CP为风能利用系数，CT为推力系数，FT为风力机推力，λ为尖速比，ω为风力机旋转角度速。
流动控制手段诱导产生的附加推力系数的计算公式为
CT(附加)＝2FT(附加)/ρAV2
其中，FT(附加)＝(P-P(额定))η/V，式中，FT(附加)为附加阻力，P为风力机功率，P(额定)为风力机额定功率，η为修正系数，V为来流风速，CT(附加)为附加推力系数，ρ为空气密度，A为风力机风轮扫掠面积。
本发明的有益效果：
本发明与现有技术相比，适用于各种型号的非常规风力机，能够获取不同风速下的非常规风力机的推力系数，进而为风力机的结构设计提供必要的气动数据。
附图说明
图1为采用扰流板控制的100kW的垂直轴H型风力机的结构示意图；
图2为图1所示风力机的风能利用系数随尖速比变化的曲线图；
图3为图1所示风力机的推力系数随尖速比变化的曲线图；
图4为图1所示风力机在不同风速下的推力系数曲线图；
图中，1—叶片，2—横臂，3—扰流板。
具体实施方式
下面结合附图和具体实施例对本发明做进一步的详细说明。
本实施例中，需要确定推力系数的风力机如图1所示，具体为采用扰流板控制的100kW的垂直轴H型风力机，风力机的旋转半径为16m，叶片高度为28m。
一种确定非常规风力机推力系数的方法，包括如下步骤：
步骤一：根据风力机的设计指标，先对风力机的启动风速、切入风速、额定风速及切出风速进行设定；本实施例中，启动风速设定为2m/s，切入风速设定为5m/s，额定风速设定为9m/s，切出风速设定为25m/s；
步骤二：在设定的启动风速下，通过数值模拟方式获取不同尖速比下的扭矩系数和推力系数，并将扭矩系数为零时所对应的推力系数作为启动风速下的推力系数；本实施例中，推力系数设为0.984；
步骤三：当风速处于设定的切入风速与额定风速之间时，通过数值模拟方式获取额定风速时不同尖速比下的风能利用系数和推力系数，并取最大风能利用系数对应尖速比下的推力系数作为切入风速与额定风速之间的推力系数；本实施例中，数值模拟的额定风速时不同尖速比下的风能利用系数见图2，数值模拟的额定风速时不同尖速比下的推力系数见图3；在“最大风能利用系数对应尖速比下的推力系数”中，最大风能利用系数为0.3012，尖速比为3.5，推力系数为0.634；
步骤四：当风速超出设定的额定风速且处于设定的额定风速与切出风速之间时，随着风速的增大，为了使风力机功率保持在100kW的额定功率，风力机的尖速比逐渐减小，通过图3所示的推力系数随尖速比的变化曲线得到相应尖速比下的推力系数；同时，为了使风力机功率保持在100kW的额定功率，通过流动控制手段对风力机功率进行调节，且流动控制手段会诱导产生风力机的附加推力系数；
步骤五：当风速超过设定的切出风速时，风力机停机；
步骤六：根据步骤三至步骤五获取的推力系数数据，拟合出风力机在不同风速下的推力系数曲线，具体如图4所示。
采用的数值模拟方式为风力机计算流体力学数值模拟方法，计算域采用结构网格，并通过滑移网格技术模拟风力机的旋转运动；采用计算流体力学软件求解非定常不可压缩的时均N-S方程，湍流方程标准k-ε二方程模型，采用有限体积法离散上述方程，对流项采用一阶迎风格式，扩散项采用中心差分格式，通过SIMPLE算法耦合压力与速度；边界条件采用速度入口边界、压力出口以及壁面无滑移边界。
风力机的扭矩系数、风能利用系数、推力系数及尖速比的计算公式为
CM＝2M/ρAV2R
Cp＝2P/ρAV3
CT＝2FT/ρAV2
λ＝ωR/V
其中，P＝2πnM/60，式中，P为风力机功率，n为风力机转速，M为风力机扭矩，CM为扭矩系数，ρ为空气密度，A为风力机风轮扫掠面积，V为来流风速，R为风力机风轮旋转半径，CP为风能利用系数，CT为推力系数，FT为风力机推力，λ为尖速比，ω为风力机旋转角度速。
流动控制手段诱导产生的附加推力系数的计算公式为
CT(附加)＝2FT(附加)/ρAV2
其中，FT(附加)＝(P-P(额定))η/V，式中，FT(附加)为附加阻力，P为风力机功率，P(额定)为风力机额定功率，η为修正系数，V为来流风速，CT(附加)为附加推力系数，ρ为空气密度，A为风力机风轮扫掠面积；本实施例中，修正系数η取0.95。
实施例中的方案并非用以限制本发明的专利保护范围，凡未脱离本发明所为的等效实施或变更，均包含于本案的专利范围中。

Claims (3)

1.一种确定风力机推力系数的方法，其特征在于包括如下步骤：
步骤一：根据风力机的设计指标，先对风力机的启动风速、切入风速、额定风速及切出风速进行设定；
步骤二：在设定的启动风速下，通过数值模拟方式获取不同尖速比下的扭矩系数和推力系数，并将扭矩系数为零时所对应的推力系数作为启动风速下的推力系数；
步骤三：当风速处于设定的切入风速与额定风速之间时，通过数值模拟方式获取额定风速时不同尖速比下的风能利用系数和推力系数，并取最大风能利用系数对应尖速比下的推力系数作为切入风速与额定风速之间的推力系数；
步骤四：当风速超出设定的额定风速且处于设定的额定风速与切出风速之间时，随着风速的增大，为了使风力机功率保持在额定功率，风力机的尖速比逐渐减小，通过推力系数随尖速比的变化曲线得到相应尖速比下的推力系数；同时，为了使风力机功率保持在额定功率，通过流动控制手段对风力机功率进行调节，且流动控制手段会诱导产生风力机的附加推力系数；最后，将得到的相应尖速比下的推力系数和附加推力系数进行求和；其中，附加推力系数的计算公式为
Figure FDA0002445704330000011
其中，
Figure FDA0002445704330000012
式中，
Figure FDA0002445704330000013
为附加阻力，P为风力机功率，P(额定)为风力机额定功率，η为修正系数，V为来流风速，CT(附加)为附加推力系数，ρ为空气密度，A为风力机风轮扫掠面积；
步骤五：根据步骤三至步骤四获取的推力系数数据，拟合出风力机在不同风速下的推力系数曲线图。
2.根据权利要求1所述的一种确定风力机推力系数的方法，其特征在于：采用的数值模拟方式为风力机计算流体力学数值模拟方法，计算域采用结构网格，并通过滑移网格技术模拟风力机的旋转运动；采用计算流体力学软件求解非定常不可压缩的时均N-S方程，湍流方程标准k-ε二方程模型，采用有限体积法离散上述N-S方程和湍流方程，对流项采用一阶迎风格式，扩散项采用中心差分格式，通过SIMPLE算法耦合压力与速度；边界条件采用速度入口边界、压力出口以及壁面无滑移边界。
3.根据权利要求1所述的一种确定风力机推力系数的方法，其特征在于：风力机的扭矩系数、风能利用系数、推力系数及尖速比的计算公式为
CM＝2M/ρAV2R
Cp＝2P/ρAV3
CT＝2FT/ρAV2
λ＝ωR/V
其中，P＝2πnM/60，式中，P为风力机功率，n为风力机转速，M为风力机扭矩，CM为扭矩系数，ρ为空气密度，A为风力机风轮扫掠面积，V为来流风速，R为风力机风轮旋转半径，CP为风能利用系数，CT为推力系数，FT为风力机推力，λ为尖速比，ω为风力机旋转角度速。
CN201811569307.8A 2018-12-21 2018-12-21 一种确定非常规风力机推力系数的方法 Active CN109505742B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811569307.8A CN109505742B (zh) 2018-12-21 2018-12-21 一种确定非常规风力机推力系数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811569307.8A CN109505742B (zh) 2018-12-21 2018-12-21 一种确定非常规风力机推力系数的方法

Publications (2)

Publication Number Publication Date
CN109505742A CN109505742A (zh) 2019-03-22
CN109505742B true CN109505742B (zh) 2020-06-30

Family

ID=65754089

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811569307.8A Active CN109505742B (zh) 2018-12-21 2018-12-21 一种确定非常规风力机推力系数的方法

Country Status (1)

Country Link
CN (1) CN109505742B (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2476509A (en) * 2009-12-24 2011-06-29 Rolls Royce Plc Turbine with reduced thrust coefficient at excessive speed
CN106812658B (zh) * 2015-11-27 2019-09-06 中国船舶重工集团海装风电股份有限公司 一种风力发电机组的控制方法及装置
KR101734757B1 (ko) * 2016-03-29 2017-05-11 한밭대학교 산학협력단 윈드쉬어에 의한 풍력터빈 블레이드의 추력 모델링 방법 및 시스템
CN106224162B (zh) * 2016-07-29 2018-10-16 电子科技大学 风电机组的载荷模型建立方法及载荷控制方法

Also Published As

Publication number Publication date
CN109505742A (zh) 2019-03-22

Similar Documents

Publication Publication Date Title
Castelli et al. Effect of blade number on a straight-bladed vertical-axis Darreius wind turbine
CN103244348B (zh) 变速变桨风力发电机组功率曲线优化方法
CN106777499B (zh) 一种双馈异步风力发电机组的整机动态建模方法
CN103617308B (zh) 一种风电场频域等效模型的构建方法
CN107045574B (zh) 基于svr的风力发电机组低风速段有效风速估计方法
Islam et al. Design and simulation of a small wind turbine blade with qblade and validation with MATLAB
CN106294959B (zh) 模型参考自适应控制与双馈风机降阶模型相结合的建模仿真方法
CN111502914A (zh) 基于线性变参数系统的风机变桨控制器设计方法
CN108757306B (zh) 一种风力发电机控制方法、装置以及风力发电机
CN109505742B (zh) 一种确定非常规风力机推力系数的方法
CN105406668A (zh) 变功率盘式发电机
CN112682258B (zh) 基于Backstepping的大型风力机最大功率点跟踪控制方法
CN109611268A (zh) 一种双叶轮水平轴风力机设计优化方法
CN109657380A (zh) 一种基于扩展卡尔曼滤波的双馈风电场动态等值方法
Rocchio et al. Development of a BEM-CFD tool for Vertical Axis Wind Turbines based on the Actuator Disk model
CN106837683B (zh) 迎风面的出流切线倾角的最优值确定方法
Nilsson Numerical computations of wind turbine wakes and wake interaction: Optimization and control
Mao et al. Theoretical studies and simulation on the wakes of wind turbines
Christ et al. Modelling of a wind power turbine
Zhou et al. Analysis of extreme operating gust influence on aerodynamic performance of wind turbine
Zheng et al. Design and Performance Analysis of Distributed Equal Angle Spiral Vertical Axis Wind Turbine
CN201461254U (zh) 垂直风力发电装置
Miguel Airfoil Improvement on Horizontal Axis Wind Turbine Suitable for Local Construction in Underdeveloped Countries
Barbier et al. Aerodynamic Study of a Small-Scale Wind Turbine
Prananda et al. Performance Analysis of Screw Turbine Design with additional flaps modification using computational fluid dynamics method

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
TR01 Transfer of patent right

Effective date of registration: 20230616

Address after: 061000 West Juxian Street and South Wenhui Road, Zoje Industrial Park, Bohai New Area, Cangzhou City, Hebei Province (in the yard of Cangzhou Lingang Xincheng Plastic Products Co., Ltd.)

Patentee after: West New Energy Technology (Hebei) Co.,Ltd.

Address before: 110136, Liaoning, Shenyang moral and Economic Development Zone, No. 37 South Avenue moral

Patentee before: SHENYANG AEROSPACE University

TR01 Transfer of patent right