CN101458735A - 一种具有高升阻比的翼型 - Google Patents
一种具有高升阻比的翼型 Download PDFInfo
- Publication number
- CN101458735A CN101458735A CNA2008102370169A CN200810237016A CN101458735A CN 101458735 A CN101458735 A CN 101458735A CN A2008102370169 A CNA2008102370169 A CN A2008102370169A CN 200810237016 A CN200810237016 A CN 200810237016A CN 101458735 A CN101458735 A CN 101458735A
- Authority
- CN
- China
- Prior art keywords
- airfoil
- profile
- drag ratio
- high lift
- functional
- 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.)
- Pending
Links
- 238000013461 design Methods 0.000 claims abstract description 23
- 230000010354 integration Effects 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 3
- 238000010248 power generation Methods 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 3
- 238000012938 design process Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
Images
Classifications
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Wind Motors (AREA)
Abstract
本发明涉及基于泛函的通用翼型型线形状设计的翼型,通过建立翼型型线的泛函集成方程,然后对该集成方程建立优化设计模型,根据所要求的目标函数选择不同的方程项数和系数来完成目标翼型的设计。本发明提供的基于泛函的通用翼型型线的高升阻比翼型与传统的风力机翼型相比,在正常工作功角范围内,不仅升力系数高,而且具有很高的升阻比,失速较晚。它的构成型线不但能够发挥通用翼型型线的优势,其构成的新型型线还能满足单一型线不能达到的技术指标,从而提高了风力机的输出功率,降低了风力发电成本。
Description
技术领域
本发明涉及风力机叶片翼型的设计,具体涉及一种具有高升阻比的翼型。它除了适应于风力机叶片外,也于适用于飞机机翼的设计制造。
背景技术
风力机叶片翼型设计理论是决定风力机功率特性和载荷特性的根本因素,一直是各国学者研究的热点所在。目前叶片翼型的研究和设计都是基于特定的几何轮廓曲面,来研究其空气动力学特性,并进行参数的设计和优化。风力机专门翼型是在航空翼型的基础上发展起来的,已经研究出的航空叶片翼型有美国航空航天局的NACA230系列、NACA44系列、NACA63-2系列、NACALS系列和美国贝尔公司的FX系列等;专门的风力机翼型有美国可再生能源实验室的NREL S系列、瑞典FFA-W系列、丹麦Ris系列和荷兰DU系列等。
翼型的性能对风力机的气动性能具有决定性的影响,高性能翼型的研究是风力机发展的一项基础性研究。高升力、低阻力的翼型一直是翼型设计中所追求的目标。对于翼型的气动外形设计目前通常采用反设计方法,这种方法是根据给定希望达到的气动特性,压力分布以及初始的基本翼型,通过几何和流动控制方程,逐步逼近给定的气动特性。这种设计方法不够直接,设计过程也比较耗时,而且往往逼近的结果不够准确。因此,直接从翼型的形状表达出发,寻找一种更为直接高效的设计方法,就显得尤为必要。
此外,由于特定的几何曲线数学模型一经确定,其固有的几何特性和数学特性无法变更,因而性能受到根本制约。对于翼型的设计如果只基于给定的翼型形状,进行参数的优化、改进,并不能对表征叶片翼型型线本质特征的根本因素—翼型型线的形函数本身进行优选,因而难以在叶片翼型型线的研究和创新中取得根本的突破。
因此,研究风力机叶片翼型型线应具有怎样的函数特性,该函数的数学模型如何表达,该数学模型表征的几何曲线具有怎样的空气动力学特性,如何对该几何曲线进行形状优化等根本问题,对于建立风力机叶片翼型型线的完整理论,并根据风电场的实际情况构造出各种具有最佳性能的叶片翼型型线,具有重大的理论意义和工业应用价值。
综上所述,本专利通过研究翼型的通用特性,建立其泛函集成方程,进而提出了一种全新的、高效的翼型优化设计方法。翼型采用泛函集成方程来表达,其几何形状及解析特性可通过级数系数的调整和优化加以控制;形状优化的泛函极值问题也可以通过对可变项数级数系数的参数优化来逐次逼近。进而实现各种翼型形状的构造。
发明内容
本发明的目的就是提供了一种基于泛函的通用翼型型线形状设计的翼型。
本发明以上目的是通过以下技术方案来实现的,首先建立翼型型线的泛函集成方程,然后对该集成方程建立优化设计模型,根据所要求的目标函数选择不同的方程项数和系数来完成目标翼型的设计。
建立笛卡儿坐标系下z平面上的一个圆,通过改变圆心的位置,利用儒科夫斯基变换式:
ζ=f(z)=z+a2/z (1)
就能够变换成ζ平面上的一个翼型,其中,a为1/4翼型弦长。
此时,ζ平面上的翼型在笛卡儿坐标系下就可以表示为:
其中,r为翼型的矢径长度。
为了表征翼型的通用性,r其实是一个关于θ的函数,表示为:
通过选取不同的就可以变换出各种不同性能的翼型。
k=1,2,3,...n (4)
对风力机翼型而言,其升力与阻力比是最重要的性能指标,对提高风轮效率和对整个风力机组工作性能的影响关系极大。在翼型设计过程中,一般期望得到较大的升力同时阻力又保持很小,所以将翼型的升阻比作为目标函数,建立广义翼型型线泛函极值模型:
f(x)=max(Cl/Cd) (5)
其中Cl位翼型的升力系数,Cd为翼型的阻力系数。
X=(x1,x2,x3,x4,x5,x6) (6)
式中变量的变化范围根据翼型的形状来确定。当各个变量的值超过一定范围时,该表征式将不再具有翼型的形状特征,所以对变量的范围进行了约束:
Xmin≤X≤Xmax。(7)
除了对形状的要求以外,其它的约束主要针对翼型的结构和强度。通常情况下,叶片主要产生功率区翼型的最大相对厚度为20%左右,选择翼型的最大相对厚度为X,即施加约束:
t/c=X (8)
其中t为翼型的最大厚度,c为翼型的弦长。
翼型的结构和强度特性除了和最大相对厚度有关以外,还与翼型最大相对厚度所处的弦长位置有关,一般翼型的转矩中心在其弦长的1/4处左右,所以翼型的最大相对厚度在弦长所处的位置控制在
0.2≤x/c≤0.3。 (9)
基于以上的翼型型线优化设计模型,以型线集成方程系数为变量,通过改变、舍取级数项数,从而使型线得到优化,设计出具有高升阻比的翼型。
有效结果:本发明提供的基于泛函的通用翼型型线的高升阻比翼型与传统的风力机翼型相比,在正常工作功角范围内,不仅升力系数高,而且具有很高的升阻比,失速较晚。它的构成型线不但能够发挥通用翼型型线的优势,其构成的新型型线还能满足单一型线不能达到的技术指标,从而提高了风力机的输出功率,降低了风力发电成本,必将促进风电产业的大力发展。
附图说明
下面结合附图给出的一个非限定性的实施例对本发明作进一步说明:
图1——相对厚度18%的翼型形状;
图2——相对厚度18%的翼型升力和功角关系图;
图3——相对厚度18%的翼型升阻比关系图;
图1为最大相对厚度为18%的翼型,该翼型最大相对厚度在弦长的27%处。
具体实施方式
在实施例给出的翼型的设计过程中,根据TAYLOR级数思想和级数函数收敛原则取前6项,选择一定的相对厚度来进行翼型形状设计。这里给出了最大相对厚度为18%的翼型,表1给出了这种翼型泛函集成方程的系数,将这些系数代入到集成方程中就能得到翼型的形状,该翼型的最大相对厚度为18%,最大相对厚度在弦长的27%处。见图1。
表1 翼型泛函集成方程系数
最大相对厚度 | X1 | x2 | x3 | x4 | x5 | x6 |
18% | 0.1132 | 0.1328 | 0.0072 | 0.0261 | 0.0001 | 0.0189 |
采用XFOIL软件对于本发明提出的这种翼型的气动性能进行了计算分析。雷诺数Re=1.6×106时,最大相对厚度为18%的本发明翼型的升阻比见表2,当功角为5°时,拥有最大升阻比161.8465,而且该翼型在风力机叶片正常工作功角0°~10°内均具有很大的升阻比。
表2雷诺数为Re=1.6×106的最大相对厚度为18%翼型的升阻比
角度(°) | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
升阻比 | 126.8564 | 138.9669 | 148.8447 | 156.1282 | 160.9059 | 161.8465 | 157.6797 | 147.1154 |
角度(°) | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 |
升阻比 | 138.2119 | 119.0821 | 98.9927 | 80.6389 | 63.9088 | 50.6272 | 40.6724 | 33.3490 |
Claims (2)
1、一种具有高升阻比的翼型,由具有翼型特性的型线平滑连接构成,其特征在于,具有翼型特性的型线按一下方式构建:
建立在笛卡儿坐标系下z平面上的一个圆,通过改变圆心的位置,利用儒科夫斯基变换式:
其中,r为翼型的矢径长度,是一个关于θ的函数,表示为:
翼型族的形状由式(4)表示的通用翼型型线来控制,上述泛函表征式集成了所有具有翼型特性的光滑曲线;集成型线系数参数通过利用Taylor级数及其收敛原则,以该系数为变量,通过改变、舍取级数项数,最终优化设计得到了具有高升阻比的集成翼型系列。
2、根据权利要求1所述的具有高升阻比的翼型,其特征是:将翼型的升阻比作为目标函数,建立广义翼型型线泛函极值模型:
f(x)=max(Cl/Cd) (5)
其中Cl位翼型的升力系数,Cd为翼型的阻力系数;
X=(x1,x2,x3,x4,x5,x6) (6)
式中变量的变化范围根据翼型的形状来确定,当各个变量的值超过一定范围时,该表征式将不再具有翼型的形状特征,所以对变量的范围进行了约束:
Xmin≤X≤Xmax (7)
其它的约束主要针对翼型的结构和强度,通常情况下,叶片主要产生功率区翼型的最大相对厚度为20%左右,选择翼型的最大相对厚度为X,即施加约束:
t/c=X (8)
其中t为翼型的最大厚度,c为翼型的弦长;
翼型的转矩中心在其弦长的1/4处左右,所以翼型的最大相对厚度在弦长所处的位置控制在
0.2≤x/c≤0.3。 (9)
基于以上的翼型型线优化设计模型,以型线集成方程系数为变量,通过改变、舍取级数项数,从而使型线得到优化,设计出具有高升阻比的翼型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2008102370169A CN101458735A (zh) | 2008-12-31 | 2008-12-31 | 一种具有高升阻比的翼型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2008102370169A CN101458735A (zh) | 2008-12-31 | 2008-12-31 | 一种具有高升阻比的翼型 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101458735A true CN101458735A (zh) | 2009-06-17 |
Family
ID=40769592
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2008102370169A Pending CN101458735A (zh) | 2008-12-31 | 2008-12-31 | 一种具有高升阻比的翼型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101458735A (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813074A (zh) * | 2010-03-25 | 2010-08-25 | 张玉福 | 简化的环形大气动力机 |
CN101615216B (zh) * | 2009-07-28 | 2011-01-05 | 重庆大学 | 对翼型尾缘进行流线形增厚的方法 |
CN101923584B (zh) * | 2009-10-30 | 2012-08-15 | 重庆大学 | 风力机专用翼型设计方法及风力机专用翼型 |
CN103136422A (zh) * | 2013-01-11 | 2013-06-05 | 重庆大学 | 翼型集成与b样条结合的中等厚度翼型设计方法 |
CN104317992A (zh) * | 2014-10-11 | 2015-01-28 | 重庆工商大学 | 风力机翼型正设计方法及风力机翼型族 |
CN104354875A (zh) * | 2014-10-31 | 2015-02-18 | 中国航天空气动力技术研究院 | 一种考虑运动轨迹的增升装置优化设计方法 |
CN104401504A (zh) * | 2014-11-19 | 2015-03-11 | 中国地质大学(武汉) | 一种固定翼航测型无人机设计方法 |
CN104612892A (zh) * | 2014-12-30 | 2015-05-13 | 中国科学院工程热物理研究所 | 一种风力机翼型的多学科优化设计方法 |
CN104834772A (zh) * | 2015-04-22 | 2015-08-12 | 复旦大学 | 基于人工神经网络的飞机翼型/机翼反设计方法 |
CN105822510A (zh) * | 2016-04-12 | 2016-08-03 | 南京航空航天大学 | 一种叶尖系列化分段叶片及其设计方法 |
CN108331708A (zh) * | 2018-01-25 | 2018-07-27 | 华北电力大学 | 风力发电叶片翼型线的确定方法 |
CN111859801A (zh) * | 2020-07-16 | 2020-10-30 | 湖北工业大学 | 失速型风力机翼型的设计方法 |
-
2008
- 2008-12-31 CN CNA2008102370169A patent/CN101458735A/zh active Pending
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615216B (zh) * | 2009-07-28 | 2011-01-05 | 重庆大学 | 对翼型尾缘进行流线形增厚的方法 |
CN101923584B (zh) * | 2009-10-30 | 2012-08-15 | 重庆大学 | 风力机专用翼型设计方法及风力机专用翼型 |
CN101813074A (zh) * | 2010-03-25 | 2010-08-25 | 张玉福 | 简化的环形大气动力机 |
CN103136422B (zh) * | 2013-01-11 | 2015-12-02 | 重庆大学 | 翼型集成与b样条结合的中等厚度翼型设计方法 |
CN103136422A (zh) * | 2013-01-11 | 2013-06-05 | 重庆大学 | 翼型集成与b样条结合的中等厚度翼型设计方法 |
CN104317992A (zh) * | 2014-10-11 | 2015-01-28 | 重庆工商大学 | 风力机翼型正设计方法及风力机翼型族 |
CN104354875A (zh) * | 2014-10-31 | 2015-02-18 | 中国航天空气动力技术研究院 | 一种考虑运动轨迹的增升装置优化设计方法 |
CN104354875B (zh) * | 2014-10-31 | 2016-08-24 | 中国航天空气动力技术研究院 | 一种考虑运动轨迹的增升装置优化设计方法 |
CN104401504A (zh) * | 2014-11-19 | 2015-03-11 | 中国地质大学(武汉) | 一种固定翼航测型无人机设计方法 |
CN104612892B (zh) * | 2014-12-30 | 2017-03-29 | 中国科学院工程热物理研究所 | 一种风力机翼型的多学科优化设计方法 |
CN104612892A (zh) * | 2014-12-30 | 2015-05-13 | 中国科学院工程热物理研究所 | 一种风力机翼型的多学科优化设计方法 |
CN104834772A (zh) * | 2015-04-22 | 2015-08-12 | 复旦大学 | 基于人工神经网络的飞机翼型/机翼反设计方法 |
CN104834772B (zh) * | 2015-04-22 | 2018-01-12 | 复旦大学 | 基于人工神经网络的飞机翼型/机翼反设计方法 |
CN105822510A (zh) * | 2016-04-12 | 2016-08-03 | 南京航空航天大学 | 一种叶尖系列化分段叶片及其设计方法 |
CN108331708A (zh) * | 2018-01-25 | 2018-07-27 | 华北电力大学 | 风力发电叶片翼型线的确定方法 |
CN111859801A (zh) * | 2020-07-16 | 2020-10-30 | 湖北工业大学 | 失速型风力机翼型的设计方法 |
CN111859801B (zh) * | 2020-07-16 | 2022-07-19 | 湖北工业大学 | 失速型风力机翼型的设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101458735A (zh) | 一种具有高升阻比的翼型 | |
CN103277245B (zh) | 一族大厚度钝尾缘风力机翼型及其设计方法 | |
CN104612892B (zh) | 一种风力机翼型的多学科优化设计方法 | |
CN107559143B (zh) | 一种大型风力机尾缘襟翼结构参数寻优及多目标襟翼优化控制方法 | |
CN101923584B (zh) | 风力机专用翼型设计方法及风力机专用翼型 | |
CN103136422B (zh) | 翼型集成与b样条结合的中等厚度翼型设计方法 | |
CN112507471B (zh) | 一种大攻角范围条件下垂直轴风力机翼型设计方法 | |
CN106777567A (zh) | 一种风力机钝尾缘翼型优化设计方法 | |
CN109145506B (zh) | 一种高气动性能低噪声水平风力机外侧翼型的优化设计方法 | |
CN104863799A (zh) | 一种利用贝塞尔函数曲线的风力机翼型设计方法 | |
CN106919749A (zh) | 一种低噪声风力机叶片设计方法及低噪声风力机叶片 | |
CN203374428U (zh) | 一族大厚度钝尾缘风力机翼型 | |
CN113742861A (zh) | 一种适用于风力发电机风洞试验的叶片模型优化设计方法 | |
CN101615216A (zh) | 对翼型尾缘进行流线形增厚的方法 | |
CN105787212A (zh) | 一种抗气动弹性变形的风力机翼型设计方法 | |
CN103711655B (zh) | 一种大厚度钝尾缘风力机叶片 | |
CN102400847A (zh) | 一种风力机叶片翼型 | |
CN102094767A (zh) | 一族用于兆瓦级风力机叶片的翼型 | |
CN104317992A (zh) | 风力机翼型正设计方法及风力机翼型族 | |
CN105868470A (zh) | 一种风力机翼型与叶片外形参数一体化设计方法 | |
CN101898644A (zh) | 用来设计翼型的系统与方法 | |
CN109783912B (zh) | 一种升力型垂直轴风力发电机叶片翼型优化设计方法 | |
Grasso | Design of a family of new advanced airfoils for low wind class turbines | |
CN106741924A (zh) | 一种全工况条件下的高升阻比低力矩特性9%厚度旋翼翼型 | |
CN106599381A (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 | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Open date: 20090617 |