CN106762409B - 一种基于β锥角的风力机叶素动量修正方法 - Google Patents
一种基于β锥角的风力机叶素动量修正方法 Download PDFInfo
- Publication number
- CN106762409B CN106762409B CN201611034031.4A CN201611034031A CN106762409B CN 106762409 B CN106762409 B CN 106762409B CN 201611034031 A CN201611034031 A CN 201611034031A CN 106762409 B CN106762409 B CN 106762409B
- Authority
- CN
- China
- Prior art keywords
- angle
- wind
- cone angle
- amendment
- conversion system
- 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
Links
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 11
- 238000002715 modification method Methods 0.000 title claims abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 12
- 238000009795 derivation Methods 0.000 claims abstract description 4
- 230000001939 inductive effect Effects 0.000 claims description 8
- 238000012804 iterative process Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011144 upstream manufacturing 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
-
- 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
- F03D1/00—Wind motors with rotation axis substantially parallel to the air flow entering the rotor
- F03D1/06—Rotors
-
- 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)
- Wind Motors (AREA)
Abstract
一种基于β锥角的风力机叶素动量修正方法，涉及一种风机叶轮动量修正方法，包括以下步骤：对叶素动量理论推导过程中的每一表达式进行相应修正，最终得到修正后风能利用系数表达式，如下；通过有限元分析软件FLUENT计算模型的发散角，并由图3确定在该发散角的条件下能达到的最大风能利用系数及对应的风轮锥角。该方法提出一种基于β锥角的针对水平轴风力发电机组的传统叶素动量的修正方法。该方法引入风轮锥角β，计算出锥角对风力机风能利用系数大小的影响，通过优化风轮锥角进一步提高风力机风能利用率。
Description
技术领域
本发明涉及涉及一种风机叶轮动量修正方法，特别是涉及一种基于β锥角的风力机叶素动量修正方法。
背景技术
目前风力发电技术多采用基于零锥角的水平轴风力机发电机组。其设计方法为，根据对功率的设计要求，选取合适的风力机叶片的翼型及尺寸，为追求最大风能利用系数C p，通过传统零锥角叶素动量理论对风轮的两个重要参数节距角θ及速度比γ进行优化计算，得到风力机的优化运行参数。传统零锥角叶素动量理论能够成功地优化零锥角风力机运行参数，但无法计算出风轮锥角β变化对风力机风能利用系数大小的影响。
发明内容
本发明的目的在于提供一种基于β锥角的风力机叶素动量修正方法，该方法提出一种基于β锥角的针对水平轴风力发电机组的传统叶素动量的修正方法。该方法引入风轮锥角β，计算出锥角对风力机风能利用系数大小的影响，通过优化风轮锥角进一步提高风力机风能利用率。
本发明的目的是通过以下技术方案实现的：
该修正方法按如下步骤进行：
步骤1：根据图1所示修正前后风轮气流模型图，对叶素动量理论推导过程中的每一表达式进行相应修正，最终得到修正后风能利用系数表达式，如下
步骤2：通过图2所示的修正后迭代流程图计算该表达式的的数值解。
步骤3：设定参数的取值范围，遍历各组参数，通过步骤2计算该组参数对应数值解。
步骤4：将具有相同参数的划分为一组，分别找出组内最大风能利用系数记录，组成随参数、变化的数据组，绘制图3所示随、变化的曲线图。
步骤5：通过有限元分析软件FLUENT计算模型的发散角，并由图3确定在该发散角的条件下能达到的最大风能利用系数及对应的风轮锥角。
本发明的优点与效果是：
本发明通过图3的随、变化曲线可看出，锥角为零时，最大风能利用系数=0.481该值与传统零锥角叶素动量叶轮计算结果一致，与发散角无关。当考虑锥角时，随锥角的增加先增大后减小，发散角为通常值12°时，最大风能利用系数=0.495，与0.481相比，提高了2.9%。
附图说明
图1(a)为修正前风轮气流模型；
图1(b)为修正后风轮气流模型；
图2为修正后迭代流程图；
图3为随、变化的曲线图。
具体实施方式
下面结合实施例对本发明进行详细说明。
针对步骤1，对叶素动量理论推导过程中的每一表达式进行修正的详细过程入下：
（a）传统零锥角叶素动量理论：
经典的叶素动量理论采用零锥角模型，即风轮旋转面为零度锥角的平面（图1（a）），此时气流的径向速度分量对风轮不产生动力作用，因此只考虑气流轴向速度对风轮产生的动力效果。由动量理论得到，气流作用在风轮根部距离处叶素的轴向推力和转矩分别为：
（1）
（2）
为空气密度；为风轮上游远端风速；为轴向诱导因子；为周向诱导因子。
同时，由叶素理论得到
（3）
（4）
为叶片数；为叶素弦长；以叶素为参照系的相对风速；为入流角，表示风轮处相对风速与风轮旋转平面的夹角，其表达式为
（5）
为速度比，其表达式为
（6）
、为风力机翼型升力、阻力系数，对于一个确定的翼型，、主要由攻角决定
（7）
为节距角，表示叶素弦长方向与风轮旋转平面之间的夹角。
分别将式（1）、（3）及式（2）、（4）建立等量关系，并引入普朗特叶尖损失因子，可计算出轴向诱导因子及周向诱导因子
（8）
（9）
其中为实度
（10）
进一步计算出风能利用系数
（11）
（b）修正后β锥角叶素动量理论模型：
β锥角叶素动量理论采用β锥角模型，即叶片与旋转轴垂直平面夹角为β（如图1（b）），风轮旋转面为圆锥面。此时，气流的径向速度分量亦会对风轮产生转矩作用，进而转换为风轮机的械能。为计算出这部分机械能，对经典零锥角叶素动量理论做如下修正。
与叶根距离为的风轮叶素，其距离旋转轴的垂直距离为
（12）
气流速度与风轮旋转锥面垂直的速度分量为
（13）
γ为该处气流流向与轴向夹角；为对轴向诱导因子的修正（本文中上角标c表示对相应量的修正值）。
将修正式（12）、（13）及β代入式（1）~（11），分别得到如下式（14）~（24）对应的修正表达式
（14）
（15）
（16）
（17）
（18）
（19）
（20）
（21）
（22）
（23）
（24）。
Claims (1)
1.一种基于β锥角的风力机叶素动量修正方法，其特征在于，所述方法包括以下步骤：
步骤1：对叶素动量理论推导过程中的每一表达式进行相应修正，最终得到修正后风能利用系数表达式，如下
；
式中，为修正后实度；为修正后速度比；为修正后轴向诱导因子；为修正后周向诱导因子；为修正后入流角；为修正后气流发散角；
步骤2：通过初始化、计算入流角、计算攻角、确定升阻系数，、计算诱导因子的迭代过程，计算该表达式的的数值解；
步骤3：设定参数的取值范围，遍历各组参数，通过步骤2计算该组参数对应数值解；
步骤4：将具有相同参数的划分为一组，分别找出组内最大风能利用系数记录，组成随参数变化的数据组，绘制随变化的曲线图；
步骤5：通过有限元分析软件FLUENT计算模型的发散角，并由步骤4中所绘制的曲线确定在该发散角的条件下能达到的最大风能利用系数及对应的风轮锥角。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611034031.4A CN106762409B (zh) | 2016-11-23 | 2016-11-23 | 一种基于β锥角的风力机叶素动量修正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611034031.4A CN106762409B (zh) | 2016-11-23 | 2016-11-23 | 一种基于β锥角的风力机叶素动量修正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106762409A CN106762409A (zh) | 2017-05-31 |
| CN106762409B true CN106762409B (zh) | 2019-04-02 |
Family
ID=58970851
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611034031.4A Active CN106762409B (zh) | 2016-11-23 | 2016-11-23 | 一种基于β锥角的风力机叶素动量修正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106762409B (zh) |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130195670A1 (en) * | 2010-04-26 | 2013-08-01 | Se Blades Technology B.V. | Rotor for a wind turbine |
| CN103207942B (zh) * | 2012-01-16 | 2016-08-17 | 北京能高自动化技术股份有限公司 | 基于动量-叶素理论的致动盘不均匀受力载荷计算方法 |
| ES2551659B1 (es) * | 2014-05-20 | 2016-09-14 | Acciona Windpower, S.A. | Método de control de un aerogenerador |
| CN205559157U (zh) * | 2016-04-12 | 2016-09-07 | 王伟亮 | 伞形分布自寻风风叶及自寻风风力发电机 |
-
2016
- 2016-11-23 CN CN201611034031.4A patent/CN106762409B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106762409A (zh) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106897486B (zh) | 考虑湍流强度影响的抛物线形风电机组尾流模型计算方法 | |
| CN107002636B (zh) | 用于估计风速，包括计算针对叶片扭转调节的桨距角的方法 | |
| CN102322407B (zh) | 风力机叶片气动外形协同设计方法 | |
| AU2020100825A4 (en) | Method for evaluating aerodynamic performance of wind wheel of wind turbine based on three-factor fitting integral method | |
| CN108488035B (zh) | 永磁直驱风力发电机组失速和变桨混合控制方法 | |
| CN113742861B (zh) | 一种适用于风力发电机风洞试验的叶片模型优化设计方法 | |
| Louw et al. | The design of an axial flow fan for application in large air-cooled heat exchangers | |
| CN109815624A (zh) | 一种考虑进气总压畸变影响的压气机稳定边界判断方法 | |
| CN109255184B (zh) | 一种风力机全尾流场风速分布的确定方法及系统 | |
| Clifton-Smith | Wind turbine blade optimisation with tip loss corrections | |
| Velázquez et al. | Design and experimentation of a 1 MW horizontal axis wind turbine | |
| CN106777525B (zh) | 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法 | |
| Wang et al. | Optimized chord and twist angle distributions of wind turbine blade considering Reynolds number effects | |
| Zheng et al. | Effect of blade number on performance of drag type vertical axis wind turbine | |
| CN106939902B (zh) | 节能型直壁前、后盘变曲率曲线元素三元叶轮及采用该叶轮的离心式风机 | |
| CN106919730B (zh) | 一种采用风速衰减因子的风电场尾流确定方法 | |
| CN106762409B (zh) | 一种基于β锥角的风力机叶素动量修正方法 | |
| CN109595040A (zh) | 一种燃气轮机叶片扭曲规律设计方法 | |
| CN117592388A (zh) | 一种基于cfd的风电场多机尾流模拟方法 | |
| CN109325274B (zh) | 基于三因子拟合积分法的风力机风轮气动性能评价方法 | |
| CN103207942A (zh) | 基于动量-叶素理论的致动盘不均匀受力载荷计算方法 | |
| CN106599381A (zh) | 通过调整风轮桨距角及转速提高风力机效率的方法 | |
| CN112761899B (zh) | 一种风力发电机前置风速管测量风向的方法 | |
| CN109635342A (zh) | 一种基于Fluent的风机尾流自动获取方法 | |
| CN106687684B (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 | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170531 Assignee: Shenyang Hetai Technology Co., Ltd Assignor: Shenyang Univ. Contract record no.: X2019210000025 Denomination of invention: Wind turbine blade element momentum correction method based on beta taper angle Granted publication date: 20190402 License type: Exclusive License Record date: 20191030 |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201113 Address after: Door 19, no.23-1, Hongze Road, Yuhong District, Shenyang City, Liaoning Province Patentee after: Shenyang Hetai Technology Co.,Ltd. Address before: 110044 Liaoning Province, Shenyang City District wanghuanan Street No. 21 Patentee before: SHENYANG University |